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TURBULENCE CHARACTERISTICS OF THE HYDRAULIC JUMP 


Hunter Rouse,! M. ASCE, T. T. Siao,2 and S. Nagaratnam3 
(Proc. Paper 1528) 


ABSTRACT 


Hot-wire measurements of the turbulence in an air-flow model of the 
tydraulic jump are described for Froude numbers of 2, 4, and 6. Results 
are analyzed and interpreted in the light of the momentum and energy 
integrals. 


SYNOPSIS 


Because the fluid discontinuities produced by entrapped air interfere with 
he measurement of turbulence in the hydraulic jump itself, the flow pattern 
4f the jump has been simulated in an air duct and the hot-wire anemometer 
ised to determine the corresponding pattern of turbulence for Froude num- 
ers of 2, 4, and 6. In an initial analytical section of this paper the modeling 
method is justified and the differential and integral forms of the momentum 
und energy equations pertinent to the investigation are explained. The ex- 
erimental apparatus and techniques are next described and the basic data 
resented in diagram form. The measured characteristics of the mean flow 
ind of the turbulence are then correlated in accordance with the foregoing 
squations, and the process of energy conversion from the mean flow through 
urbulence to heat is followed in detail. Various aspects of the hydraulic 
ump long subject to conjecture or misunderstanding are thereby clarified. 


Jote: Discussion open until July 1, 1958. A postponement of this closing date can be 
obtained by writing to the ASCE Manager of Technical Publications. Paper 1528 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 1, February, 1958. 
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INTRODUCTION 


Scattered through the twelve discussions of the Bakhmeteff-Matzke paper, 
“The Hydraulic Jump in Terms of Dynamic Similarity,” (1) are three sorts of 
information pertinent to the present investigation: first, a well-annotated ac- 
count of the gradual development of knowledge on this subject; second, a 
somewhat contradictory mixture of factual matter, hypothesis, and error; 
third, and by inference only, a guide to further research. At the time the 
paper and discussions were published (1935), the tools for such research had 
only recently come into use in other fields, and it is hardly to be expected 
that they would already have been applied in hydraulics. In fact, the investi- 
gation herein described had the purpose of demonstrating the efficacy of the 
new approach as well as providing greater insight into the phenomenon under 
consideration. 

As was evident from the aforementioned paper and discussions, past ad- 
vancements had been concerned with the vertical elements -i.e., depths and 
heads - of the mean flow pattern as correlated by the one-dimensional forms 
of the continuity, momentum, and energy relationships. Some attention had 
been given to the empirical evaluation of the longitudinal elements as well, 
and to the expression of the functional relationships in terms of what is now 
called the Froude number. Concern with the accompanying pattern of turbu- 
lence, however, was almost entirely speculative. For lack of even an ap- 
proximate evaluation of the turbulence characteristics, their importance to 
the phenomenon was both over- and underestimated. The first of the present: 
writers, for example, claimed that the computed loss in head at the jump 
represented merely a transfer of energy from the mean motion to the turbu- 
lence, and that the total head would hence be found to remain practically 
constant through the jump if the kinetic energy of the turbulence could be 
measured and included as velocity head. Others, on the contrary, under- 
estimated the importance of eddy generation along the surface of dis- 
continuity between the secondary flow of the roller and the primary flow 
beneath - some even going so far as to claim that the roller played no 
essential part in the phenomenon. 

The difficulty lay, as is now apparent, in the general lack of clarity which: 
then prevailed in nearly all matters involving fluid turbulence. Not only did 
existing turbulence theories seem to have little connection with hydraulic 
reality, but the only practicable means of measuring turbulence - the hot- 
wire anemometer - was restricted to use in air. During the intervening 
years, progress has been made in the interpretation - if not the quantitative 
prediction - of the transformation from energy of the mean motion through 
energy of turbulence to that of heat, and the hot-wire technique itself has 
been extended to use in water. Means would therefore seem to be at hand 
for the measurement of the primary and secondary patterns of flow in the 
hydraulic jump, and at least the evaluation of the various forms at represen- 
tative sections. Unfortunately, the jump will probably be one of the last 
hydraulic phenomena to prove susceptible to exploration with the hot-wire 
instrument, because of the presence of countless fluid discontinuities - 
bubbles of entrained air - in the region of greatest interest. 

To overcome this rather formidable obstacle (yet without taking recourse: 
to the extremely tedious process of statistically analyzing motion pictures 02 
suspended particles), it was resolved at the Iowa Institute of Hydraulic Re- i 
search to simulate the flow pattern of the hydraulic jump in an air duct, and h 
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hen use the hot-wire anemometer for turbulence surveys just as in previous 
nstitute studies of other diffusion phenomena. The analog will seem at first 
slance a rather artificial one, to say the least. However, it must be realized 
hat few model flows portray perfect dynamic similarity with their proto- 
ypes, and that for all practical purposes only the pertinent part of the 
Ihenomenon need be closely approximated. In the present instance it is 

lot the gravitational action that is to be evaluated but the viscous. If the 
nean flow patterns are geometrically similar, and if the changes in energy 
ire comparable, then it would seem safe to assume that the patterns of tur- 
yulence are also similar - and it is specifically these with which the present 
nvestigation is concerned. 


Analytical Bases of Experimentation 


Air-Water Correlation 


With reference to the schematic representation of Fig. 1 and the list of 
symbols at the end of the paper, application of the simplified continuity, 
nomentum, and energy relationships for one-dimensional flow leads to the 
ollowing well-known equations for the hydraulic jump in homogeneous liquid 
lowing through a horizontal, resistance-free channel of rectangular cross 
section: 


find (/reeuin) (1) 


a, Vv 
Rival dhoobd, tof peerdioy wake 2 
sae maak | el (2) 


30th the depth ratio dg/dj and the relative loss of total head Hy, /a1 are seen 
0 be explicit functions of the Froude number of the approaching flow 


fy = Uj/Vgd1 , as are also such longitudinal characteristics as the relative 
ength of the jump Lj/d1 and the shape of the flow profile. 


Now it is not too difficult to visualize the application to the flow of an up- 
ver boundary having the same profile as the mean free surface. Like the 
spillway profile that is shaped according to the lower surface of the cor- 
“esponding weir nappe, the presence of the solid boundary may be expected 
o have some effect upon the flow. In the case of the spillway, the predomi- 
lant effect is the addition of shear. In the case of the jump, it is rather the 
‘limination of both the freedom of the surface to fluctuate and the presence 
f entrained air. However, measurements of the pressure head along the 
loor of a flume under a jump yield values which differ from the local depths 
y no more than the experimental error, indicating at once that the pressure 
listribution can be assumed to be hydrostatic and that the quantity of air in 
suspension is actually insufficient to change appreciably either the density 
wr the specific weight. As for the surface fluctuations, measurement of the 
‘esulting wave amplitudes just downstream and approximation therefrom of 
he resulting wave energy has shown this to be less than a percent of the 
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energy of the flow itself. It would thus seem reasonable to assume that the 
presence of the upper boundary would have at the most only a secondary ef- 
fect upon the flow pattern. 

If the foregoing premise is accepted, then it must also be granted ( in ac- 
cordance with the general mechanics of confined flow) that the same pattern 
would obtain between the upper and lower boundaries at other rates of flow, 
if only the Reynolds numbers did not differ greatly. However, no longer 
would the mean pressure intensity along the upper boundary be atmospheric ~ 
at all points, for it must change in accordance with the Bernoulli equation. 
On the other hand, if the linear quantities in the foregoing equations are 
considered to be the piezometric heads h (i.e., the sum of pressure head and 
elevation, which is presumably constant over any vertical section), the 
following relationships must continue to prevail between the enclosed model 
and its free-surface prototype. 


Wee ny -_ d2- a 
2/2 U,?/2g (3) 


= he-h, ae d,* = Ay 
|- Ug dt Uz72% ie 


Again in accordance with the principles of fluid mechanics, it now makes 
no difference whether the fluid is a liquid or a gas - just as long as the 
Reynolds numbers are still comparable and the velocities of the gas are not 
so high (i.e., two or three hundred feet per second) as to cause compressi- 
bility effects. Hence the same patterns of mean flow and turbulence should 
be obtainable through use of air of density p as the experimental fluid. 
Moreover, since the elevation term could then be neglected in the Bernoulli 


and related equations, the foregoing model-prototype expressions would re- 
duce to the following pressure forms: 


Pe= Pic pe plana 
pU?/2 be /24 (5) 


pU,?/2 a, U,2/2g (6)) 


Although the Froude number cannot be expressed by replacing dj with 
either hy or pj/Y (since the pressure load on the system is no longer 
determined by free-surface conditions), it will be seen that the right-hand 
sides of Eqs. (3) and (5) are still completely governed by the prototype 
magnitude of F; - and hence the left-hand sides as well. Moreover, for 


any Froude number the pressure-distribution terms (p - P1)/(Pg - py) and 
(h - hy)/(hg - hy) should be numerically equal to the profile dimension 
(d - dy)/(dg - dy) at homologous sections. To the extent that these equalities 
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Fig. 1. Definition sketch of hydraulic jump. 


Fig. 3. Location of stream lines from velocity distribution. 
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are actually realized in prototype and model, the flow patterns themselves 
may be considered comparable. 


Momentum Relationship 


To obtain a more precise form of the momentum equation than that cus- 
tomarily used in one-dimensional hydraulics, recourse must be had to the 
differential equations of motion! as adapted to turbulent flow by Reynolds 
from the form derived for laminar flow by Stokes, and as presented in 
standard texts on hydrodynamics. (2) This adaptation involved, in brief, the 
replacement of every instantaneous velocity magnitude by the sum of its 
mean value and its deviation therefrom, and then the retention of only those 
terms having a temporal average other than zero. The local vector magni- 
tude V is thus considered to have in general the mean components u, v, and 
w inthe directions x, y, and z, respectively, with the corresponding instan- 
taneous deviations u', v', and w'. For the direction x the Reynolds 
equation has the typical form (bars denoting temporal means) 


r = mA? 1 > |! 
G24 4 yok. wol 4. yt a aye ee 
us Y Z x oy dz 


ay UD ae “(2 du pats 
a p Ox? Oy? dz? 


(7) 


Herein yw is the dynamic viscosity and X_the x-component of the body force 
(such as the attraction of gravity) per unit mass. 

This is fundamentally an equation of acceleration, for the terms to the left 
of the equality sign represent the rates of change of the velocity components 
of the mean motion and of the turbulence, and the terms to the right represent! 
the pressure, body, and viscous forces per unit mass by which the velocity — 
changes are produced. Through introduction of the two equations of continuity; 


ou, N tw _ (8) 


and 


Cur oe Ee 
ax Gy Oz (9) 


1. Space limitations prevent the inclusion in full of the derivations outlined 
in the following pages. However, they are readily available in the litera- 
ture, and they can probably be checked by those with advanced training 
even without reference to the papers cited. What is important to the 
general reader is the implicit demonstration that (a) the customary 
jump relationships are one-dimensional approximations (b) a more re- 
fined analysis of the role of turbulence must stem from ‘the basic equa- 
tions of motion, (c) the derived expressions must in turn be simplified, 


but only in accordance with experimental evidence i 
and (d = 
sults are broadly significant. ; (d) the final re 
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he acceleration terms can be written in the alternative form 


(ae) , tev", uw) out , Ouv', du'w' 
ex dy vz ex Oy Jz 


which is used in a subsequent derivation. 

Since there is an equation of acceleration for each of the three coordinate 
lirections, a considerable saving in space, without sacrifice of clarity, is 
gained by the use of tensor notation, in which the quantities u, v, and w are 
represented by u,; or Uj, and x, y, and z by x; or xj. In shorthand form, 
the indices i and j denote in turn each of the three coordinate directions. 
The foregoing equation and its counterparts for the other two directions are 
hus represented at one and the same time by 


a.) , wis _ 10 (Ky, weH 
se iitaet aed a AR A te OOM: (10) 


wherein i refers to the direction characterizing each equation and j to the 
jirection of each independent term within it. Repetition of indices thus de- 
1otes a summation of terms. 

Any equation of acceleration may - like the Newtonian equation itself - 
also be regarded as an equation of momentum, and it will be seen that multi- 
lication of both sides of the foregoing relationship by the mass density will 
yive each term the dimension of rate of change of momentum per unit volume 
Ir impulse per unit volume per unit time. The physical significance of this 
2xpression becomes the greater as it is integrated over a given space through 
which the fluid moves. By the Gaussian rule of the calculus relating volume . 
mtegrals and surface integrals 


eon ie ons 45 (11) 


he result can be put in the following general form: 
OX. Na ee ORS w 
a Yay 5 aS + [om = ds 
o% 45 + /oXd¥% + wi ae 
ae ie wf Ox (12) 


Jerein S denotes the surface of the region over which the integration is per- 
ormed, n the outward normal to the surface, and ¥the enclosed volume. 
[The terms at the left (each representing nine different quantities) embody 

he net flux of momentum (three components in each of three directions) of 
he mean flow and of the turbulence out of the region in question - i.e., the 
lifference between the rates of efflux and influx. The first term at the right 
lenotes the three components of the mean normal force exerted externally on 
he surface of the region; the second, the three components of the weight of 
he fluid contained within the region; and the third, the three components of 
he méan tangential force exerted on the surface. 


a 


7 
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For the case in which the hydraulic jump is formed on a horizontal bed of 
great width, the pertinent momentum equation is that for the direction x. If — 
it is assumed that the turbulence is negligible at the initial section, that the 
pressure distribution is hydrostatic at any arbitrary section, and that both the 
viscous and the turbulent stresses are negligible over the free surface, the 
equation reduces to 


d d,  — 
[peau - pur dy + [euay 
fy 0 0 


_ wd? _ yd? , f* fu) g | 
= Wh Es false) dy (13) 


For the limiting case that x = L, and d = do, this expression will be seen to. 


differ in the following ways from the simpler version normally used in study- 
ing the jump. First, the change in the mean momentum flux, expressed as the’ 
difference between the first two integrals, per-force includes the effect of 
local departure from the average velocity over the vertical section. Second, 
the change in the momentum flux of the turbulence is represented by the next 
integral, which embodies the sole effect of turbulence in the force balance. 
Finally, the effect of variable bed shear Tg = yp (0 0/8 ¥)y-0 over the length: 


of the jump is included in the integral at the right. It goes without saying 
that application of the equation requires the determination of each integral 
from known distributions of the variables in question - an evident reason for 
the more customary use of the simplified relationship. Only through the 
quantitative evaluation of the integral terms, however, can the magnitude of 
the error arising from the neglect of the distribution details be appraised. 


Equations of Energy 


An essential characteristic of the Reynolds equations of motion, and hence: 
of the momentum equations derived from them, is the absence of all viscous 
stresses except those involving the mean velocities. However, although the 
viscous terms containing the velocity fluctuations are invariably eliminated 
by the averaging process, there remain certain inertial terms consisting of 
fluctuation products which can conveniently be considered to have the nature 
of a stress and hence to take their place. These are the momentum-transpor 


terms of the secondary motion, which may be combined with the terms for 
mean viscous stress to yield 

Se =. | OU Oe ee 

0 als i Ox; f anels 


(14) 


The momentum-transport terms pu! iU' j are known as Reynolds stresses. 


If, before averaging, each of the equations of motion is multiplied by the 
corresponding component of the instantaneous velocity, and the three equa- 


tions are then added, there will result a differential equation of work and — 
energy(3) of the form 
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(eV? , pV? ) MeV 22 iy lev , 2d, pity) 
J dx; Z. af 2 t - x; 0x; 


(15) 


= z = Poa 2 
Be hy Dees OP 5 eh a Se FN Or fn ee 
‘ : Ox, ox; 
Inspection of this equation will reveal not only that many of the fluctuating 
stresses which disappeared from the momentum equations during the averag- 
ing process are now retained, but also that there is a fairly consistent parallel 
between the terms for the mean flow and those for the turbulence. In fact, it is 
instructive to segregate the terms of the two categories in the following work- 
energy relationships for the mean flow and for the turbulence (which are also 


directly derivable by combining the pertinent equations and velocity compon- 
ents): 


A: a(oV 7/2) ne duu, ive Mine op = = — 07u 
2 Ox; vege BAe 4 ast etek eerie Ox, Ox (16) 
= ov'/2) , 2) cee Oe mney arya 9 BRE (17) 
mee ie : Ox, ox 


As was done with the momentum form of the equations of motion, it is now 
in order to integrate the differential work-energy equations over a given re- 
gion of space. After the appropriate volume integrals have been changed by 
the Gaussian rule to surface integrals, the work-energy equation for the mean 
flow assumes the form 


ee Ox; Ox, 4a _ ; a 
Meee a + [eS nite EEO G ae 


. Ba S00 + [ ORD, av | ale +S a OAS -[n( oe (18) 


The successive terms have the following significance: That at the extreme 
left represents the net flux of kinetic energy of the mean motion out of the re- 
gion in question (i.e., the rate of increase in kinetic energy as the fluid passes 
through the region). The second and third terms, while also representing en- 
ergy transport, can best be explained as the rates at which work is done by the 
Reynolds stresses over the surface and throughout the interior of the region, 
respectively. The first and second terms at the right of the equality sign are 
the rates at which work is done by the external pressures and the body forces. 
The third and fourth are evidently the rates at which work is done by the vis- 
cous stresses of the mean flow over the surface and throughout the interior of 
the region, respectively. 

Now the work done externally by the viscous stresses (the third term at 
the right) is wholly conservative (since, as is evident from the Gaussian 
transformation, it is derivable from a potential or space derivative), whereas 
that done internally (the last term) is wholly dissipative. Because of the 
analogy between the turbulent and viscous stresses, one might seek a further 
parallel between the types of work done in each instance. In fact, whereas 
the second term at the left, like the third at the right, is conservative, the 
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third at the left can be considered dissipative like the last at the right if one 
assumes that energy transferred from the mean flow to the turbulence can 
never be recovered. The third term at the left, in other words, represents — 
the rate at which turbulence is produced at the expense of the mean flow. 
(The integral itself is inherently negative, so that the negative sign before it 
denotes a positive rate of production.) 

The corresponding work-energy equation for the secondary motionis as 
follows: 


&, Jaq 


Gee ae yi2 y os. 
Ve a ae fae uw, i do + | oe, OU gay 
a lees on 4 OX; 


__ fag Oh a, Duy Oe WTA VIE 
= fae dS + [a(S tpt an Ge [ul + 2s) Gu dX (19) 
r vy fay J 


The first two terms at the left, by analogy to the first in the equation for the © 
primary motion, represent the net flux of kinetic energy of turbulence out of © 
the region by convection and diffusion (i.e., by the mean flow and by the 
turbulence), respectively. The third is again the rate of production of 
turbulence, now a negative quantity so far as the work accomplished by the 
turbulence is concerned. The first term on the right of the equality sign is 
the rate at which work is done externally on the surface of the region by the 
fluctuating pressures. The last two terms represent the rates at which work 
is done by the viscous stresses of the turbulence over the surface of the re- 
gion and throughout the interior, respectively, the former being conservative 
and the latter dissipative. : 

On the previous assumption that the initial section is free from turbulence, 
that both the mean and turbulent stresses can be neglected over the free sur- 
face, and that the pressure distribution is hydrostatic throughout, reduction 
of the work-energy equation for the mean motion to the case of the hydraulic 
jump on a level bed of great width yields the following result: 


d 72 dV nd — 3 soa 
fe wage +f wut VU, 
<Q J, 25 e g 4 


ae Pie fou Oe, eae a 
[PUES )s 22 aude = gd, 9d 


[oR Blo (PA iS Be al 


For the limiting case that x = L and d = do, the first two terms on each 


side of the equation correspond to those of the usual one -dimensional open- 
channel relationship, with due account taken of the variation in velocity over ; 
the cross section. The third terms represent work done on the end section _ 
by the turbulent and viscous stresses. The last terms combine to form the _ 
quantity Hj q, since they indicate the rate at which turbulence is produced | 
by the mean motion and the rate at which energy is dissipated by the mean 


ASCE HYDRAULIC JUMP 1528-11 


viscous stresses. 


Reduction of the energy equation for the turbulent motion to the same con- 
ditions yields the following result: 


Sore rd Ne as a ‘dx a 2 
ig a / 1) , | & fo ENE 72) Ou Y ony Mid 
fy u’ dy ay ar Ge US on ay {] Ke eu) dys (21) 


Herein terms appear only for the arbitrary section, since it was assumed that 
no turbulence existed at the initial section. The proportionality factor K 
permits what are actually 15 terms to be written as one; its use is further 
justified by the fact that (a) the dissipation is most intense in the smallest 
eddies, (b) the smallest eddies tend to be isotropic, and (c) in isotropic tur- 
bulence this factor has the constant magnitude of 15.(4) Only if the cumula- 
tive production (the third term on the left) and the cumulative dissipation 

(the third term on the right) are equal will the turbulence at the final section 
me.,for d=do at x = Lj) be negligible. The extent of the actual inequality 
of these terms can evidently be ascertained only through detailed measure- 
ment and analysis of the mean and fluctuating patterns of flow. 


Experimental Equipment and Procedure 


To simulate the mean and fluctuating patterns of the hydraulic jump in air, 
a special duct having a length of 9 feet, a width of 1 foot, and a maximum 
depth of 1 foot, was built of plywood and plexiglass. Both the front wall and 
the curved surface were made transparent to simplify the visual setting of 
instruments, and for added convenience in supporting them the plane surface 
corresponding to the bed of the channel was placed on top (see Fig. 2). The 
form of the curved surface was controlled by successive pairs of adjusting 
screws at 6-inch intervals over the entire length of the duct. Both the plane 
and the curved surfaces were provided with piezometer holes along the 
centerline, and a grid with a 0.02-foot mesh was ruled on the rear wall as a 
suide in controlling the surface profile. At the downstream end the duct 
underwent a smooth transition to a 12-inch circular section connected to the 
intake of a 5-horsepower centrifugal blower, which discharged into the sur- 
rounding space through a baffled diffuser. The intake end was provided with 
a2 rectangular bell having carefully rounded sides. The bottom side of the 
bell was adjustable in position, and the curved surface of the test section was 
hinged to it after a 2.5-inch tangent. 

Because the velocity of the flow varied greatly in magnitude and direction, 
sven undergoing full reversal in the roller, it was considered desirable to 
utilize a single instrument that would respond in the same manner to flow in 
sither direction. This precluded the ordinary type of Pitot tube, but either 
a Pitot cylinder with opposing holes or a simple hot-wire anemometer would 
lave sufficed. The instrument actually used was a double Pitot, symmetrical 
ipstream and down, formed of two L’s of 0.04-inch tubing placed back to back 
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with stagnation openings 0.6 inch apart. The tube was calibrated in an air 
stream for both magnitude and inclination of the velocity vector, differential 
pressures being indicated on a precision manometer reading to 0.001 inch 

of alcohol. In operation, the tube was inserted through each of a series of 
centerline openings at 6-inch intervals in the plane top of the duct and a 
velocity traverse made with the tube in its normal position. On the assump- 
tion that the velocity indication corresponded to a zero inclination, the loca- 
tion of the stream lines for equal increments of the stream function W was 
evaluated from plots of the velocity distribution (see Fig. 3) according to the 
relationship 


p = [tay (22) 


0 


With the approximate angularity thus indicated by the stream-line pattern, 
the velocity readings were corrected by means of the angularity calibration, 
and the determination of stream-line location was repeated. A single cor- 
rection of this nature usually proved sufficient. 

Supplementary measurements for the determination of the local intensity 
of boundary shear were made with a 0.04-inch stagnation tube in contact with 
the plane boundary. The tube had been calibrated in connection with other 
boundary-layer measurements, and Preston’s method of shear calculation(5) 
was followed in both instances. The various characteristics of turbulence 


y u'2, V vi2y wie u'v', and (Ou'/@t) ) were measured with a constant- 


temperature hot-wire anemometer,(6) the operation of which is also too 
specialized for further comment here. It should be remarked, however, that - 
though designed for a precision of +2% under ideal conditions - the hot wire 
‘like any means of turbulence measurement) must be expected to have a 
orobable error of +5-10% in practice and to range as high as +20% in zones 

of abnormally great intensity of fluctuation. ce 

Since the efficacy of the air-model method depends entirely upon the 
Jegree of similarity between the model and its water prototype, rather ex- 
ensive preliminary tests were conducted to determine both the mean-flow 
characteristics of the jump itself and the accuracy with which they could be 
simulated in the air duct. The jump was produced at various Froude num- 
sers in the 1-foot glass-walled flume of the Institute at about the same scale 
as the air duct (i.e., dg = 1 foot). The form of the free surface was deter- 
mined by point gage, the pressure distribution on the bottom by water mano- 
meter, and the end of the roller by surface float. The several length ratios 
hus obtained are plotted in Fig. 4. For the Froude number F =4 the 
velocity distribution was also measured by means of the double Pitot tube 
‘using the technique of flushing before each reading to eliminate errors due 
© entrapped air). The stream lines are plotted in Fig. 5, at distorted scale 
© exaggerate the vertical displacement. Superposed on the latter diagram 
ure the corresponding stream lines for flow in the air duct at a comparable 
jepth ratio, surface profile, and Reynolds number, together with the pres- 
sure distribution on plane and curved surfaces. The agreement between the 
model and prototype results is seen to be quite satisfactory. 

Inspection of Fig. 4 will show the usual slight deviation of the measured 
lepth ratio from that computed by means of the simplified momentum rela- 
ionship of Eq. (1). This is attributable primarily to the neglect of boundary 
shear in the course of the derivation. In fact, it was the measured depth 
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Fig. 5. Comparison of mean patterns of flow in 
model and prototype for F = 4. 
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ratio rather than the computed that had to be used in setting the duct profile. 
Moreover, in order to reduce wall shear - and hence its influence upon the 
flow pattern - to a minimum, after introductory measurements the original 
duct was rebuilt to a width of 2.5 feet. All further data presented herein were 
obtained with the wider flow section. Even then, however, wall effects upon 
the centerline flow were still noticeable, for the unit discharge q = fa u dy 
gradually increased with distance downstream. As a first-order correction 
for this variation, the local rather than the initial unit discharge was used in 
computing the reference velocity; in no case did the change exceed 6%. Un- 
fortunately, correction for the volume flux is not the same as that for either 
the momentum or energy flux, the only means of handling all together being 
to flare the duct laterally in proportion to the boundary-layer growth. 
Following the preliminary tests, systematic measurements were made of 
the primary and secondary flow patterns for Froude numbers of 2, 4, 6, and 
8 - though without attaining a satisfactory degree of approximation at the 
highest of these. In each instance the downstream dimension was held con- 
stant at 1 foot, and the upstream dimension and the surface profile were 
varied in accordance with the prototype requirements. Reproduction of the 
necessary conditions was found to become increasingly difficult as the Froude 
number increased. In no event, of course, was it possible to obtain exact 
agreement of surface profile, pressure distribution, and flow pattern (as 
indicated by the length of roller) simultaneously, and hence each case in- 
volved slight inaccuracies in all three factors to the end of obtaining the best 
overall duplication. Barely perceptible inclination of the short tangent sec- 
tion of the bell was of considerable use in this regard. However, no such 
artifice sufficed to produce a state even roughly approximating similarity at 
a Froude number of 8, for despite the use of intermediate splitter vanes the 
rapidly divergent flow could not be stabilized either laterally or longitudinally. 
The basic results for the other three Froude numbers are plotted in Figs. 
6 and 7, the former diagram embodying the mean-flow data and the latter the 
turbulence data. (The w' component of turbulence is not shown, for reasons 
of simplicity; it did not differ greatly from v', and the measurements, more- 


—— i ae ' a 
over, indicated that a ul? ay = (v!2 4 w'2) - i.e., e v'® ay = 
: A , 


dd 6) 
2 4 ul® ay Values of the time derivative of u' are introduced later.) 
0 


Wheras the initial depth is commonly used as the length parameter in deal- 
ing with the jump as such, in comparing the results for various Froude num- 
bers it proved to be more convenient to use the final depth, as this maintains 
essentially the same longitudinal scale. Purely geometric ratios are hence 
given in terms of dg. The Froude number necessarily involves dy, of 


course, and all other parameters are also expressed in terms of the initial 
depth and velocity. 


Interpretation of Results 


Figures 6 and 7, even without further manipulation of the data, display 
trends of great significance in both the primary and the secondary aspects 
of the flow. The roller, first of all, is seen to become progressively longer 
as the Froude number increases. From the curves of velocity distribution 
three pertinent loci can readily be traced: that of zero velocity, passing 
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through the middle of the roller; that of maximum velocity gradient near the 
roller edge; and that of maximum velocity, essentially midway between the 
roller and the bed. Examination of the curves of turbulence intensity will 
show that this is a minimum where the velocity gradient is a minimum, and 
vice versa. The intensity of turbulent shear is likewise highest in zones of 
maximum velocity gradient, and the distribution curves of -u'v' are seen to 
undergo a change in sign with the velocity gradient much in accordance with 
the Boussinesq-Prandtl approximation,(7) 


OU | 
oy 


The data presented in these figures may be utilized directly to evaluate the » 
several terms in the momentum nelaHouse? of Eq. (13), reduced to non- 
dimensional form through division by oU,% dy =U 4q. The relative magni- 
tudes of the terms can be shown most conveniently by calculating and plo 
the sum 


OU ~ 1520U 


wy a oy (23) 


as a function of longitudinal distance, as shown in Fig. 8. Since the sum 
should be a constant at all successive sections, departure from its original 
value is a measure of the experimental (and perhaps computational) error. 
Of principal significance in connection with these results is the fact that the 
contribution of the turbulence to the momentum flux, although of considerable 
importance at intermediate sections, is small at each final section. The con-- 
tribution of the bed shear, though likewise small, necessarily increases with 
both distance and Froude number (and, it should be remarked, with the vis- 
cous and roughness effects as well). 

A comparable analysis of the energy relationship for the mean flow - 
Eq. (20) - proceeds in much the same fashion. On the assumption (verified 
experimentally) that the terms involving the mean veo stresses are of 
much smaller order than the others, and with dU 43/2g = qU12/2g as com-: 
mon denominator, it is now the non-dimensional sum 


that should be constant from section to section. Calculation of the various 
quantities from the data presented in Figs. 6 and 7 yields the curves shown 
in Fig. 9. The gradual change that characterized the corresponding sum in 
the momentum calculations is now even more pronounced - particularly for 
F =4. No explanation is evident beyond experimental error, including the ef- 
fect of wall retardation upon the central flow. 
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Fig. 8. Dimensionless plots of momentum balance. 
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Since it reappears in the energy equation of the turbulence, it is the 
turbulence -production term - the last member of the foregoing sum - that 
remains of primary interest. In fact, its variation in comparison with that of 
the turbulence -dissipation term is particularly significant - not only from 
section to section but over each section as well. For this reason the non- 
dimensional form of the local rate of production 


(the double integral of which was used in preparing Fig. 9) is shown in detail. 
in Fig. 10. Superposed on the same figure is the corresponding distribution 
of the local rate of dissipation 


the measured mean-square temporal gradient (@u'/Ot)2 having been used 
> Sh areeamagey the required mean-square spatial gradient through division 

. The integrals of these distribution curves over the vertical sections 
are plotted against longitudinal distance in Fig. 11, together with the corres- 
ponding integral for the rate at which the kinetic energy. of the turbulence is 
convected past each section: 

\ a 
Ud, f V“ady 


Use of the isotropic value K=15 would have yielded dissipation curves 
differing only slightly from those plotted. However, estimation of the re- 
maining terms of Eq. (21) for the final section of the jump showed them to be 
so small that the cumulative production and dissipation of turbulence energy 
through the jump (the areas under the curves of Fig. 11) could be considered 
essentially equal. This permitted bulk values of K to be calculated for each 
of the Froude numbers, the results being 12.2, 13.5, and 15.0 for 2, 4, and 6, 
respectively. How much significance can be attached to the trend in magni- 
tude is a moot question. In any event, the curves plotted in Figs. 10 and 11 
correspond to these values. 

In studying the energy transformation of the turbulence, the non- 
dimensional sum that should be constant (i.e., essentially zero) from sec- 
tion to section is 
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The corresponding diagrams for the three Froude numbers are presented in 
Fig. 12. Not included in either the foregoing sum or the diagram itself is the 
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Fig. 10. Distribution of turbulence production and dissipation. 
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rate at which work is done externally by the viscous stresses, for calculation 
of this term (the fifth) in Eq. (21) showed it to be of much smaller order than 
the others. It was impossible, on the other hand, to evaluate either the rate 
at which work was done by the fluctuating pressures or the rate of diffusion 

of the turbulent energy, since neither the velocity-pressure correlation in- 
volved in the former nor the triple velocity correlation in the latter could be 
measured with available equipment. The small difference between the com- 
puted values of the other terms might hence be assumed to represent - except 
for experimental error - the sum of these two. 

Several related characteristics of the turbulence are to be noted in Figs. 
10, 11, and 12. In the first of these, the locus of points of maximum produc- 
tion coincides closely with the border of the roller, whereas that for maxi- 
mum dissipation lies appreciably higher, and the vertical spread of the dissi- 
sation curves is likewise greater. This emphasizes the fact that the 
turbulence is not wholly dissipated at its point of formation but is convected 
by the mean motion (including the secondary flow of the eddy) and diffused 
oy its Own mixing action. Much the same situation is seen from Fig. 11. 
Whereas the maximum rate of production per section occurs even ahead of 
the middle of each roller, the maximum rate of dissipation per section - as 
well as the maximum cross-sectional energy of the turbulence itself - occurs 
shortly before its end. As shown graphically by Fig. 12, the difference be- 
tween the cumulative production and cumulative dissipation rates up toa 
Ziven section is represented with close approximation by the rate of convec- 
tion of turbulence past that section. 

The most essential of the foregoing results are summarized in the three 
yeometric diagrams of Fig. 13. Whereas these correspond to the usual one- 
Jimensional representations of the jump, not only are they drawn to proper 
scale longitudinally as well as vertically, but they are quantitatively signifi- 
cant within the non-uniform zone of transition as well as in the uniform zones 
up- and downstream. They are, to be sure, approximations rather than exact 
renditions of the phenomenon, for the least important of the variables have 
been ignored and the discrepancies evident in Figs. 9 and 12 have been dis- 
tributed among the remaining terms. All in all, however, the inaccuracies 
are judged to be no greater than the variable effects of bed roughness and 
slope encountered in the field. 

Lowermost among the curves plotted in Fig. 13 are those showing the 
variation of depth - i.e., the profile of the free surface - for each of the 
Froude numbers. Next are those representing the addition of the velocity 
head of the mean flow to the depth, much as in the usual manner. However, 
whereas the correction factor 


d 
«=e | Viady 


0) 


sompensating for use of the mean velocity U =q/d is usually assumed to 
9e unity, its magnitude within the non-uniform zone was found to rise as high 
us 4; at the end of the transition, on the other hand, it was within 5% of unity 
m every case. The two sets of curves just discussed, it should be noted, are 
he only two included in the customary diagrams. 

Third in the sequence illustrated are the curves showing the increase in 
1ead resulting from the inclusion of the velocity head of the turbulence. For 
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he effect that the turbulence has upon the flow, the relative magnitude of this 
erm is surprisingly small, even at the sections where it is a maximum; at 
he beginning of the jump, of course, it was purposely made insignificant; at 
he end it is small compared to the maximum value that it attained, though 
still an appreciable fraction of the final mean velocity head. Next in order 
ire the curves indicating the addition of the terms for the work done by the 
teynolds stresses - quantities that are very small in zones of normal turbu- 
ence but too great to be completely ignored in zones of pronounced mixing. 
fhe only other quantity with a magnitude comparable to those already plotted 
s the head lost through viscous dissipation, and this is represented by the 
lrop in the sum of all the other terms from section to section. As has al- 
“eady been emphasized, this drop is essentially complete at the end of the 
ump. 


CONCLUSION 


Through the analysis of turbulence measurements conducted in an air- 
low model of the hydraulic jump under conditions simulating three repre- 
sentative Froude numbers, significant information has been obtained on the 
nergy transformation associated with the phenomenon. In particular, the 
art played by the roller in producing the turbulence by which the transforma- 
ion is effected has taken on new clarity. Although the measurements were 
lecessarily lacking in both precision and completeness, their analysis pro- 
‘ceeded in accordance with the basic equations of motion, and conclusions 
yased thereon can be considered to have both qualitative and quantitative 
vorth. 

As in any other type of flow past a zone of discontinuity, at the juncture 
yetween the oncoming stream and the return flow or roller of the hydraulic 
ump there is a pronounced velocity gradient, and the resulting shear gives 
‘ise to the rapid generation of turbulence. The effect of the turbulence, as 
isual, is twofold. On the one hand, the mixing process diffuses all charac- 
eristics of the flow - momentum, energy, and even the turbulence itself - 
oth laterally and longitudinally. On the other hand, the increased viscous 
shear produces a rapid conversion of mechanical energy to heat. Were it 
ot for the roller, the formation of the turbulence - and hence the dissipation 
f energy - would be minimal, as can be demonstrated in a closed duct. Ex- 
‘ept for Froude numbers well below 2, however, the hydraulic jump without 
| breaking front is physically impossible. The roller is thus an inseparable 
art of the phenomenon. 

In view of the customary practice of considering the length of the jump to 
ye the distance between the limiting sections of uniform flow, it should be 
joted that the roller extends about half this distance, and that half the energy 
f the turbulence is produced within the first half length of the roller. Be- 
ause of both the convective effect of the mean flow and the diffusive effect 
f the turbulence itself, this energy is not wholly dissipated at the point that 
t is produced, but rather some distance away both vertically and longitudinal- 
y. In other words, there is a distinct lag between the loss of energy to the 
nean flow in the form of turbulence and its ultimate dissipation in the form 
f heat. 

Of interest in this connection is the surprisingly small kinetic energy of 
he turbulence at any section, even in the zones of maximum production and 
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dissipation between the mid-section and end of the roller. At the final sec- 
tion of the jump, moreover, it is already small in comparison with the maxi- 
mum value that it attained. Hence the belief that the energy lost to the mean 
flow persists for a considerable distance in the form of turbulence is : 
definitely in error. In fact, the hydraulic jump is not only an effective means — 
of reducing mean-flow energy but a very efficient mechanism for restoring 
conditions of uniformity to the flow. 
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depth of flow 

Froude number 

acceleration of gravity 

piezometric head (sum of pressure head and elevation) 
total head (sum of velocity and piezometric heads) 


head lost in jump 
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tensor subscript, corresponding to coordinate direction of equation 


of motion 


tensor subscript, corresponding to each coordinate direction in suc- 


cession 

dissipation coefficient 

Prandtl mixing length 

length of jump 

length of roller 

pressure intensity 

rate of flow per unit width 

Reynolds number 

surface area 

time 

instantaneous longitudinal component of velocity at a point 
average of u across vertical section 

instantaneous vertical component of velocity at a point 
instantaneous magnitude of velocity vector at a point 
volume 

instantaneous lateral component of velocity at a point 
longitudinal coordinate 

longitudinal component of body force per unit mass 
vertical coordinate 

lateral coordinate 

velocity-head coefficient for energy equation 
specific weight 

Boussinesq mixing coefficient 


dynamic viscosity 


1528-30 HY 1 February, 19: 
p = mass density 

w = stream function 

A bar (~) denotes the temporal mean value of the quantity beneath it; 


a prime (') denotes an instantaneous deviation from the temporal mean. 
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ABSTRACT 


Relationships are proposed which give both the quantity and quality of the 
tal, suspended, and bed loads as functions of the stream and sediment 
aracteristics. In the process of empirically defining the relationships, an 
\couraging correlation of laboratory and field data (including some original 
periments) was obtained. 


SYNOPSIS 


Although a rigorous analysis of the general problem of sediment transpor- 
tion is not yet possible, by means of a descriptive analysis the various 
ctors involved have been isolated and the relationships among them quali- 
tively indicated. Parameters linking the hydraulic characteristics of the 
yw and the characteristics of the bed material were then formed through the 
e of appropriate approximations. The relationships between these para- 
eters, however, could only be defined empirically. 

The experimental data which were used for this purpose included original 
periments conducted at the Iowa Institute of Hydraulic Research and 
iblished data from other sources. Empirical curves were drawn for the 
tal, suspended, and bed loads, and a computation procedure devised to con- 
der the non-uniformity of the bed material. The correlation of laboratory 
ta which was finally achieved was good considering the probable errors of 
€ measurements. The proposed relationships were also used to predict the 
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sediment-transporting characteristics of three natural streams with encourag: 


ing results. ? 
‘ 


INTRODUCTION 


That the flow in a channel will exert a tangential force on the boundary of : 
the channel is well known. If the boundary is composed of discrete particles, 
it is readily apparent that those particles may thus be set in motion. Not as 
apparent is the functional relationship which must exist between the flow and 
the particle movement. The paramount example of the general problem is the 
transport of sediment in alluvial streams. Despite the fact that rivers and 
canals have been the concern of man since the dawn of history, even today 
their behavior cannot be predicted with complete satisfaction or certainty. 

Most investigations in the past have been concerned with the particle 
movement at the bed—i.e., the bed load. As a result, many formulas for the — 
rate of bed-load transportation have been proposed. None, however, has 
achieved general acceptance. More recently, the phenomenon of sediment : 
suspension has been investigated both experimentally and analytically. Aside | 
from certain secondary, but nevertheless material, considerations that are © 
as yet unresolved, the distribution of suspended sediment has thereby been — 
satisfactorily formulated. The latest studies, several contemporaneous with 
the one of which this analysis is a part, have now progressed to the most 
general phase of the sediment transport problem—i.e., the combined bed load | 
and suspended load, or total load. The experimental phase of the investigation 
of the total sediment load conducted at the Iowa Institute of Hydraulic Re- 
search is only included briefly herein. It is reported in detail, including all 
pertinent experimental data, in a final report to the Office of Naval Research, 

_the sponsors of the study. 


é 


Qualitative Analysis 


Bed Load : 


The total sediment load can be divided into two parts: the bed load,in 
which the particles move essentially in contact with the fixed (or semi-fixed) | 
boundary, and the suspended load, in which the particles move entirely sur-_ 
rounded by and at essentially the velocity of the water. Although under some: 
conditions the bed load may be only a small fraction of the suspended load, — 
the bed load is fundamentally the more important, because it is necessary to 
the existence of the suspended load. 

In order to set any individual particle of bed sediment into motion, the 
flow will have to exert upon it a finite force of some certain magnitude which} 

‘depends on the shape, size, and density of the particle and its placement 
among the particles that surround it. The force exerted on this particular 
particle will depend on the average tangential force per unit area exerted by 
the flow, the position of the particle with respect to other particles in the 
neighborhood, and, if the flow is turbulent, on the particular instant of obser- 
vation. Unless, however, one subscribes to the notion that turbulent fluctu- _ 
ations can be infinitely large, there is a limiting mean-flow condition below 
which this particle will not move. More importantly, there is a limit below 
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lich none of the particles making up this semi-fixed boundary will move. 
sually this limiting flow condition is characterized by a “critical tractive 
rce,” T, = (YyQS),, (see List of Symbols), which, although real, is difficult 
define except subjectively. 

If the tractive force on the boundary is greater than the critical, some of 
e particles on the surface will move; the number, the size, and the velocity 

the particles in motion will then determine the rate of bed-load transpor- 
tion. The mode of movement will be rolling and sliding over other station- 
y particles, although, because the surface is rough, the moving particles 
uld conceivably lose contact with the boundary briefly. The force exerted 
the flow on the moving particles must be transmitted to the stationary 
ticles forming the fixed boundary. This notion does not ignore the acceler- 
ion and deceleration experienced by the particle, but recognizes that ina 
atistical steady-state condition of bed-load transport the total tractive force 
the flow must eventually be transmitted to the stationary boundary. 

Many equations have been suggested for the prediction of the rate of bed- 
ad transportation as a function of the flow conditions and the sediment 
operties. In several of the equations the total tractive force, To = YyQS, 
the only flow characteristic considered; in almost all of the others either 
e total tractive force or the slope is included, together with either the ve- 
city or the discharge. Although ae degree of similarity can be demon- 
rated among these various equations 1) (see References), the difference in 
te of transport as predicted by them is fully as great as one might expect 
om one’s first ay ten of the dissimilarity of the formulations of the 
fferent equations. 2) Unfortunately, none of the equations can be shown to 
} better than the others on grounds of either expediency or principle. All 
-some data, but none fits all data. Little or no basis in theory is claimed 
r most of the suggested equations. Rather, they are relationships which 
ve been found by cut-and-try curve fitting or more or less arbitrary 
uation forms together with experimentally determined coefficients. 

Two equations for the rate of Ga transportation, proposed compara- 
rely recently by A. A. Kalinske 3) and H. A. Einstein(4,5) have the distinction 
being considered theoretical. To the extent that both could be derived, 
gically, from explicit, but different, sets of assumptions, they can be so 
msidered. Both, however, can be criticized on a number of counts, especial- 

as to the implications of some of the necessary assumptions. For ex- 
nple, Kalinske at an intermediate step in his presentation equated the 
mber of particles in motion to the number of particles on the bed. This as- 
mption implies that the number of particles in motion is constant no matter 
at the rate of transport. The principal objections that can be raised in re- 
ect to Einstein’s development cannot be as simply stated (elsewhere they 
ve been reviewed in detail), (6) but they are equally damaging to the confi- 
nce that can be placed in his equation. Unfortunately, therefore, the “theo- 
tical” equations are no more acceptable than the “empirical” equations. In 
subsequent section still another bed-load equation will be presented which 
basically empirical, although it has some rationality—its merit being that 
is an integral part of the total load relationship. 


spended Load 


A suspension of particles heavier than the fluid (and too large for Browni- 
movement) is possible because of the mixing action of the turbulent flow. 
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In a turbulent flow there is a constant exchange of fluid masses, or volumes, © 
across planes everywhere in the field of flow. For reasons of continuity equa 
volumes of fluid must move up and down past any horizontal plane. Itis 
readily apparent that, since the particles are falling with respect to the fluid 
surrounding them whether the fluid is moving up or down, the mixing action — 
can offset the action of gravity only if there are more particles in the fluid ~ 
volume moving up than there are in the fluid volume moving down. That is to 
say, that there must be a concentration gradient such that the concentration ~ 
is greater at lower levels if the mixing action is to maintain a statistically 
steady-state condition. ; 

The equation expressing the equilibrium state of a sediment suspension is 

MS aes a (1) 

the left-hand side of the equation representing the rate at which sediment is — 
falling per unit area across a horizontal plane under the influence of gravity, 
and the right-hand side the net rate at which sediment is being lifted by the © 
mixing action of the turbulent flow at the same elevation y. The mixing coef-: 
ficient €, is obviously a measure of the rate of fluid exchange, at least as a : 
first approximation. ; 

Similarly, the equation for the apparent shear in turbulent flow can be 
written as 

dv 


ee @) 


If one disregards secondary effects such as viscous stress and virtual mass, 
it would seem that the two mixing coefficients €g and € m should be the same, 
It is possible, however, that the mixing concept as represented by these two 
equations is oversimplified and that in moving a finite vertical distance the © 
fluid does not transfer, on the average, the temporal mean conditions ofits 
point of origin. That is to say that a natural-selection process may actually 
occur in the mixing action. 

The usual expression for the distribution of sediment in the vertical is ol 
tained by assuming proportionality of the mixing coefficients (€g = = Bem; 
linear variation of the shear, and logarithmic velocity distribution.(7) The 
differential equation for the sediment suspension can then be integrated to 
give 


$ 


: : 
: if ae ) . 
Cc. Ti tig aie (3) 
where z = w/Bx Vrop . Only a relative distribution is thereby obtained, as the 
concentration at any level y is expressed in ratio to the concentration Cg at 
some level a. Obviously, if the concentration c, could be specified by some _ 


other means, the average suspended-load concentration Cg could be obtained , 
by integration of the expression 
eo 
) 


lm @ 


aes 
oO. SS = 
q 


The two factors which will dominate in determining the average concentration 
of the suspended load are cq and \r0/P'/w. The range of variation of other 
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ctors such as § and X is small in comparison to the possible variation in 
ese two. 

As can be seen from Eq. (3) the ratio of the shear velocity to the fall ve- 
city is the primary factor determining the degree of uniformity of the 
mcentration. Given a value of cg at the lower-most level of suspension, the 
ore uniform the dispersion of the sediment (i.e., the larger the value of 
To/p /w) the greater will be the value of the mean concentration Cg. Al- 
ough perhaps not observable, there should be a level, or zone, near the bed 
here the concentration becomes less dependent on the mixing action of the 
rbulent flow than on the rate at which particles are cast up, or entrained, 
om the bed. 

The assumptions made by Einstein in expanding his bed-load function to 
srmit the calculation of the total load(5) were that the reference level could 
> taken as twice the diameter of the sediment particle and that the concen- 
ation in this zone was proportional to the rate of bed-load transport divided 
; the shear velocity associated with the sediment particle and the distance 
om the reference level to the bed. Implicitly he also assumed that the 
echanism of entrainment was the mixing action of the turbulent flow which 

so effective at higher levels, and that the logarithmic velocity distribution 
nd, therefore, the concentration distribution) was valid as close as this to 
e bed. 

Earlier Lane and Kalinske(8) had suggested that, for the sediment particles 
a size class represented by a mean fall velocity w, the concentration at the 
yottom” was proportional to the percentage of the bed of that particular size 
id a function of the shear-velocity /fall-velocity ratio, Experimental investi- 
tions under their direction resulted in modifications of the originally pro- 
ysed relationship. Pien(9 simply adopted a distance from the bottom of 0.1 
e depth of flow as the reference level. Hsia, 10) pecause he used a very 
ne sediment which was almost uniformly dispersed in the flow, considered 
e bottom concentration equal to the mean concentration, but found that a 
cond parameter, dYT9/p/y, proportional to the ratio of the particle diame- 
r to the thickness of the laminar sublayer, was needed. The basic concept 
both the original and modified relationships was that the vertical turbulent 
ictuations “picked up” material from the bottom. 


terchange Between the Bed and the Flow 


A sediment particle which is suspended in the turbulent flow will follow an 
‘ratic path which will depend on the velocity of the fluid about it from instant 

instant. Each of the many particles which pass through any point in the 
eld of flow will, of course, follow a different path. Although the paths of the 
dividual particles cannot be predicted, in the aggregate a pattern of diffusion 

the particles passing through the given point will obtain. Even if the parti- 
es were of the same density as the fluid, within a nominal distance some of 
em would reach the bed of the stream because of the turbulent mixing action. 
order to maintain a steady state of suspension within the field of flow for 
ery particle that returns to the bed, a particle must, by some means, be 
‘moved from the bed and injected into the flow. 

Considering the importance of the interchange phenomenon to the under- 
anding of the phenomenon of the suspended load, there has been surprisingly 
tle speculation as to the mechanism whereby the sediment particles are re- 
oved from the bed. The two notions that have been advanced are (1) that the 
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flow around the particle results in a lift force greater than the weight of the — 
particle and that the particle consequently moves up from the bed into the — 
flow, (11,12) and (2) that the mixing action of the turbulent flow is sufficiently _ 
strong at the bed level to remove the particle from the bed, (12,13 1 

It cannot be disputed that there might be lift forces on a particle on the 
bed, because the particle is not in a uniform or even symmetrical flow field 
but in a velocity gradient, because the pressure may approach the stagnation 
magnitude under some particles, or because the velocity possesses locally a 
component in the upward direction. No matter how this possible lift force 
may come into being, however, there will also be a drag force on the particle: 
such that the particle will begin to move or, if already moving, will move ' 
faster parallel to the bed—i.e., in the direction of the drag force. Ifthe lift — 
became almost equal to the weight of the particle, any small drag force would’ 
be sufficient to roll or slide the particle. Since the forces exerted by the flow 
on the particle are the result of the relative motion of particle and fluid, they | 
could be expected to decrease once the particle is accelerated. Howa lift 
force greater than the weight of the particle could develop on a particle that | 
is free to move is, therefore, difficult to envision, except possibly through 
the action of vertical components of large scale turbulence. 

This possibility is, of course, the same as the notion that the mixing action 
of the turbulent flow extends to the level of the bed. For a smooth, solid bed 
it is readily apparent that at least the vertical turbulent fluctuations will de-— 
crease with proximity to the boundary to the limit of zero at the boundary. ~ 
Although a rough, porous boundary such as is formed by the sediment might © 
either result in, or permit, larger turbulent fluctuations than have been found 
in the vicinity of smooth boundaries, one would still expect the mixing action 
of the turbulent flow to grow weaker as the boundary is approached. Althoug 
the process cannot be categorically ruled out, it would seem that the turbulent 
exchange which can be quite effective in the interior of the flow would be : 
minimal and rather inadequate at the boundary. 

A more plausible explanation of how a sediment particle can leave the bed 
can be based on Newton’s first law of motion—that a body in motion will move 
uniformly in a straight line until acted upon by some force. Consider nowa 
particle moving up the weather slope of a dune. The forces on the particle 
are the propelling forces of pressure and shear due to the flow around the _ 
particle, the force of gravity, and the reaction force of the bed that canbe 
resolved into a normal, support force and a tangential, resistance force. Just 
before the crest of the dune all of these forces will be acting on the particle; 
just after the crest of the dune the bed forces supporting the particle and re- 
sisting the motion will be lost, but the fluid forces and the gravitational force 
will continue to act on the particle (although the fluid forces may change be-_ 
cause the flow around the particle changes). The motion of the particle will 
then depend on the velocity (both direction and magnitude) of the particle as it 
leaves the crest of the dune and the forces which thereafter act upon it. ) 

If the velocity of the particle and the propelling forces are small and the 
weight of the particle under water is large, the particle will merely roll over 
the crest and down the lee slope of the dune. This, of course, is a descriptio! 
of bed-load movement. If the gravitational force is small in comparison to 
the momentum (mass times velocity) of the particle, the particle will only 
gradually deviate from its initial direction of motion and will tend to move in 
the parabolic path of a projectile. This ideal path, of course, will be modifie 
by the effect of the fluid forces on the particle. 


ee, Sn 
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Although the prediction of the path of the particle may not be possible, it 
hould be clear that, if the velocity of the particle moving as bed load is great 
nough, the particle can leave the immediate vicinity of the bed in what might 
e described as a self-launching action. Moreover, the factors which will 
overn the motion of the particle, other than the properties of the fluid and 
he particle itself, are the form of the bed, the flow pattern (including the ve- 
ocity magnitude), the gravitational force, and the velocity of the particle at 
he crest of the dune. The amount of material involved in this action will be 
‘elated to the number of particles in motion and the size of the particles. 
ince the product of the number, size, and velocity of the particles is the 
olume rate of bed-load transportation, it is immediately apparent that the 
ed load will play a primary role in the determination of the suspended load. 

If the bed is covered with dunes, this self-launching action is easy to visu- 
Jize. Even more to the point, it can be observed. If there are no dunes— 
.€., if the bed is plane—the same action is still possible, because the bed- 
oad particles must move over roughness elements of their own size. Asa 
esult, the direction of motion of the particle can be away from the bed and, 
fits velocity is great enough, the particle can escape from the immediate 
eighborhood of the boundary. If the moving particle in the case of the duned 
ed can be likened to a ski-jumper, in the case of the plane bed it can be lik- 
med to a hot grounder in a sand-lot baseball game. 

The secondary motion, or large-scale vortices, of the turbulent flow can 
Jso result in a movement of the particles in contact with the bed whereby 
hey can be self-launched into the flow. This can be demonstrated simply in 
. glass jar filled with water with a bed of fine sand on the bottom. By moving 
_pencil in a small orbit an almost irrotational vortex can be induced. Be- 
ause of the secondary motion, the sand particles move inward in spiral paths 
nd a small, sharp-pointed, cusped dune is formed. If the velocity of flow 
ear the bed, and consequently, the velocity of the particles, is sufficiently 
reat, the same self-launching action from the peak of the dune is observable. 
The importance of macroturbulence such as this in the movement of sedi- 
nent has been noted by Matthes$14) although the mechanism whereby the 
articles are removed from the bed is not described. A similar phenomenon 
as also been noted by Lane).(15) In a stream or a laboratory flume this ac- 
ion is quite impossible of observation. It can be readily observed every 
utumn, however, with dust and fallen leaves as the sediment particles. Ona 
arge scale, of course, it is the tornado of the western plains. 


Observations 


iquipment and Technique 


Since most previous investigations of sediment transportation have been 
oncerned with bed-load movement, very few data on the rate of transport of 
\tal load, or even suspended load, were available at the time this study was 
nitiated. An experimental phase of the investigation, therefore, was essential 
) the determination of quantitative information relative to this most general 
ase of sediment movement. Because the data which was gathered were so 
mportant in the development of the relationships for the sediment load, a 
rief description of the experimental phase is included herein. hen 

The principal item of experimental equipment was a recirculating, tilting 
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flume with an overall length of 105 feet, shown in Fig. 1. The test section 
was 90 feet long and had a cross section 36 inches wide and 18 inches deep. 
In order to permit visual observation of the flow, the sides of the test sectio 
were made of 1/4-inch plate glass. 

Because the resistance to flow will depend in large measure on the rough 
ness of the self-formed alluvial bed, the discharge and the depth of flow (and 
hence the velocity) are usually fixed in experiments such as these, and the 
slope is allowed to develop as required. In order to speed up the deposition- 
erosion process whereby the slope is established, it is expedient to be able — 
to tilt the entire flume during operation. The mechanism which was used to : 

} 


: 
: 


accomplish this was a system of interconnected cams. 

The rate of flow was controlled by the speed of the pumps and measured ~ 
by calibrated, segmental orifices in each line in conjunction with an air-water 
manometer. The depth of flow was determined by the volume of water in the” 
system and the water-surface elevation was measured at 10-foot intervals ] 
along the test section by piezometers connected to an open manometer board ; 
The water-surface and sand-bed elevations relative to the flume could be de= 
termined by means of point gages mounted on carriages which in turn slid on 
stainless-steel rails which were parallel to the flume bottom. The slope of © 
the flume itself was obtained by measuring the travel of the lower end with a 
vernier and scale made from a standard point gage. : 

Values of velocity and concentration at a point were obtained by means of © 
the combination instrument shown in Fig. 2. Two Pitot tubes, modified to 
produce a larger differential reading, were mounted on each side of a rec- _ 
tangular sampler nozzle. The differential-head indication was obtained with © 
an air-water manometer. A small variable-speed pump was used to establish 
any desired rate of flow into the sampler-intake nozzle. The rate of flow was’ 
measured by a calibrated Venturi meter and air-water manometer. The 
pumped sample was collected at first in quart milk bottles; later large glass” 
tubes were obtained to simplify the handling of the sample. The amount of 
sediment in the samples was determined by weighing in calibrated pycno- 
meters. For small concentrations a standard 10-ml pycnometer with a fitted 
stopper having a capillary hole was used and for large concentrations a 500-— 
ml flask with a narrow, graduated neck. ; 

In order to measure the total load, a low submerged weir was installed at 
the downstream end of the test section. The weir was formed of aquarter ~ 
cylinder cut from a nominal 4-inch brass pipe, and placed so that the convex 
side was upstream and a tangent to the top edge, or crest, of the weir was . 
parallel to the flume bottom. The crest of the weir was approximately an 
inch above the mean elevation of the sand bed and slightly above the crest of 
the dunes. 

In a study such as this the operation and the measuring techniques must 
be designed first to establish normal equilibrium transport conditions and 
then to obtain a reliable and complete description of those conditions. Al- _ 
though equilibrium here refers specifically to a state of flow in which there — 
is neither aggradation nor degradation of the bed, it also implies that there is 
no change in state with time. Because the rate of change from a non-equi- i 
librium to an equilibrium state will depend primarily upon the degree of non= 
equilibrium, the desired steady state will tend to be approached asymptotic- 
ally. D 

Arbitrary, but judicious, operation of a sediment-transport flume can ma- 
terially reduce the time required to attain a normal equilibrium condition. 
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9 at 5'=45' 
Fig. 1. Flume. 
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To confirm the establishment of equilibrium, however, repeated measure- 
ments over considerable periods of time are necessary. To verify the uni- 
formity of a reach, measurements at successive sections are required. 

Given a width of section (which is fixed by the construction of the test 
section) and a bed material (which can be chosen arbitrarily) only two of the 
three primary flow characteristics (discharge, depth, and slope) can be con- 
sidered to be independent variables. In a natural stream the slope and the 
discharge are the imposed conditions and the depth is conditional upon these 
and upon the configuration of the stream. Operation of a recirculating, tilting 
flume is simpler if the discharge and depth are chosen as the imposed con- 
ditions and the other characteristics of the flow are allowed, or helped, to at- 
tain their natural value or form. 

In the first runs measurements were taken of all pertinent characteristics 
of the flow as the condition of equilibrium was approached. The slope proved 
to be the most sensitive index, as it was affected by both the non-uniformity 
of the flow and the roughness of the bed. As deposition and scour occurred, 
remolding the overall bed so that the rate of sediment transport at every 
section was the same, there was a concomitant change in the water-surface 
slope. As the natural ripple or dune pattern developed, the roughness of the 
ved changed. Therefore, the resistance to flow changed, and, consequently, 
the slope of the water surface, The velocity distribution and concentration 
listribution were affected by the position of the dunes relative to the measur- 
ing section—especially, of course, near the bed. As a result, the spatial vari- 
ation of these measurements was of the same order as the temporal variation, 
inless the flow was adjusting itself quite rapidly toward the normal equilibri- 
im condition. 

For the coarser of the two sands used, the slope proved to be a sufficient 
ndicator of the attainment of the normal equilibrium condition. When the bed 
1ad been remolded and its configuration established, the slope of the water 
surface no longer varied with time. (Tilting the flume, in effect, simply re- 
juced the time necessary for remolding the bed.) A reach of approximately 
30 feet, from about 35 feet to about 85 feet from the entrance, would then ex- 
1ibit uniformity—as defined by a constant slope and constant depth of flow. 

For the finer sand the water-surface slope alone was not found to be an 
wiequate measure of the state of flow. This was because the length of the 
reach displaying uniformity was only about 30 feet, and, therefore, the 
neasurement of the slope could not be as precise. Measurements of the 
s3and-bed slope and elevation were found to be necessary also. In order to 
neasure the sand bed the flow had to be stopped and the run interrupted. 
Therefore, the sand-bed measurements were only used as a final check after 
he water-surface slope indicated stability. 

Fortunately, because the total time required for a single complete run 
nade continuous operation impractical, interruption of a run did not affect 
naterially either the attainment or the state of equilibrium. If care was taken 
n stopping and then resuming a run, the state of flow at the time of inter- 
uption was reestablished within an hour at the most. As has been mentioned 
efore, the criterion for normal equilibrium flow is that the characteristics 
f flow do not change with time. Since this state is approached asymptotical- 
y, it is obvious that the final determination of whether the desired state has 
een reached must rest on the judgment of the operator. When normal equili- 
rium had been attained—in the operator’s judgment—measurements were 
nade of water-surface elevations, velocity distribution, concentration distri- 
ution, and sand-bed elevations. 
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For the coarser sand, only centerline distributions of velocity and concen 
tration were determined. For the finer sand, complete traverses during | 
several runs were made across the flow section in the uniform zone and over 
the low quarter-circle weir. Measurements of the dune size could be obtained | 
in the case of the coarser sand. However, for the very fine sand the material 
settling out of suspension almost obliterated the dune forms. ‘ 

Water-surface and sand-bed elevations could be read to 0.001 foot. Al- — 
though errors in any single determination could be greater, the error which ~ 
might be expected in the average of repeated measurements should not exceed 
this value. The depth of flow ranged from about 0.25 to about 1.00 foot. At 
the smallest depths the percentage error, therefore, could be of the order of | 
1%. The slope ranged from about 0.0004 to about 0.0018 over the 50-foot 
reach of the coarser sand and from about 0.0008 to about 0.0012 over the 30- ; 
foot reach of the finer sand. Errors as large as 10%, therefore, are possible, 

Few of the individual velocity determinations in the interior of the flow 
should be in error by more than 2 or 3%. Near the sand bed, however, the 
possible error in the measurements became larger as the boundary was ap- 
proached, because the velocity became smaller. Moreover, the position of 
the dune relative to the measuring section was unknown, so that the meaning © 
of the measurements near the sand bed was somewhat indefinite. : 

The samples which were taken for the determinations of point concen- 
trations were large enough to give a good time average. Fig. 3 shows the : 
long-period fluctuations of the concentration at a point as a function of time. ~ 
The short-period fluctuations have been averaged out because each successive 
determination of concentration was made from a 1-quart sample, each quart 
being equivalent to a stream tube almost the length of the test section. In the 
first runs 5-quart samples were taken to reduce this source of error inthe — 
concentration measurements. The slightly larger 5000-cc tubes, later substi- 
tuted for the milk bottles, were equivalent in volume, and the procedure of 
transferring the sediment sample to the pycnometer was simpler. 
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Experimental Results 


The size-frequency distribution of the two sands which were used is shown 
in Fig. 4. Despite the fact that the finer sand is within what is generally con- 
sidered the silt range, it is a true sand. Both sands were fractions of the 
Ottawa deposit (the coarsest fraction of which is the standard concrete test- _ 
ing sand) and are almost pure silica. | 

The experimental results for the coarser sand are summarized in Table I. 
A more detailed summary is included in the report to the sponsor, (16) In re- 
gard to the basic problem being investigated—the correlation between the 
sediment load and the flow and sediment characteristics—little can be learned 
from comparing the pertinent values for different runs, even with the aid of 
simple plots. An empirical relationship describing this correlation was 
found only after a qualitative analysis of the problem had indicated parameter 
which could be taken as expressing the factors governing the phenomena in- 
volved. This qualitative analysis and the relationships finally obtained are | 
the subject of the following section. Several general characteristics of inter- 
est, however, are immediately apparent from the tabulated summary. | 

The recorded dune heights h and 1 show that within the range of conditions 
represented here, the size of the dunes did not vary greatly—except that in | 
one run a plane bed was formed. (A number of other attempts to obtain a . 
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Fig. 4. Size-frequency distribution of sands. 
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plane bed were unsuccessful. The combinations of velocity and depth giving 
Froude numbers of about 0.8 resulted in very unstable flow conditions in the 
flume.) As can be seen in the respective tabulations, the relative roughness 
2yo/h of the duned bed exhibited a greater variation than did the absolute 
roughness h, : 
In Fig. 5 the Weisbach resistance coefficient is plotted as a function of the _ 
relative roughness together with the relationship for fully rough-pipe flow 
based on the Nikuradse experiments with sand-roughened pipes. Despite the © 
scatter of the experimental points a very good agreement is evident between H 
the two kinds of roughness in the two flow systems. In fact, one should proba= 
bly interpret the superposition of the points on the curve as due to a fortui- : 
tous choice of a typical roughness length until further evidence is available 
as to the effect of roughness shape and spacing on the resistance coefficient. — 
One can conclude, however, that the bed configuration is of major importance” 
in the resistance to flow of alluvial streams. ‘ 
Similarly the roughness has an effect on the velocity distribution, which— ~ 
except for the plane bed—is in accord with the Nikuradse experiments. Fig. 7 
6 shows the power-law exponent as a function of the relative roughness to- } 
gether with the Nikuradse relationship. The values of the exponent m were 
obtained from log-log plots of the velocity distribution. Semi-logarithmic 
plots of the velocity distribution were also made and the values of K obtained © 
therefrom are listed in the table. No systematic correlation for the variation 
of K was found with either the concentration or the roughness. : 
Similarly, no pattern of variation was found for the values of 8 listedin ~ 
the table. The value of 8 was taken as the ratio of the theoretical exponent of 
the concentration distribution z; to the measured exponent z as obtained from 
log-log plots of the concentration distribution. In the theoretical evaluation 
of z, the fall velocity obtained in a bottom-withdrawal-tube determination of © 
the size-frequency distribution and the nominal value of 0.4 for K were used. 
The effect of the concentration would be to decrease the fall velocity and, 
therefore, z; by less than 10%. Use of the experimental values of x would 
increase z, and, therefore, 8. For most of the runs this increase would be 
of the order of 30%. 1 
Table II summarizes the results for the finer sand. From the tabulated f 
values in this table it can be seen that the resistance to the flow for these 
runs with the finer sand is about the same as for those with the coarser sand. 
Although meaningful measurements of the roughness could not be obtained be= 
cause of the deposition of the suspended load, it could be observed that the 
size of the dunes was approximately the same as for the previous series of 
runs. The velocity distribution characterized by the exponent m as obtained 
from log-log plots is also indicative of a rough boundary similar to the runs © 
with the coarser sand. The values of kK from semi-logarithmic plots listed in 
Table I are slightly smaller for this series. A comparison of theoretical i 
and measured concentration distribution by taking a ratio of the appropriate 
z values results again in B values greater than unity, and about the same as 
those found for the coarser sand. Correction of the fall velocity for the 
concentration effect would decrease B by as much as 40%, but use of the 
experimentally determined values of k would almost double the value of B. : 
Of perhaps greater significance is the fact that the total load measured 
over the weir was almost the same as the suspended load measured in the ; 
normal sections. In fact, the errors in measurement of the suspended load 


were larger than the bed load. It is estimated that the bed load for these runs 
was probably not more than 1% of the total load. | 
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Development of the Relationships 

7 


Given a flow condition (mean velocity, depth, and slope), and a sediment ' 
(mean size, frequency distribution, specific gravity, and shape), there will be 
a corresponding rate of sediment transport, partly as bed load, partly as sus-_ 
pended load. The qualitative discussion of the mechanism of sediment ? 
movement in the preceding sections was for the purpose of selecting the f 
factors which determine the rate of transport. If the hydraulic and sediment } 
characteristics can be linked in parameters descriptive of these factors, then 
at least empirical relationships for the sediment load may be obtained. (At — 
the present time any attempt at a rigorous analysis is defeated at the outset ~ 
because the fluid forces on the individual particles cannot be specified in suf- 
ficient detail.) One of the desired parameters is readily apparent—the shear= 
velocity /fall-velocity ratio, YT 9/p/w. This ratio expresses in large 
measure the effectiveness of the mixing action of the turbulence and, of a 
certainty, should be included in the relationship for the suspended-sediment 
load. 

The other desired parameters cannot be as confidently stated. Especially 
needed is a clear understanding of the forces acting on the particles in motion 
as bed load. For the present it is practically necessary to relate the forces ' 
causing motion to the tractive force, and the forces resisting motion tothe — 
critical tractive force for the beginning of movement. Only in the case of the 
plane bed, however, can the total tractive force be considered significant. If 
there are dunes on the bed (the more common condition) the part of the total © 
tractive force which is resisted by the tangential component of the integrated 
pressure distribution along the duned boundary cannot be considered effective 
in causing the sediment particles to move. An approximate measure of the 
remainder of the tractive force can be obtained by the use of the Manning ; 
formula and Strickler’s expression for n as a function of the sediment diame- 


ter :(17 
ve q/s | 
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Although this expression may be only a crude approximation for the shear in 
the rather complex flow along the weather slope of the dune, it has the ad- 
vantage of being computationally simple. 

An expression for the critical tractive force can be written as 


{ 
é 
' 
; 
: 
i 
tT, =HCd oF 
where C is a coefficient dependent on the sediment characteristics andthe 
flow conditions near the boundary. Both T9 and T, are representative of the : 
_ average shear per unit area of the bed, rather than the force on the particle 
which is moving or about to move. Although neither is a true measure, even | 
in approximate form, of the desired forces, one might conjecture that each is. 
related to these forces in much the same manner. Thus, the ratio T9/Te can 
be assumed to be a significant parameter descriptive of bed-load transpor- _ 
tation. Moreover, if the rate of bed-load transport is a primary factor in the 


interchange phenomenon, as seems likely, this ratio would also be important 
to the suspended-load movement. 
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A third parameter, which is also suggested by the discussion of the pre- 
ceding sections, is the ratio of the velocity of the sediment particles moving 
as bed load to their fall velocity. This ratio would be a measure of the 
tendency for self-launching and would be important for suspended-load rather 
than bed-load movement. Similarly, the form of the bed might be descriptive 
of the opportunity for self-launching. In the final analysis of the data, the 
influence of such a self-launching parameter could not be assessed—partly 
because of the scatter due to experimental errors, and partly because it would 
necessarily resemble the other parameters. 

To test the significance of the ratio 23 /Te the data obtained in this study 
were plotted as ¢ = 265 q./q against Rava and an almost linear variation 
was found for each sand. For reasons more intuitive than rational, the factor 


tip 4 ( a 


Vig *\% 


Yo (7) 


was included in the relationship and a plot made of 


c ‘Je 
ree sy A 
bah ant, Ply (8) 
Jo Se 
using the data of this study and other published data. The resulting relation- 
ship is shown in Fig. 7. The solid curve drawn through the points represents 
the total load and the dashed curve the bed load. Table III summarizes the 
source and the pertinent information for the various data plotted in Fig. 7. 
Only those runs for which the particle shear approached or was less than the 
critical shear were rejected. A simple ratio like 1 o/ Tes of course, cannot be 
expected to describe conditions at this limit. 

For the data of Toch and Hsia, C values for the critical tractive force of 8 
and 16, respectively, were used; otherwise a value of 4 was employed. These 
values are plotted together with Shields’ curve in Fig. 8. The values used are 
approximately in agreement with the recommendations of White, 18) and for 
values of d/6' above 0.2 with Shields (19) as well. For smaller values, where 
the laminar sub-layer completely envelops the particles, the variation of the 
critical tractive force with the diameter of the particle is not clear. It should 
be noted that Shields’ data did not include values of d/6' much less than 0.2. 

It was noted, especially on the larger work sheet from which Fig. 7 was 
prepared, that there was a systematic scatter to the points for any series as 
identified by the different symbols. It was surmised that this tendency for 
the sets of points to cross the mean curve might be due to the non-uniformity 
of the sediments, and a computational procedure was devised to consider the 
composition of the bed material. This computational procedure is represent- 


ed by the equation 


sea) (2) (47) ow 


Yo 
wherein the contributions of each fraction p of size d are summed to obtain 
he mean contribution. Values of T, and w were determined for the mean di- 
ameter d of the fraction, but the same values of Tg and ¥79/p'were used for 
all fractions, TO being determined from the mean diameter of the total sedi- 
ment sample. 
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Fig. 9 shows a comparison between the measured mean concentration and 
the concentration as computed in this manner. A tendency for the plane beds 
to have a measured concentration less than computed could be noted, as well 
as a hint that there might be a systematic scatter with the ratio To/ To. Itis 
readily apparent that this ratio could be taken as a measure of the roughness 
of the bed; together with the other parameters which have been introduced it 
might also be taken as descriptive of the ratio of the velocity of the bed-load 
particles to their fall velocity. The unsystematic scatter due probably to 
errors in measurement, was too great to assess the possible influence of 
these two factors. Actually a rather good correlation is indicated over the 
ten-thousand-fold range in concentration represented in this plot. 

Equation (9) also allows the prediction of the composition of the sediment 
load (bed load, suspended load, or total load depending on the function of the 
shear-velocity/fall-velocity ratio employed). Fig. 10 shows a comparison 
with measurement of such a prediction for several runs. 

Although the data used in Figs. 7 and 9 represent a large variation in sedi- 
ment size and concentration, the measurements were all made in laboratory 
flumes under controlled conditions. Therefore, an attempt was made to pre- 
lict the sediment load under field conditions using the data published by H. A. 
Einstein for Mountain Creek in South Carolina and West Goose Creek in 
Mississippi, (20) and by the USGS for the Niobrara River near Cody, Nebras- 
ca.\21) The results are shown in Figs. 11, 12, and 13, respectively. 

In the upper part of the figures the velocity and slope data are plotted to- 
zether with the mean curves which were used in conjunction with Eq. (9) to 
obtain a predicted sediment load. On the left a curve, or curves, for the pre- 
licted load and the data of the measured load are plotted as a function of the 
lischarge. On the right are plotted the predicted composition of the sediment 
load and the measured composition of the bed material. 

In the case of Mountain Creek the predicted load was in general less than 
measured. It may be noted, however, that artificial flash floods were used 
© obtain much of these data (Fig. 11), with a difference in conditions between 
rising and falling stage as indicated in the figure. The measured points of 
he largest flood are connected in the figure by light lines indicating the 
emporal continuity. The measured difference between the rising and falling 
stages is roughly comparable to the predicted difference. Of perhaps greater 
significance is the discrepancy between the predicted and measured compo- 
sition of the bed load. The predicted mean size of the bed load was 0.6 mm 
iS compared to a measured value the same as the bed material, or 0.9 mm. 

f a smaller value for T, had been used, the predicted load would, of course, 
lave been greater and the composition of the bed load would have more close- 
y approached that of the bed material. 

The predicted bed load for West Goose Creek (Fig. 12) agrees very well 
vith that measured for lower discharges, but it is considerably too large for 
ligher discharges. It is noteworthy that the discrepancy is pronounced only 
it velocities of flow above 2.5 fps. Based on the experience of the writer, it 
seems doubtful that a slot 2 feet in width such as was used would capture - 
nuch more than half the load at this high a velocity with so fine a sand (0.3 
nm). The predicted size of the bed load for West Goose Creek is only slight- 
y smaller than the size of the bed material, and is well within the scatter of 
he measured composition of the bed load. 

The excellent prediction of the sediment load for the gaging-station section 
f the Niobrara River (Fig. 13) does not include some fifty-nine measurements 
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the range of 200 to 2000 tons per day which would fall on the computed 
Ive or on points next to the curve. In order to compare the predicted and 
easured composition of the suspended load, seventeen measurements which 
cluded a large amount of fine material (the result of surface erosion) were 
jected. It should be noted, however, that the measured values of concen- 
ation have not been corrected for the inherent sampling error due to the 

ct that the sampler does not traverse the entire depth. 


plicability of the Proposed Relationships 


The proposed rate-of-transport relationships, Eq. (9) and Fig. 14, give 
th the quantity and the quality of the total, suspended, and bed loads as 
actions of the basic hydraulic characteristics of the stream and the charac- 
ristics of the bed sediment. Thus, they can be applied directly and simply 
the prediction of the sediment-transporting characteristics of a stream 
om other measurable or computable characteristics. However, it should be 
pt in mind that the relationships are basically empirical and hence can be 
ed with confidence only within the range of conditions for which they have 
en tested against actual measurements. 

In order to increase the reliance that can be placed on their applicability, 
© proposed relationships must be tested against reliable measurements 
er an ever wider range of conditions. Special field studies would be par- 
ularly valuable—not only to extend the scale to include large rivers, but 
so to assess the effect of such stream characteristics as channel shape and 
ignment. Such continued testing could result in the addition of secondary 
rameters as well as some modification of the present relationships between 
e basic parameters. 

Because approximations had to be made in the formulation of those para- 
eters, continued study of the factors entering into the sediment-transport 
lationships is also indicated. For example, from Fig. 8 it is quite evident 
at the critical tractive force relationship for fine materials is not adequate- 
understood. An investigation of the magnitude and distribution of the forces 
a boundary of composite roughness such as a duned bed is also needed. It 
readily apparent that studies such as these would permit better formu- 
‘ions of the parameters governing the sediment load. Other investigations 
at would be desirable concern the mutual interference of particles of vari- 
S sizes and the formation of dunes and ripples. 

Unfortunately, the need for a means of predicting the sediment-transport- 
s characteristics of a stream cannot wait on the completion of all desirable 
search. Since considerable rationality can be ascribed to the parameters 
the proposed relationships, and since their ability to correlate measured 
oratory and field data has been demonstrated, it would seem reasonable 
d proper to attempt to apply the relationships to field problems—so long as 
sre is full realization that they are new and comparatively untried. 

For any stream which has been gaged for a period of years the hydraulic 
stors needed for the prediction of the sediment load should be available. 
ormation as to the characteristics of the normal bed material should then 
rmit the construction of load-discharge relations such as those in Figs. 11, 
,and 13, Sheet erosion during heavy rains will contribute a temporary fine 
iction to the normal bed material. The “wash” load which results there- 
ym cannot be predicted from the proposed relationships unless the tempo- 
ry fine fraction can be estimated. For streams of which the gaging 
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Fig. 14. Relationships for sediment load. 


SCE SEDIMENT LOAD 1530-33 


‘ogram has included suspended-sediment measurements it may be possible 
ith the aid of the proposed relationships to correlate the wash load with the 
atershed and rainfall characteristics such as soil type, cultivation, basin 
opes, season, intensity of rain, antecedent rain, etc. The utility of sucha 
rrelation of unmeasured streams is obvious. 

In the case of degradation studies, the influence of the wash load should be 
nall; since the bed material can be considered known, the rate and location 
the degradation should be predictable. The case of aggradation studies is 
ore difficult, because the source of the excess sediment load will determine 
S composition. In this respect the aggradation problem is similar to the 
coblem of wash load on poised streams. 


CONCLUSIONS 


Through the use of (1) a qualitative analysis, (2) original experiments of a 
ecialized nature, and (3) supplementary data from other sources, empirical 
slationships for the primary aspects of the sediment load have been obtained. 
hese relationships permit the composition and the rate of transport of the 
ital load, the suspended load, and the bed load to be evaluated from the hy- 
saulic characteristics of the stream (mean velocity and depth of flow and 
1ergy gradient) and the characteristics of the bed sediment (frequency distri- 
ition of size and fall velocity). 

In the process of defining the relationships between the parameters which 
ere found to govern the sediment loads, a correlation of laboratory data 
presenting a ten-thousand-fold range in rate of transport was obtained with 
scatter probably not much greater than the error in the experimental obser- 
itions. The extent to which the proposed relationships could predict field 
mnditions was demonstrated for three natural streams. The degree of ap- 
-oximation was especially encouraging, since the ultimate goal of the search 
r a general sediment-transport function is its application to practical engi- 
sering problems. 
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List of Symbols 


a elevation of reference concentration Cg, ft 
c sediment concentration at a point, percent by weight 
é mean concentration, 265 q,/q (subscripts s, b, t refer to suspend- 


ed, bed, and total load, respectively), percent by weight 
Ca concentration at reference level a, percent by weight 


coefficient relating critical tractive force to sediment size 


d diameter of sediment particle (mean diameter of fraction p of bed 
material), ft 
qa mean diameter of bed material, ft 


f Weisbach resistance coefficient 
FE Froude number, VNeyo” 
h dune height, ft 

1 dune length, ft 

m exponent in power-law velocity distribution 
n Manning roughness coefficient 

p fraction of bed material of diameter d 

q rate of flow per unit width, Vyg, cfs/ft 


ds volume rate of sediment transport per unit width (subscripts s, b, 
t refer to suspended, bed, and total load, respectively), cfs/ft 


hydraulic radius, ft 


R 

R Reynolds number, 2Vyo/ Y 
Ss energy gradient, ft/ft 

Vv 


velocity of flow at a point, fps 


Vg surface velocity obtained by extrapolating power-law velocity 
distribution, fps 
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Vv mean velocity, a/Yo, fps 
Ww fall velocity of sediment particle, fps 


y elevation above stream bed, ft 
Yo depth of flow, ft 
Z measured exponent for concentration distribution 


Z 1 theoretical exponent for concentration distribution 


B coefficient of proportionality between €, and € ,, 
uf specific weight of water, lb /tt? 
6' thickness of laminar sub-layer, ft 
Sam mixing coefficient for momentum 
s mixing coefficient for sediment 


coefficient in logarithmic velocity distribution. 
kinematic viscosity of water, ft2/sec 


€ 
K 
V 
p density of water, y/g, lb sec2 /ft4 
Bs 
7 


intensity of shear, Ib /ft2 
c Critical tractive force for beginning of sediment movement 
Tg boundary shear, or tractive force, at stream bed, VyoS 


T> boundary shear associated with sediment particles 
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SYNOPSIS 


Lake Bloomington, an impounding reservoir, has been subjected to detailed 
jurveys in 1948, 1952 and 1955 to determine the deposition of sediment. Dur- 
ng each of these surveys samples of the sediment were obtained. Particle 
ize distribution analyses of 30 of these sediment samples were utilized to 
letermine the tons of sediment deposited in the lake during each of these 
edimentation periods. Postulating that sediment particles which had a di- 
meter greater than 50 microns had been moved into the lake as bed material 
oad, the total tons of this sized material was calculated based on the sedi- 
nent samples. 

An hydraulic study was made of the 2-1/2 mile reach of Money Creek im- 
nediately upstream from Lake Bloomington to determine its sediment-carry- 
ng capacity. A series of sediment samples were taken of the bed material of 
he Money Creek channel. Utilizing these data, curves of water discharge 
ersus sediment discharge were computed utilizing three different methods: 
he Einstein procedure, the Schoklitsch formula, and the DuBoys formula. 

At the lower end of this stream reach, immediately upstream from Lake 
sloomington is a stream-gaging station for which flow records were available 
or each of the three sedimentation periods. Utilizing flow duration infor- 
aation from this stream gage, the total quantity of bed material moved 
hrough the Money Creek reach was calculated utilizing each of the three sedi- 
1ent transport relationships developed. The actual bed material-size sedi- 
1ent in Lake Bloomington is compared with the sediment transport as com- 
uted by the three methods. 


fote: Discussion open until July 1, 1958. A postponement of this closing date can be 
obtained by writing to the ASCE Manager of Technical Publications. Paper 1531 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 1, February, 1958. 

. Associate Engr., Illinois State Water Survey Div., Urbana, Ill. 

. Formerly Asst. Engr., Illinois State Water Survey Div., Urbana, Ill. 

. Formerly Asst. Engr., Illinois State Water Survey Div., Urbana, Ill. 
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INTRODUCTION . 


Purpose 


The purposes of this investigation are: (1) to attempt to compute the sedi 
ment moved as bed material load in Money Creek, Illinois, for a reach im- 
mediately upstream from a stream-gaging station and from Lake Bloomi . 
(2) to provide a comparison of three well-known bed-load formulas including : 
the most recent one proposed by H. A. Einstein(1) and (3) to compare the re- 
sults of these formulas with bed material measured by actual survey in Lake — 
Bloomington. 

Most sediment transport formulas are derived from laboratory flume 
studies. The importance of this study is believed to be the check of three 
formulas under natural conditions. 


General 


The city of Bloomington is located in the central part of McLean County, 
Illinois. A public water supply derived from wells was installed for the town ~ 
in 1875. The wells were utilized until 1929 when Lake Bloomington was 
formed by the construction of an earth dam across Money Creek, a tributary 
to the Mackinaw River about 15 miles northeast of Bloomington. Since that 
time Lake Bloomington has been used for the public water supply. The lake 
has a drainage area of 61.0 square miles, a surface area of 487 acres, and 
had an original storage capacity of 2.17 billion gallons. Fig. 1 shows the lo- 
cation of Lake Bloomington and Money Creek watersheds. 

A detailed survey was conducted on Lake Bloomington in August 1948 to 
determine the volume of sediment deposition. At that time a series of 19 
cross sections of water depth and sediment thickness were taken on the lake 
as shown in Fig. 2. In August 1952 and in July 1955 soundings were repeated 
along these same cross sections to measure the further sediment deposition 
during the intervening periods. These cross sections were plotted and the 
total quantity of sediment (including the bed-load and wash-load) deposited in 
the lake was calculated by the method devised by the Soil Conservation 
Service. ‘2 

Table 1 is a summary of the results of the three surveys. It will be seen 
that the capacity of the lake has been depleted every year at an average rate 
of about 0.46 percent of the original capacity. The 1955 survey showed a tot: 
loss in capacity of 791 acre feet or 258 million gallons. 

During a period of low inflow and low lake level during 1954 less sediment 
was carried into the lake and a portion of the deposited sediment bed in the 
upper portion of the lake was exposed to air drying and consequently compact= 
ed. The specific weight of the sediment deposit in each segment of the lake 
was determined by a series of sediment samples obtained during each of the 
three surveys. A total of 30 samples were utilized to determine the tons of 
sediment deposited in the lake during each of the sedimentation periods. 


Choice of Reach 


Lake Bloomington and its watershed have been the subject of an hydrologic 
study by the State Water Survey since 1933,(3) The principal tributary to 
Lake Bloomington is Money Creek. Immediately upstream from the head- 
waters of the lake is located a stream-gaging station. This gage is sponsor 
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Figure 1 - Location of Money €reek and Lake Bloomington 
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Figure 2 - Sedimentation Survey of Lake Bloomington 
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Table 1, Summary oi: Sedimentation Data 


Take Bloomington, Tllinols 


AGE 
Dec. 1929 - Aug. 1948 - 18.7 yrs. 
Aug. 1945 - Aug. 1952 - 4.0 yrs. 
Aug. 1952 - July 1955 - 2.9 yrs. 
Dec. 1929 - July 1955 - 25.6 yrs. 
WATERSHED 


Total area - 61.9 sq. miles 
39,00 acres 

Land area - 60.2 sq. miles 
38,553 acres 


RESERVOIR 


Area at spillway crest - 87.2 acres 


1929 1948 1952 1955 __Units _ 
Ape a ae ome Dee tat 
cope catiabe ance 109 99 97 96 Acre-feet 
SEDIMENTATION 1929-1948 1948-1952 1952-1955 1929-1955 
Total 592 US7 2 791 Acre-feet 
Averaze Annual Accumulation 
Sotecshed 3087 39.3 142 30.8 acre-feet 
Per sq.mile! 0.53 0.65 0.234 0.51 Acre-feet 
Per acrel S6eu h.2 1972 34.7 Cubic-foot 
Tons per acre 0.7h 0.91 0.34 0.72 Tons 


DEPLETION OF STORAGE 


Loss of orig- 


inal capacity 2 
Total period 8.90 2.36 0.635 11.89 Per cent 
Per Year 0.48 0.59 O21 0.46 Per cent 


diana area only. 


-Volune compacted due to drying. 
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by the State Water Survey and is operated by the United States Geological 
Survey. It records the drainage from 51.9 square miles of the total lake ' 
watershed. Records are available at this station from June 1933 to date. Be-_ 
cause of the availability of these discharge records, the stream reachim- ~ 
mediately upstream from this gage was given consideration for the present — 
study. ( 

In selecting a river reach for sediment transport calculations, it must be — 
remembered that such a function can be applied only to a river reach of uni-~ 
form flow. This means that the length of the channel must be sufficient to 1 
permit adequate determination of the over-all slope. Also the channel itself — 
must be sufficiently uniform in shape, sediment composition, slope and out- 
side effects such as bank vegetation, that it is possible to treat the reach as a 
uniform channel characterized by an over-all slope and by an average repre- 
sentative cross section. Practically, it is difficult to realize such an ideally 
uniform channel. After a field inspection of Money Creek, however, it was ~ 
decided that this reach was sufficiently uniform to be utilized for bed-load 
calculations. 


Hydraulic and Hydrologic Determinations 


Hydraulic Properties of Channel 


A field investigation was made of the 2-1/2 mile reach immediately up- 
stream from the gaging station to determine hydraulic properties. A series 
of 13 cross sections was taken of the stream at approximately 1000-foot inter- 
vals. The slope of the water surface was determined to be 0.000905 by utiliz- 
ation of these 13 stations. The slope measurement was taken during a time 
when the average discharge was 160 cubic feet per second. Flow duration 
studies showed discharge equalled or exceeded this amount for six percent of 
the period of record. 

To determine the stage-area relationship for the reach, each of the 13 
cross sections was plotted in actual position in elevation. The average cross 
section for the reach was then determined by sliding all cross sections down 
the channel along the slope 0.000905 into the plane of the section at the lower 
end of the reach. This downstream cross section was at the stream-gaging 
station, giving a means of comparison between the mean cross section for the 
reach and the actual cross section of the gaging station. The stage-area 
curve for the reach determined in this manner is shown in Fig. 3. 

In a similar manner the average stage-wetted perimeter curve was ob- 
tained and is shown in Fig. 4. In a wide channel like Money Creek the width © 
of the channel was considered as the wetted perimeter for a known elevation. — 
Consequently, corresponding values of the two curves made it possible to 
compute the hydraulic radius for the same elevation. These are shown in F 
Fig. 5. The curve of stage versus hydraulic radius with bank friction was 
computed by means described by Einstein, (1) The stage-discharge relation- © 
ship with and without bank friction is shown in Fig. 6 as compared to the actu- 
al stage-discharge relationship as measured at the gaging station. 

In order to determine the width of the stream bed along which transpor- ) 
tation takes place, the average widths of the stream at various stages of all 
of the 13 cross sections wre plotted and are shown in Fig. 7. As reported by 
Chang(4) and by Einstein 5) the movement of bed material is reasonably ex- — 
pected to take place only along the bed portion of the channel. From Fig. 7a | 
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Figure 3 - Stage-Area Relationship for Money Creek 
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Figure | - Stage-Wetted Perimeter Relationship for Money Creek 
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Figure 5 - Stage-Hydraulic Radius Relationship for Money Creek 
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Figure 6 - Stage Discharge Relationship for Money Creek 
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width of 30 feet was arbitrarily chosen to represent the width of bed for the 
reach of Money Creek under consideration. 

Computations of hydraulic properties need to be made only up to a stage 
-orresponding to the highest flood that had occurred in this creek. This stage 
of seven feet was the maximum utilized in determining all hydraulic proper- 
‘ies. 


Bed-Sediment Samples 


Since the bed-load formulas used in this study relate the grain-size compo- 
sition of the bed material with the flow of the channel it is necessary to ob- 
ain representative samples of the bed material. A total of 18 samples of 
sediment was taken along the active channel and the flood plain of Money 
creek. Samples were obtained by means of an auger or a pipe sampler and 
were taken to a depth of about 1.5 feet, the estimated depth of scour or active 
Jed movement. The flood plain samples indicated 90 percent by weight to be 
finer than 50 microns, and it was concluded that these finer particles were 
leposited during the recession of flood flows. Consequently the flood plain 
samples were abandoned and five samples were chosen to represent the bed 
material in the active channel. The size distribution of this bed material 
Jased on these samples is shown in Table 2. 

The data from Table 2 are plotted in Fig. 8, from which can be noted the 
sharacteristic grain sizes. The size which enters the Einstein equation of 
Transport is Dgs = 0.195 mm = 0.000639 feet (grain size of which 35 percent 
s finer). The size characteristic for friction Dg5 = 1.22 mm = 0.0040 feet. 

The sediment transport was calculated for grain sizes between 9.4 milli- 
neters and 0.050 millimeters which represents 67 percent of the bed materi- 
ul. It is important to recognize however that 12 percent of the bed material — 
Ss coarser than 9.4 millimeters. These gravel-size particles do not move for 
iormal flows although an extremely small rate of transport may occur at high 
lood stages. Twenty-one percent of the bed material is finer than 0.050 
nillimeters. As much as 15 percent of this size may be expected to be a part 
yf wash-load particles lodged behind the coarser grains. This can generally 
ye neglected. No adjustments were made since a large percent of these finer 
naterials were found in the bed. Though included in the bed material size 
listribution, no bed-load function exists for these finer particles. Calculations 
yere made for individual sieve fractions using as representative the average 
‘(rain sizes varying from 7.3 to 0.073 millimeters as shown in Table 2. 


Jalculation of Bed Material Deposited in the Lake 


One of the most critical phases of the present investigation was the de- 
ermination of the total quantity of sediment deposited in Lake Bloomington 
vyhich was of the size moved as bed-material-load through the tributary 
reek, Money Creek. The quantities of sediment measured in the lake survey 
ontained principally fine material which was undoubtedly moved into the lake 
is wash load. The 16 sediment samples taken during the 1948 and 1955 lake 
urveys were considered to be sufficient in number to indicate the sediment 
ature in the various segments of the reservoir. These samples were sub- 
ected to size distribution analyses. The locations of these sediment samples 
re shown in Fig. 2. 

All particles of sizes 0.050 millimeters (50 microns) and more were con- 
idered as bed-load and particles of sizes smaller than 0.050 millimeters, in 
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Figure 7 - Money Creek Average Cross Section 
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Figure 8 - Average Size Distribution of Bed Material Samples 
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Table 2. Average Size Distribution of Money Creek 
Bed Material Based on Five Samples 
Grain Size Per Cent Mean Diameter 
mm by Weight mn Feet 
Less than 0.050 cle Oe Beta Gssetercs 
0.050 to 0.097 4.0 0.073 0,000239 
0.097 to 0.19 9.5 OvLiiee O.00072 
0.19 to 0.37 10.7 0.28 0.000918 
0.37 to 0.72 153 0.545 0.00179 
OpFPavtockelh: TOT 1.06 0.00348 
AIT eye 8.8 2.05 0.00672 
Sy te sie 7.0 3.95 0.01295 
Sse) tetoy Deilh 5.0 703 0.02393 
Greater than 9.) 12.0 we weeenae 


Not 
tp 


Table 3. 
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Sand Content of Lake Sediment Samples 
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the range of silt and clay, were considered wash load. Table 3 presents the 
results of these analyses, showing for each sample, the total percent by ; 
weight of material having a diameter greater than 50 microns. In only two of 
these samples was there material having a diameter in excess of 9.4 milli- 
meters and in each of these cases the percentage was extremely small. 

The sediment samples at the 16 locations in the lake were utilized to com= 
pute the total quantity of sediment in tons in each segment of the lake. The — 
locations of the sediment samples and the lake segments are shown in Fig. 2. 
Table 4 shows the results of these calculations. } 

Reference to Figs. 1 and 2 shows that Money Creek is tributary to that 
portion of the lake containing segments 5 through 15 and Hickory Creek con-— 
tributes the flow to segments 16 through 21. The two arms of the lake come 4 
together to form that portion of the lake containing segments 1, 2, 3 and 4. H 
The tonnage of sediment as calculated in Table 4 for these four segments (1 
through 4) was calculated by proportioning the total tonnage coming from t 
Hickory Creek and from Money Creek. The drainage area of the Money Creek 
above the lake is 51.9 square miles and the drainage area on Hickory Creek 
above the lake is 10 square miles. The proportion of these two drainage : 
areas was used to calculate the tonnage of bed material sediment deposited i! 
in these lower four segments of the reservoir. It will be noted from Table 4) 
that the total quantity of bed material, which can be assumed to have moved | 
down Money Creek and into. Lake Bloomington during the period 1929-1955, fi 
amounts to 60,527 tons. i 


Flow Duration Data 


Computation of the bed load by any of-the three formulas developed a re- : 
lationship between sediment discharge in tons per day and water discharge in | 
cubic feet per second. To determine the total quantity of material moved — 
through this reach of the stream as bed load, it was necessary, therefore, J 
construct flow duration curves. It was desired that such data be available fol 
each of the three periods during which sediment deposition was measured in © 
Lake Bloomington. Water discharge data from the stream-gaging station on ~ 
Money Creek were complete in this respect except for the period 1929 to 1933 
and for the year 1941. To complete the flow record, the discharge at Money 
Creek was synthesized for these missing years. 

Stream-gaging records were available for the neighboring Mackinaw River 
for the missing periods of record as well as for the complete period of recoré 
for Money Creek. The two drainage basins were assumed to be homogeneous 
in regard to their general flow characteristics. Flow duration curves for the 
two rivers were drawn based on the same period of records, namely 1934 
through 1940 and 1942 through 1954. These were used to synthesize the mis 
ing records at Money Creek based upon the actual flow measurements on the 
Mackinaw River using a method described by Mitchell.(6) The actual flow 
records of Money Creek which is in the Mackinaw River basin have been 
published by the United States Geological Survey in their Water Supply 
Papers. \7 

To determine the bed-load quantity, the flow duration data for each of the 
three different periods were compiled as described above and shown in Tabl 
5. Only flows above 110 cubic feet per second have been considered. This 
assumes that the quantity of bed material moved by lesser flow is negligible, 
Table 5 shows the duration of the high flows which have occurred at the Mon 
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Table 4. Total Tonnage of Sand and Larger Material 
Deposited in Lake Bloomington 
Total Sand 
Segment Sediment Samples Per Cent Tons 
Tons Sand Hickory Cr. Money cr. 
(>50 microns) Arm of Lake Arm of Lake 

1 10,881 12a,12b 2.75 Sh 245 

2 32, a0 12a,12b 2.75 163 737 

3 46,828 vt veo 85 383 

4 h2, 238 10,11 U.75 134 605 

5 26 5,157 10 2.5 -- 703 

6 33,945 9a,9b 4.5 -- 1528 

7 49,29) 9a,9b 4S -- 2218 

8 41, 323 8 0.5 -- 207 

9 56,2h5 7a,7o,8 0.5 -- 281 
10 71,130 6 211 os 8594 
11 99,630 5 RT -- 16,937 
12 4138 364 27.5 = 11,506 
1 27,527 2,3 37-5 -- 10, 323 
1h 5,061 1,2 1 -- 3796 
15 2,541 1 97 -- 26h 
16 45,285 16 1.0 453 See 
1 16,407 16 220 16L ie 
1 11,960 15a,15b 2.75 329 --- 
19 15,182 1h,15a,15b 2.37 329 ee 
20 16,477 13,14 10.3 1697 py 
21 10,027 13 18 1805 — 

60,5 27 Tons 
Table 5. High Flow Duration During Lake Sedimentation Periods 


Money Creek Gaging Station 


Mean 


Discharge Duration in Days 
f ist Period ond Period 3rd Period Total 


cfs 


Dec. 1929 Sept. 1948 uae 1952 


to to fo} 
Aug. 1948 Aug. 1952 July 1955 


6 i 59 559 
ig9 24 20 277 
10 34 has 148 

5 a 6 2 
31 10 3 yh 
11 2 -- 13 

3 1 cs 4 
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Creek gaging station during the three periods for which sedimentation was 

measured in Lake Bloomington. Discharges considered are mean values of 
incremental discharge ranging from 110 cubic feet per second to 1100 cubic 
feet per second. 


Sediment Transport Calculations 


The Einstein Procedure 


H. A. Einstein developed and published in 1950 a complex procedure for 
computing the quantity of bed material transported by a stream.(1) This b 
load function was applicable to an alluvial channel in an equilibrium state, 
which moves the material through which it flows. 

The approach was based on the probability of movement of a particle of 
particular diameter in the “bed layer.” Movement in this layer was con- 
sidered to occur by rolling and sliding on the bed or by making a series of 
short hops and was termed “surface creep.” The thickness of the bed layer 
was postulated to be twice the grain diameter. The movement of particles 
was considered to be governed by statistical laws of probability and was so 
related to the flow. The average distance traveled by a particle betweenpe 
ods of deposition was assumed to be 100 diameters. 

The concentration of particles having a particular diameter at the top of 
the bed layer is assumed to be equal to the concentration of suspended parti- 
cles of the same diameter at this same boundary. This concentration is then — 
related to the concentration of similar particles at any elevation in the verti- 
cal. By integration of the function the total sediment load of this diameter, 
per unit width, was determined at a representative vertical in the stream 
cross section at a given discharge. Load was calculated for a number of 
grain-size categories based on the samples of bed material. In calculating 
the total load of a mixture of particles, corrections were introduced for the 
“hiding factor” or interference of the larger grains with the smaller. A later 
publication by Einstein\8) improves this correction. . 


The principal relationships utilized by Einstein in the bed load function 2 
as follows: 


d 
qs7f Cy dy dy 


HOEY ad 
C, =, Veda 


Vs 


Z +d 
Uy 
Where, 


4s = Rate of transportation of suspended load 
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d = Water depth 
Uy = Velocity at distance y above bed 


Cy = Sediment concentration, weight per unit volume, of the fluid-sedi- 
ment mixture at distance y above the bed 


C, = Sediment concentration at distance a above bed 
Vg = Settling velocity of a sediment grain in still water 
k = Universal constant of von Karman 

» = Shear velocity at the bed 


Einstein‘) considered the velocity distribution in open-channel flow over a 
diment bed as being best described by the logarithmic formula based on 

n Karman’s similarity theorem with the constants as proposed by 

ulegan, 9) He gave the vertical velocity distribution including the transition 
tween the rough and smooth boundaries as: 


u 


reg SI Few LO Srey 


(4) 


= 5.75 09,9 (30.2 — ) 


erein x is given as a function of kg/g 
Uy = the average point velocity at the distance y from the bed 


ul = V To / S¢ = V Se -R°- g (the shear velocity) (5) 


sgt =the density of the water 
Se =the slope of the energy grade line 
= the hydraulic radius 


R 
g = the acceleration due to gravity 
y =the distance from the bed 

k 


= the roughness of the bed 


x =a correction parameter 
A= k Ss coax the apparent roughness of the surface (6) 
| 1.6 ; the thickness of the laminar sublayer 
d= of a smooth wall (7) 
Uu 


v =the inneraatic viscosity of the water 
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The value of k, for uniform sediment equals the grain diameter as de- 
termined by sieving. The representative grain diameter of a sediment mix- 
ture is given by that sieve size of which 65 percent of the mixture by weight 
is finer. 

The total rate of sediment transportation is the sum of the suspended and 
bed-load transport rate and is given by Einstein in Equation (63) of Referenct 


(1) as, 
= in dp( PI, + I, + |) 


where: 
iy = Fraction of total load in a given size range 


ay = Total transport rate, weight per unit time and width 


ip = Fraction of bed load in a given size range i 
ap = Bed load transport rate 
Ij = Integral value (Evaluation tables furnished by author) iy 
I, = Integral value (Evaluation tables furnished by author) } 
P = Parameter of total transport ; 
and, 
} 


| = (SORES | 
0.434 'Si0\kK./d o) 


P= 


considered in this paper, the hydraulic character of the channel was computed 
in accord with the methods reported by Einstein. Detailed computations are 
not presented here but the resultant effects of bank friction are shown in Figs. 
4 to 6 of this report. The computations of sediment movement are based on 
these hydraulic computations including the bank friction. . 

The relationship of sediment discharge to water discharge for Money 
Creek as determined by the Einstein procedure is tabulated in Table 6 and is 
shown graphically in Fig. 11. Table 6 shows the utilization of this relation _ 
and the flow duration data to determine the total bed material movement into 
Lake Bloomington during the sedimentation periods under consideration. 
Total transport calculated by this means amounts to 196,477 tons. 


Schoklitsch Bed-Load Formula 


The sediment movement in Money Creek was calculated by utilization of 
the Schoklitsch formula for uniform sand, 


feo O..f, 1.5 
G ” Aare ss) Bq: —qyyam . (19) 
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~O.000 52d 


do = (11) 
S PS i 


id the bed load for a mixture 
Gy =aGq + bG, + CGc ---+mGm (12) 


here 
G = Bed load in tons per day 
Gt = Total bed load for a mixture of particles 
Ga = Quantity of bed load of a particular diameter 


a = Percent weight of a particular diameter in a mixture 
m = Number of size-gradation divisions in a mixture 
d = Diameter of particle, inches 


S = Hydraulic slope 
B- = Bed width, feet 
q = Discharge, cfs 


Qo = Critical discharge at which movement of particle of diameter d, 
begins 


The Schoklitsch(10) formula serves to compute that portion of the total load 
‘solids in the river which is transported (not in suspension) along the river 
sd by the tractive force of the stream. The Schoklitsch formula is based 
ainly on the classic flume experiments of G. K. Gilbert besides additional 
perimental data collected by Schoklitsch. It was developed for uniform 
‘ain material but there can be no valid objection to its being applied to mix- 
res as well. It has been verified and found to agree closely with the 
easurements in the River Danube and the Terek River. 

In the application of this formula to a natural stream it was stated by 
hoklitsch that the reach studied be relatively straight and the depths of 
ater as uniform as possible in order that the width of the stream change as 
tle as possible with stage. Table 7 shows the relation of area and width of 
oney Creek at the various discharges considered. 

Table 8 summarizes the movement of bed-load material in Money Creek 
calculated by the Schoklitsch formula. The relationship of sediment dis- 
arge to water discharge is plotted in Fig. 11. Table 8 shows the product of 
diment discharge rating and flow duration information converted into total 
antity of bed-load movement in tons for the various periods. Total 
ansport by this method amounts to 79,065 tons. 


iBoys Formula 


The sediment transport in Money Creek has been calculated by the DuBoys 
rmula. 


Gea 6 eels (= — "7 0) (13) 
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Table 6, Money Creek Sediment Discharge by 
Einstein Bed Load Function 


For particles 0.05 mm to 9. mm 


Water lst Period 2nd Period rd Period 

Dis- Sediment Dura- Sédi- Dura- Sedi- a- = Sedi- 

charge Discharge tion ment tion ment tion ment ment 
cfs Tons/Day Days Tons Days Tons Days Tons Tons 
110 94 376 35,344 12h 11,656 59 5,546 2,546 
165 145 199 28,855 58 8,410 20 2,900 40,165 
2h5 232 103 23,896 3h 7,888 11 epee 3, 336 
355 335 58 19,430 18 6,030 6 2,010 27,470 
520 550 31 17,050 10 5,500 3 1,650 24,200 
760 920 ll 10,120 2 1,80 “= ewer 11,960 
1100 1450 3 4, 350 1 1,450 -- ----- 5,800 

Total 139,045 42,774 14,658 196,477 


Table 7. Money Creek Area-Width Relationship 


Discharge Stage Area Width 
SoA Ss Feet Sq.Ft. Feet 
110 1.98 49.5 25.0 
165 2.38 63.0 26.5 
2h5 2.91 83.0 28.5 
355 3-49 106.0 30.4 
520 4.13 133.5 32.3 
760 4.8 171.0 35.3 


1100 5263 226.5 0.2 
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d the bed load for a mixture 


G = (aqqg+ bqy+ cag+ --- nan) YS, us) 


T = vys (15) 
1ere 


q = Transport rate of a particular diameter particle in volume per 
second per foot of width 


Cg = Sediment parameter 
= Intensity of bed shear 
= Unit weight of water 


T 
? 

y = Depth of flow 
s = Hydraulic slope 
Te 


= Value of 7 for which qg is zero 
Ss = Specific gravity of sediment particle 
a = Percent weight of a particle diameter in a mixture 
n = Number of size-gradation divisions in a mixture 
G _ =Bed load total for mixture, pounds per second per foot of width . 


The DuBoys formula was one of the earliest published to determine the 

‘iid transport of sediment. A great number of other mein have been de- 
loped subsequently and have a similar nature. Johnson 11) tested a number 
these and concluded that all formulas fitted equally well, thus indicating 

at the choice could be made on the basis of convenience. In order to utilize 
s formula, it was necessary to evaluate the parameters Cg and Tc. The 
lues summarized by Straub and published in Engineering Hydraulics(12) 

re utilized. It was necessary, however, to extrapolate these relationships 
shown in Fig. 9 for the relation of Cg, to particle diameter, and in Fig. 10 
: the relation of J, to particle diameter. 

In Table 9 is summarized the results of the calculation of bed-load 
yvement in Money Creek by the DuBoys formula. The relation of sediment 
s;charge to water discharge has been plotted in Fig. 11. Table 9 shows the 
mputation of the total sediment movement throughout this reach based on 

» flow duration of Money Creek for the three sedimentation periods. It will 
noted that the total quantity of sediment moved calculated by this means 
jounts to 529,944 tons. 


Discussion of Results 


neral 


In Fig. 11 is shown the sediment discharge versus water discharge for 
ney Creek as determined by the three different methods. Table 10 shows 
comparison of the quantity of sediment measured in Lake Bloomington 


1531-20 HY 1 February, 1 
Table 86. Money Creek Sediment Discharge by 
Schoklitsch Bed Load Formula 
For particles 0.05 mm to 9.) mm 
lst Period 2nd Period ged Period 
Water Sediment ura= edi- an edi- ura= edi- 
Discharge Discharge tion ment tion ment tion ment 
cfs Tons /Day Days Tons Days Tons Days Tons 
110 33.6 376 12,622 12h 4,163 59 1981 
165 59.8 199 11,894 58 3,467 20 1195 
2k5 98.0 103 10,090 3k 3,331 11 1078 
355 152 58 8,797 18 2,730 6 910 
520 232 SH 7,492 10 2,320 696 
760 39 11 3,8h2 2 698 a —— 
1100 515 3 1,544 a 515 a Bok 
Total 55,981 17,224 5890 
Table 9. Money Creek Sediment Discharge by 
DuBoys Bed Load Formula 
For particles 0.05 mm to 9.4 mm 
lst Period 2nd Period rd Period 
Water Sediment ura = Sedi- Dura- Sedi- Siege Boar. 
Discharge Discharge tion ment tion ment tion ment 
» efs Tons/Day Days Tons Days Tons Days Tons 
110 2b2 376 105,855 12h 34,910 59 16,610 
165 08 199 81,180 58 23,661 20 8,159 
245 657 103 67,647 3h 22,330 11 7,22h 
355 921 58 53,402 18 16,573 6 5,52h 
520 1295 31 40,130 10 12,945 3 3,884 
760 1679 11 18,65 2 3,357 Soh Bi ea ar 
1100 2022 3 6,066 a 2,022 -—- wwe 
Total 372, 745 115,798 41,401 529, 944 
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with the amount moved through the Money Creek reach, as computed by eac 
of the three methods. 

As shown in Table 10, the Schoklitsch formula gives results most nearly ~ 
in agreement with the actual sediment measured in Lake Bloomington. For 
the total period of record, this formula gives results in tons only 31 percent 
greater than the measured quantity. In comparison the Einstein procedure 
gives results 225 percent too great, and the DuBoys formula, 776 percent too” 
great. 

The usefulness and limitations of the three methods utilized in this paper 
to compute sediment movement have been discussed in detail by Chien, (13) 
Recognizing the limitations, these approaches merit continued study, trial, 
and improvement. These approaches were utilized on the Niobrara River 
near Cody, Nebraska 14) in 1955. Results showed severe limitations to the — 
Schoklitsch and DuBoys approaches and excellent results from the Einstein 
approach, 

One recognized source of error in the use of a bed-load formula is the us¢ 
of the water surface slope instead of the slope of the energy gradient. An ac 
curate determination of the slope of the energy gradient requires the measure 
ment of the velocity distribution at each end of the experimental reach. This 
observation is often eliminated and the resulting error involved in the slope” 
determination is fairly large for the usual experimental conditions. It is 0 
interest to note that Gilbert{15) was undecided as to the proper value of slog 
to use and stated “I do not find it easy to decide which slope should be regz 
ed as the true correlative of capacity for traction but as all of our laborator 
data include the debris slope while the determinations of water slope are 
relatively infrequent the discussion of the results has adhered almost ex- 
clusively to the former. If the water slope is the true correlative then the 
use of the debris slope involves a systematic error.” 

Professor O’Brien and Lt. B. D. Rindlaub{16) support Gilbert’s selection 
in the statement “It is to be noted that the slope at the bottom is more near]j 
equal to the slope of the energy gradient than is the slope at the water sur x 
and partly for this reason the data of G. K. Gilbert show less scattering thai 
the data of more recent experimenters who have criticized Gilbert for not — 
measuring the slope of the water surface in all of his experiments.” 


The Einstein Method 


Sediment movement and river behavior are inherently complex since natt 
ral phenomena involve a great many variables. In applying the unified methe 
presented by Einstein, 1) questions may arise such as: (1) Is it possible to 
obtain a truly representative bed material sample for size distribution 
characteristic curves? (2) What is the average (representative) diameter ¢ 
the entire sediment mixture? (3) Is it possible to select an ideally uniform 
channel in nature? (4) Can the formulas for the hydraulics of the open 
channel be applied to such a complex problem as that of sediment transport 
phenomena? (5) If only one point of the size distribution curve is to be use 
for roughness height kg, (Dg5), how much confidence can the engineer have i 
its practical use? (6) Is the evaluation of bar resistance accurate and suf- 
ficient? (7) What is the lower limit in integrating suspended load? Many 
other questions may also arise. However, the engineer is warned against bi 
ing discouraged by the absence of a better solution. In answer to the above 
questions it should be pointed out that the available information on the subje 
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able 10, Bed Load Material Deposited in Lake Bloomington 


Compared to Computed Bed Material Movement in 


Money Creek, 


Tons 
jiment : s erio n erio r er 
asured 
9osited in 
1ey Creek Arm 
Lake Bloomington 48,176 10,725 626 


nputed 
(Per cent error shown below each value in tons) 


istein iets 2,774 1,658 
159 299 802 
10klitsch 55,981 Bly payer lh 5,890 
16 61 262 


30ys 372,745 115,798 41,401 
67h. 980 2,446 


? 


60,527 


196,477 
225 
79,065 
31 


529,944 
776 
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of hydraulics of open channel flow has been applied as closely as possible. 
Other information is rationalized through practical experience and field 
measurements. 


Schoklitsch and DuBoys Formulas 


Although these two formulas are essentially the same in structure, thei 
application to Money Creek gives results at great variance. It will be see 
that the results of the Schoklitsch formula are more nearly in agreement with 
the survey data. One reason for the high values in the case of DuBoys may ; 
be the limitations to the evaluations of the parameters Cg and T;,. Since the 
sediment parameter C, expresses the relative susceptibility of a givens 
ment to movement and since the shear terms Tand7T ¢ involve the compl 
system of forces exerted by the flow upon the bed, the evaluation of these 
parameters by means of suitable experimental methods determines the reli: 
ability of the results. The adaptation of the values, utilized as shown in Figs. 
9 and 10, to conditions in Money Creek is therefore somewhat questionable, — 


The Manning formula permits DuBoys relationship to be written in the 4 
following alternative form. 


i 
q =CsT (T-Te) (16 


re g |-4 


G pie AS De 
$(1.49/n)/-2 


The exponent of the slope is 1.4 in this relationship as against 1.5 in the ; 


(q L.2ac8 q ol 


5 
Schoklitsch formula while the exponent of qg is 1.2 as compared with 1.0. A 
though this shows general agreement the disparity between the results of 
two formulas is well accounted for. In addition to this, it is a recognized 
fact(12) that the DuBoys formula was based on an incorrect assumption as to 
the sliding motion of the sediment particles in movement. 


CONCLUSION 


By utilizing the measured quantity of sediment of bed-material size in 
Lake Bloomington as a check on the computed sediment which has moved 
through Money Creek it is concluded that the Schoklitsch formula gives the 
most reliable results, being only 31 percent in error for the over-all period 
of the study. The Einstein procedure gave results 225 percent high and th 
classic DuBoys formula gave results 776 percent high. 
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ABSTRACT 


The Corps of Engineers hydraulic field testing program is coordinated by 
e U. S. Army Engineer Waterways Experiment Station, Vicksburg, Missis- 
opi. Included are flood-control and multiple-purpose projects, river and 
rbor works, estuarine and wave problems. Test data discussed include 


illway, conduit-intake, and gate-leaf pressures; air demand; and concrete 
nduit friction factors. 


INTRODUCTION 


This paper is limited to hydraulic field investigations of spillways and out- 
- works constructed and operated by the Corps of Engineers, particularly 
9se in which the U. S. Army Engineer Waterways Experiment Station, locat- 
at Vicksburg, Mississippi, has played a part. We call the structure where 
ld tests are made “the prototype” as opposed to a working model on which 
sts are made at a smaller scale. Usually more variables are studied in the 
Jdel tests than in field tests because of the ease in control and lesser time 
d costs involved. 

The main purpose of field tests is to obtain data for use in producing more 
onomical designs. Frequently data are needed for conditions where model 
idies are not applicable, or the scale ratios are unknown. For example, 
nulating air-water mixtures is difficult in a model, and the scale ratios for 
vitation and those involving high Reynolds numbers are not clearly under- 
yod. Another purpose of field tests is to investigate and evaluate the 


te: Discussion open until July 1, 1958. A postponement of this closing date can be 
obtained by writing to the ASCE Manager of Technical Publications. Paper 1532 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 1, February, 1958. 
Presented at Jackson Convention of the American Society of Civil Engi- 
neers, February 18-21, 1957. 

Chief, Prototype Section, Hydraulic Analysis Branch, Hydraulics Division, 
U. S. Army Engineer Waterways Experiment Station, Vicksburg, Miss. 


1532-1 


1532-2 HY 1 February, 1958 


performance of the prototype, usually for some phase of performance in w 
operating difficulties are occurring. Also tests have been made to verify . 
model indications, and in some instances field tests have been necessary to 
obtain data for designing or operating the model. : 

Within the Corps of Engineers field tests are made either by the District — 
operating the project or by the Waterways Experiment Station in conjunction ~ 
with District personnel. District personnel make tests involving convention 
equipment such as current meters and stage recorders, take discharge : 
measurements for rating gates, and make river and harbor soundings. Som 
extensive tests involving pressure measurements, vibration, and hoist force 
are also made by Districts. However, most of the extensive tests requiring 
the use of special equipment are made by the Waterways Experiment Station 
and Districts cooperating. 

The Experiment Station is the coordinating agency for hydraulic prototy 
testing in the Corps of Engineers. In this capacity recommendations are ~ 
made to the Districts concerned regarding collection of field data and insta 
lation of testing facilities such as embedment of piezometers or conduit for 
electric leads during construction. The Waterways Experiment Station is in 
a favorable position to carry out this assignment because of early associatio 
with new projects—sometimes through model tests during which the need for 
field tests may become apparent. Also the Experiment Station in dischargin 
its responsibility of disseminating hydraulic design criteria to the Corps 
frequently sees the need of field tests. The presence of an Instrumentation 
Branch with a reservoir of electronic measuring equipment and instrument- 
ation engineers adds much to the capabilities of the Station. 

The scope of the Corps’ hydraulic field testing program includes flood- 
control and multipurpose projects, river and harbor works, and estuary and 
wave-action problems. In recent years the program has been concentrated ¢ 
multipurpose projects, particularly the spillways and outlet works of concre 
and earth dams. The spillway tests have been primarily confined to the 
measurement of pressures and the collection of data for determining spillwa 
discharge coefficients. Tests on outlet works have been of broader coverag 


Spillway Tests 


In addition to the problems involved in the installation of facilities during 
construction, the problem of water availability for testing is more acute for 
spillways than for outlets. Obviously spillways cannot be operated until the 
reservoir level rises above the spillway crest. In many cases this does not 
happen for years. Some uncontrolled spillways on projects over 15 years ol 
have not yet gone into operation. A severe limitation in testing controlled — 
spillways is that large flows cannot always be tolerated downstream. This 
restriction hampered testing at Pine Flat Dam in June 1956. Nevertheless, : 
limited number of spillway tests have been made. This include velocity 
measurements in the chute of Fort Peck Dam in 1946; ogee crest pressure _ 
measurements at 19 per cent of the design head and the taking of water- 
surface profiles at Arkabutla in 1953; recording of crest and baffle pier 
pressures at 75 per cent of the design head at McNary2 in 1954; and record: 
ing of crest, tangent, and flip-bucket pressures at Pine Flat at 57 per cent of 


2. Walla Walla District, Corps of Engineers, Model-prototype Conformity, 
McNary Dam Spillway. June 1956. 


SCE FIELD INVESTIGATIONS 1532-3 


e design head for gate openings up to 8 ft in 1956. 

The most extensive tests yet made at a Corps spillway were at Chief 
»seph Dam (see photograph 1) during the summer of 1956. Chief Joseph Dam 
a 235-ft-high, concrete-gravity structure located on the Columbia River in 
istern Washington some 50 miles downstream from Grand Coulee. Pressures 
9th upstream and downstream from the tainter gate were measured with 23 
ezometers, and pressure pulsations near the gate seat were measured with 
ee electric pressure cells. At the time the field tests were made, only 
ur power units were in operation; therefore, most of the 300,000-cfs dis- 
large was passed over the 19-bay spillway. For test purposes the reservoir 
901 elevation was altered by changes in discharge during the preceding night. 
hus testing at heads of 32.5, 37.5, 41.6, and 46.5 ft or 0.8, 0.9, 1.0, and 1.1 
mes the design head of 41.6 ft was possible. The piezometer lines were 
srmanently installed in bay 9 during construction, as were inserts and con- 
lit for the pressure cells. The cells themselves were installed just prior to 
sting. Piezometer pressures were observed from gages attached to a 
ilve-bank manifold in the access gallery. Pressure-cell data were record- 
1 on an oscillograph, also located in the gallery. 

A section through bay 9 showing piezometers and pressure cells located 
ong the bay center line is shown in Fig. 1. The crest conforms to the 
andard Corps of Engineers shape for overfall spillways. The rising portion 
this shape is a compound curve with radii 0.2 and 0.5 times the design 
sad. The falling portion of the curve conforms to the equation x*°°?= 
0 Hg9-85y, where Hg is the head for which the crest shape is designed. In 
e case of Chief Joseph the crest shape is designed for a head of 41.6 ft 
hereas the head at maximum pool would be 55.4 ft. The practice of design- 
g the crest shape for a head as low as three-fourths of the maximum ex- 
scted head was advocated after model tests indicated that discharge coef- 
cients would be substantially increased and negative pressures would not be 
ccessive. 

Pressure profiles at the design head of 41.6 are shown in Fig. 2 for gate 
yenings of 5 ft, 10 ft, 20 ft, and free overflow. Pressures on the area down- 
ream from the gate are about the same regardless of gate opening, being 
yout minus one foot of water for the 10-ft opening and plus one foot for the 
)-ft opening and for free overflow. These profiles were obtained with the 
ljacent gates closed. Tests were also made with the adjacent gates open full 
id at openings equal to those in the test bay with resulting changes in 
‘essures being less than a foot or two of water. Furthermore, there was 
tle difference in pressures measured along the center line of the spillway 
y and along a parallel line one foot from the pier. However, for the free 
‘erflow condition positive pressures on the rising portion of the crest were 
msiderably higher along the center line than near the pier. 

Comparisons of pressure profiles from model tests of the dam conducted 
the Bonneville Hydraulic Laboratory and the prototype at 46.5-ft head, or 
1 times the design head, are shown in Fig. 3 for 20- and 30-ft gate openings 
d for free overflow. Negative pressures did not exceed minus two feet of 
iter. It is worthy of note that good correlation between model and prototype 
curred for all comparable tests as would be expected from the Froude law. 

The pressure cells measured the magnitude and frequency of the pulsations 
‘the crest. Oscillograms of the tests at design head are shown in Fig. 4, 

11 A was located flush with the concrete surface five feet upstream from the 
te seat. Cell B was one foot upstream and cell C was 9 feet downstream 
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Photograph 1. Chief Joseph Dam 


Photograph 2. Pine Flat Dam 


FIELD INVESTIGATIONS 1532-5 


SCE 


6 AVE AVM 11IdS HONOYHL NOILOAS 
GNV SNOILVD071 17139 
BYNSS3Yd ‘Y3ILIWOZAld LSAYD 


a1vos 
6 AVG AVMT1IdS HONOYHL NOILD|S 


p- 


f- 4 


WVd Hd3sor 43IHD 


SNOILVDO71 11399 3YNSS3Yd GNV Y3LIWOZIld LS3YD 


SNOILVLS 
09+0! Or+ol OzZ+Ol 00+0! 


SYaLINOZIld 


$1139 3YNSS3ud 


TSW 1334 NI NOILWA313 


f— 


ais 


o0c8 77 
AWITIV9 SS399V 


OGE8 77 
s39v9 is3l 


14917 =PH 


P 
seo 


& 
$10 


6 


HO'Z=¢9X 


a 


CEM 

677 LS749 
eee ae 

ae SaaS EE 


i 0096 77 1 


1o1s 9071-dO01S 


FIGURE | 


February, 195 


HYe 


1532-6 


9EPE 100d 


S3YNSS3Yd AVMT1IdS 
Wvd Hd3SOP 43IHD 


SNOILVLS 


oo+0! 06+6 Ov+0l O€ +01 02+0! O1+01 00+01 06+6 


ie 


MO14483A0 


— 


Ov+0l O€ +0! 02+! O1+01 


' 
EE 


: aan ue 
ONINSdO 3LV9 14-01 
— a fa SS 


~~ 
oS 


TSW 1344 NI NOILWAa13 


FIGURE q 


FIELD INVESTIGATIONS 1532-7 


ICE 


SONIN3dO 3LV9 174 
GNv ‘L4-O0€ ‘14-02 
S3AYNSSAYd AVMT1IdS 


Wvd Hd3SOF J3IHD 


ONIN3dO 3LV9 11N4 
SNOILVLS 
OS+Ol_Ov+O] OE+O! OZ+0! IFO] OO+0! O6+6 


876 =NO|LVA313 100d 
A11VNOA N3dO S3LV9 11¥ :3LON 


ONIN3dO 31V9 14-0€ 
SNOILWLS 
OS+0! OF+0l O€+0! OZ+0l Ol1+01 OOF! 06+6 


adALOLOYd 
AGN ——— 


GN39371 


ONINAdO 3Livd 14-02 
SNOILVLS 
OS+0I OF+Ol OE+OI OZ+0l Ol+0!_00F01 


ASW 1334 NI NOILWA313 


FIGURE 3 


February, 195 


HY 1 


1532-8 


(PH) €%6 100d 
SNOILVWS1Nd 3YNSS3Yd LS3YD 
WvO Hd3SsOfr 43IHD 


13-01 SONIN3SdO 31V9 


SONA INGYNL 3WN14 


14-S ONIN3SdO 3LV9 


FIGURE 4 


SCE FIELD INVESTIGATIONS 1532-9 


‘om the gate. The pressure pulsations were of small magnitude and irregu- 
rr frequency varying from 0.1 to 1.1 ft of water at frequencies of 0.1 to 100 
ycles per sec. It is interesting that the turbulence at cell C, downstream 
‘om the gate, was much greater at small gate openings than that at A or B, 
stream of the gate. This is believed to be associated with the development 
‘the turbulent boundary layer. As the gate opening increased, turbulence 
stream from the gate seat also increased. At full gate opening (free over- 
ow) turbulence was about the same at all locations. The pattern clusely re- 
smbles velocity pulsations observed in flume studies. 

The only condition at which any tendency towards uniform, sinusoidal pul- 
ations occurred was at a 20-ft gate opening with the adjacent gates closed 
hile testing at a head of 32.5 ft or 0.8 design head. At this particular head 
ie jet broke loose from the gate at an opening of 25 ft. Therefore, flow con- 
tions at the crest were extremely turbulent. In fact a violent roller alter- 
ately formed and broke against the gate. It is doubtful that operation would 
Irmally occur under these conditions since it could be eliminated by rather 
mall changes in gate openings. It is to be noted that the maximum magnitude 
the pulsations shown in Fig. 5 was about 1/4 psi at 8-sec periods. The 
sriodic oscillations were damped after three or four cycles. 

As far as the author knows the tests at Chief Joseph spillway are the first 
) be conducted at a head approaching the head for which the crest was de- 
gned. 


Outlet Works Tests 


Outlet works have been tested more extensively than spillways because of 
ore installed facilities and greater availability of water for test purposes. 
ested outlets vary from 4-ft by 6-ft sluices to 25-ft-diam tunnels and tests 
iclude measurements of pressures in intake curves, gate liners, conduits 
roper, exit portals, buckets, and stilling basins. In addition pressures and 
bration have been measured on a variety of control gates including hy- 
raulically operated slide gates, cable-operated tractor gates, and tainter 
ites. Rather extensive field’tests were made on the Fort Peck$ cylinder 
ite and control shaft to investigate pressures, vibration, and air demand. 
ur most extensive installation for testing sluices at high heads is at Pine 
lat Dam (see photograph 2) located on the Kings River near Fresno, Califor- 
a. 

This dam is a 440-ft-high, concrete-gravity structure with an upper tier 
id lower tier of five sluices each. A section through the spillway and sluices 
shown in Fig. 6. Sixty-two piezometers have been installed in one of the 
-ft-wide by 9-ft-high lower sluices and 71 piezometers in an upper sluice. 
ourteen series of pressure measurements have been made in the lower 
uice at heads varying from 110 to 300 ft. Data from the closely spaced 
ezometers in the entrance curve have been reduced to pressure drop coef- 
cients and compared with model data for similarly shaped elliptical intake 
irves based on the equation X2 + y2 = 1, where D is the dimension across 

D2 (2 
3 


Corps of Engineers, Waterways Experiment Station, Hydraulic Prototype 


Tests, Control Shaft 4, Fort Peck Dam, Missouri River, Montana. Techni- 
-cal Memorandum No. 2-402, April 1955. 
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the conduit in the direction concerned. The prototype coefficients shown in ~ 
Fig. 7 are averages obtained from the 14 tests just mentioned. The model — 
coefficients are interpolated from model data on vertical and 10° sloping up- 
stream face to correspond with the 2°52' sloping face at Pine Flat. The 
prototype data indicate slightly greater coefficients than the model. 

The superiority of the 45-degree lip over the old flat bottom for slide gat 
was demonstrated in model and prototype tests of the Norfork Dam4 gates. 
Tests were made at Pine Flat® to measure vertical and flexural vibration of 
the 5- by 9-ft, 45-degree lip, slide gates. Acceleration measurements indicat 
ed stable vibration conditions for all gate openings. These gates have been 
operated for long periods of time at partial gate openings with no apparent ~ 
damage. Piezometric measurements of pressures on the gate leaf have been 
made for heads ranging from 100 to 256 ft. Pressures for gate openings of : 
ft and 5 ft observed at heads of 105 and 256 ft are shown in Fig. 8. 

Air is admitted to the sluices through a vent pipe joining the roof of the 
gate chamber just downstream from the gate. Velocities in the 30-in.-diam 
air vent have been measured for a range of heads varying from 92 to 370 ft. 
A plot of air demand vs gate opening in Fig. 9 shows peak demands at small 
openings and at about 60 to 80 per cent opening. This substantiates similar 
measurements at other projects including Norfork and John H. Kerr. 

One of the earlier installations of prototype testing facilities is in a 20-ft- 
diam, cut-and-cover type conduit at Denison Dam.6 Piezometric pressure 
measurements have been made using the 47 piezometers installed during 
construction. Piezometer locations and the results of piezometric measure- 
ments at a 98-ft head are shown in Fig. 10. It is believed that the flat, uni- 
form slope of the grade line in the downstream portion of the conduit is the 
friction gradient. The Darcy-Weisbach friction coefficient computed from — 
this slope is 0.00644 at a Reynolds number of 120 million. These values whe 
plotted on a chart for friction factors in concrete conduits fall about one- 
fourth of the distance between the smooth- and rough-pipe curves developed 
by Von Karman and Prandtl (see Fig. 11). Pressure-cell tests at the Fort 
Peck 24-ft-8-in.-diam concrete conduit indicated an “f” value of 0.0069 at a 
Reynolds number of 60 million. 

Discharges through tunnels of earth dams are frequently controlled by 
cable-operated tractor gates. Usually 2, 3, or 4 gates separated by piers ar 
used. Pressure and vibration measurements were made on the 45-degree 
lip, 11- by 23-ft intake gates controlling flow into a 22-ft-diam tunnel at For 
Randall Dam,’ on the Missouri River. Vibration tests at heads of 79, 97, ane 
110 ft revealed no unstable oscillations that would indicate resonance with 


4, Corps of Engineers, Waterways Experiment Station, Slide Gate Tests, 
Norfork Dam, North Fork River, Arkansas; Model and Prototype Investi- 
gations. Technical Memorandum No. 2-389, July 1954. j 

5. Corps of Engineers, Waterways Experiment Station, Vibration, Pressure 
and Air-demand Tests in Flood-control Sluice, Pine Flat Dam, Kings . 
River California. Miscellaneous Paper No. 2-75, February 1954. 

6. Corps of Engineers, Waterways Experiment Station, Pressure and Air- 
demand Tests in Flood-control Conduit, Denison Dam, Red River, Okla- 
homa and Texas. Miscellaneous Paper No. 2-31, April 1953. 

7. Corps of Engineers, Waterways Experiment Station, Vibration and 
Pressure-cell Tests, Flood-control Intake Gates, Fort Randall Dam, 


Missouri River, South Dakota. Technical Report No. 2-435, June 1956. 
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with a forcing frequency. Pressures along the gate center line and near the 
side measured with ten pressure cells were positive for heads of 97 and 110° 
ft. The ratio of the average pressure on the gate bottom to the height of pool 
above gate sill is shown in Fig. 12. The average pressures were determined 
for that part of the gate bottom within the conduit proper by integrating con~ 
toured plots of the measured pressures. The pressure-height ratios can be | 
used to estimate pressures on the gate bottom for other pool elevations. 
Water-surface elevations in the gate well were measured with an electrode 
sounding device simultaneously with the pressure-cell tests. The ratio of the 
height of the column of water in the gate well above the tunnel roof to the poo 
tunnel roof depth is shown in Fig. 13. The data in Figs. 12 and 13 are, of 
course, applicable only to the Fort Randall project or to systems geometrica 
ly similar with respect to gate shape and gate-well clearances. 

A discussion of spillways and outlets would not be complete without 
mentioning stilling action. Pressures in spillway and outlet buckets have bee 
measured for nominal flow conditions at a few projects. Pressures were ~ 
measured at two points on a baffle block at McNary. Prototype pressure-ce 
measurements on the elliptical-shaped steps of the Bull Shoals8 stilling basit 
revealed high-frequency pressures whose averages compared favorably with 
those obtained in the model. However, fast pressure drops to as low as -34 
ft of water were indicated by the pressure cells. These low pressures oc- 
curred too fast to be indicated by model piezometric measurements. Con- | 
siderable cavitation damage to the prototype concrete steps occurred. Per= 
haps the lesson to be learned is to be wary of piezometer indications where 
high-frequency, pulsating pressures might be expected such as in gate slots, 
baffle piers, and sudden changes in boundary shapes for high-velocity flows. 

Without going into the details of testing equipment, the types of measurir 
equipment used in the field tests will be discussed briefly. Piezometer 
pressures are usually measured with a pressure gage or a manometer filled 
with mercury or some other fluid. Pulsating pressures are measured with 
electric pressure cells. Vibration is measured with accelerometers. Direc 
writing recorders are used with pressure cells or accelerometers when fre 
quencies are less than 100 cycles per sec. Recorders using a photographic 
process are used for higher frequencies. Air demand has been measured wi 
pitot tubes, velometers, and orifice plates. Water-surface elevations have - 
been observed with the customary stage recorders and electrical-probe de- 
vices. Adaptations of some special testing equipment include deflection gag 
electrical-resistance thermometers, and air-pressure cells. 

The field testing program is a continuing one. Future tests will be made 
at many of the projects already mentioned as additional water for test purpoi 
es becomes available. Facilities have been installed or are now being in- — 
stalled at a number of projects where either water has not yet become avai 
able or construction is not completed. 

In addition to extending the range of measurements that have just been de 
scribed, investigations will be made of crest pressures and pressure pul- | 
sations on a low, ogee crest; velocity profiles and air entrainment will be 
measured on chute-type spillways and at a high overfall dam; pressure pul- 
sations will be measured in the gate slots of a large tunnel; velocity 


8. Corps of Engineers, Waterways Experiment Station, Pressure-cell Tests 


Bull Shoals Dam Stilling Basin. Miscellaneous Paper No. 2-77, Februar’ 
1954, . 
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stribution and turbulence phenomena will be studied with a removable strut 
stalled in a 22-ft-diam, steel-lined tunnel at Fort Randall Dam; and pro- 
sions for measuring sluice discharge by an ultrasonic method are being 
nsidered at a high-head concrete dam. 

In conclusion, the hydraulic field investigations of the past few years have 
rown considerable light on the understanding of many hydraulic problems. 
1e increasing interest of design engineers in providing facilities for making 
sts is encouraging. The development of new equipment and new techniques 
ll facilitate gathering data heretofore considered impossible. Our program 

hydraulic field investigations should continue to be fruitful. 
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SEVEN EXPLORATORY STUDIES IN HYDRAULICS2 


Closure by Hunter Rouse 


HUNTER ROUSE, | M. ASCE.—Because of the evident limitations of class 
ie, none of the seven subjects reported upon in this summary paper was 
npletely covered even in an exploratory manner. Nevertheless, it was felt 
t some knowledge of value to the profession was contained in each study. 

e tone of the discussion of the two studies supervised by Mr. Laursen would 
tify this belief; it is hoped that the five which were not discussed have 

ved equally interesting to the profession. The discussion by Messrs. 
Pherson and Dittig serves to illustrate the complexity of the overfall 

blem. It is readily apparent that the transitions from one regime of flow 
another are especially in need of further investigation. Messrs. McPherson 
| Dittig are correct in assuming that the nappe was not ventilated; in answer 
their direct question, the brink depths which were measured were not con- 
ered precise enough to be significant. The discussion by Mr. Paderi is 
icome because of the information that it provides regarding his own work 
the free overfall. 


Proc. Paper 1038, August, 1956, by Hunter Rouse. 
Director, Iowa Inst. of Hydraulic Research, State Univ. of Iowa, Iowa City, 


Iowa. 
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MEASURING STREAMFLOW UNDER ICE CONDITIONS 


Closure by A. M. Moore 


A. M. MOORE, | A.M. ASCE.—The writer should have covered more fully 
ae of the points brought out by Professor Behlke. Obtaining a good record 
stage is difficult during cold weather even when the temperature does not 
below 0° F. Inlet pipes may freeze or ice may form in the well. Where 
ctricity is readily available engineers have prevented the formation of ice 
leaving a 50- or 75-watt bulb burning in the gage well or by installing soil- 
ting cable in inlet pipes and wells. At many installations a ready source 
ower is not available and the engineer must rely on placing the inlets be- 
the frost line and seating the float that controls registration of stage ina 
e containing 5 or 10 gallons of oil. This pipe must be of such length and so 
ated that there is little danger of the oil escaping at the bottom or top of 
Pipe as the stage fluctuates. The layer of oil must be thick enough so that 
underlying water surface will be below the bottom of the ice in the well. 
an added precaution subfloors are sometimes installed in the wells, below 
frost line. At some stations, ground water several degrees above freezing 
ps into wells and supplies enough warmth to keep well and inlets free from 
. Usually such an occurrence is accidental but in at least one instance 

und water was purposely piped into the well. But as Professor Behlke 

its out, in parts of Alaska and Canada the cold is so intense and prolonged 
- little or nothing can be done to insure a reliable record of stage during 
winter months. The Geological Survey is currently experimenting with a 
ically different type of gaging station requiring no stilling well and with 
recorder driven by a battery-operated motor. This may prove to have 

at merit for cold-weather operation. 

Sold weather can also cause the oil in recorder bearings to congeal, and 
added friction may cause clock stoppage. At least one manufacturer of 
er-stage recorders is now trying out plastic bearings that require no oil, 

$ may minimize one major cause of lost records during cold periods. 

[he writer agrees with Professor Behlke that during prolonged periods 

n the temperature remains many degrees below zero, correlation on a 

y basis between air temperature and stream discharge is impossible. 

The showing of a dashed line representing effective gage height in Fig. 2 

r have led Professor Kolupaila to believe that the author uses the “Stout” 
ffective gage-height method exclusively. The dashed line was included in 
figure merely to show the amount of backwater, and was obtained by 

king backward from the adjusted discharge. The writer does use the 


Proc. Paper 1162, February, 1957, by A. M. Moore. 
Hydr. Engr., U. S. Dept. of the Interior, Geological Survey, Portland, Ore. 
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effective gage-height method, described in the reference quoted by Profess 
Kolupaila, for days when the effect of anchor ice is clearly visible on the re 
corder chart and sometimes for the first or last day of a period of backwz 
from surface ice. Some engineers of the Geological Survey use this me 
throughout the entire period of backwater from surface ice but merely as 0 
of the tools. Actually neither the “Stout” method nor the “Lithuanian” meth 
as such as used by the Survey. Both methods are given consideration along 
with many other factors and perhaps the most reliance is placed in record: 
for adjacent gaging stations and weather records. When adjusted discharge 
and discharge corresponding to gage heights are plotted on the same loga- 
rithmic-scale graph the linear distance between the two graphs actually ist 
cologarithm of “K”, the discharge ratio. Therefore, if during a period of © 
relatively stable ice conditions the adjusted discharge graph is drawn to pa 
through the measurements and to follow the general shape of the unadjustet 
discharge graph, the Lithuanian method is essentially being used without th 
need for actually computing a value of “K”. 

The intent of the Paper was to identify the factors that are considered wW 
winter records are computed by the survey, rather than to compare the rel 
tive merits of various technical details. 9 

The writer does not comprehend the empirical formula for “K” ratio fe 
Elbe River credited to Mr. H. A. Klein. This formula, as shown by Profes 
Kolupaila, is K = 0.308 + .02T, where T = number of days since ice cover 1 
formed. Possibly the formula is intended as a guide for determining adjus 
discharge after the ice has reached its maximum thickness; but it is diffic 
to imagine a hydrometeorological pattern that closely repeats itself each y 
The formula indicates that after ice cover has formed the amount of back- 
water decreases with time and is completely gone 35 days later. 

Professor Kolupaila clearly scored a point in his listing of references. d 
ing with the effect of ice on stream-flow measurement. 

The writer is greatly indebted to Professors Behlke and Kolupaila for # 
cogent comments and discussions that undoubtedly have gone far toward cli 
fying and supplementing the paper. 
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FREQUENCY ANALYSIS OF STREAMFLOW DATA2 


Closure by David K. Todd 


DAVID K. TODD,! J.M. ASCE.—The discussions are a welcome supplement 
his paper. 

Mr. Foster’s statements emphasize points brought out in the paper. His 
tributions and active interest in frequency analysis over many years have 
‘ouraged the proper application of these methods by engineers. 

Mr. Benson’s point that mathematically-fitted frequency curves are not 
erior to graphically-fitted curves is well taken. The basis for his 

tement follows from the fact that it is not possible to establish that a series 
‘lood events must conform to a predetermined statistical distribution. 

hin the range of recorded floods on a stream, the two methods generally 

l yield similar results. Which method, then, is preferable? From an engi- 
Tring standpoint it is the author’s opinion that the mathematically -fitted 
quency curve is superior. By fitting observed data with a curve based upon 
pecified statistical procedure, a consistent result, eliminating the factor 
luman judgment and completely reproducible, is obtained. This is im- 

tant for planning and design work. 

Mr. Benson questions fitting two straight lines to frequency data where 

yds are influenced by two distinct causes. It must be granted that in situ- 
ms where the two types of floods extend over the same discharge range, 

| frequency curves can only be fitted by separating the data and fitting 

ves to each population. More common, however, are situations where the 
/causes produce floods of different ranges which overlap, or merge, to 

y a small extent. Here it is possible to fit two curves to the data, one for 
higher floods and one for the lower floods, but at the same time recogniz- 
that a transition between the two types exists in the vicinity of the inter- 
tion of the curves. Examples have been reported within the last year by 
ter? for a group of 70 watersheds and by Williams? for trade-wind and 
ricane floods in the West Indies. 
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BUTTERFLY VALVE FLOW CHARACTERISTICS2 
Closure by M. B. McPherson, H. S. Strausser and J. C. Williams, Jr. 


M. B. McPHERSON, ! A.M. ASCE, H. S. STRAUSSER,2 A.M. ASCE, AND 
C. WILLIAMS, JR.,3 J. M. ASCE.—The outstanding classic mathematical 
velopment by Mr. Turgut Sarpkaya has more than offset the disappointingly 
aall amount of additional data furnished by manufacturers. 

The profile of the Rockwell 6-inch valve (at centerline) and measurements 
tained with the two-dimensional apparatus appear in Table X. Values of 

calculated by means of these measurements compare quite well with the 
yw test data given by Mr. Henry Voltmann. 

Table XI compares free discharge values of Cg from the 4-inch CDC valve 
yw tests, from Mr. Sarpkaya’s Fig. 5, and as calculated from the two-di- 
ensional apparatus. Mr. Sarpkaya’s theoretical treatment for free discharge 
is restricted to a blade of infinitesimal thickness and presumed that Zy, 
uld not extend beyond the valve body. The two-dimensional results, where- 
the centerline valve profile was used, are quite similar to those determined 

Mr. Sarpkaya. This similarity, despite a lack of correlation of blade pro- 
€, suggests that the empirical adjustment of Eq. (8) by the authors might 
ve corrected more for measurement errors in y than for disparities be- 
een the two and three-dimensional cases. 

When this research was initiated the effects of installation, blade shape 
d closure angle were unknowns and consideration of a theoretical treatment 
y others) was abandoned in favor of more general two-dimensional tests. 
ice the results of the program indicate a relative independence from blade 
ape, the theoretical analysis by Mr. Sarpkaya is now especially appropriate. 
te in Tables VI, VII, and VIII of the paper the comparative constancy in 
ea ratios for all types of installations. Therefore, the results obtained by 
¢, Sarpkaya in considering a thin disk and a free jet should also be applica- 
2 to estimates for other installations. In Table XII are values of Cg from 
ww tests, compared with those determined by means of two-dimensional tests 
lotted in the authors’ Fig. 1) and those calculated from Mr. Sarpkaya’s Fig. 
hrough use of Eqs. (7-a) and (7-b). 

Although time did not permit of further study, it seems quite likely that the 
e of values of o and £, for the actual leading blade tips of the various 
lves, in Mr. Sarpkaya’s equations should yield even closer estimates of 
tual coefficients, although some approximation would be inherent in such a 
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TABLE XI 


FREE DISCHARGE 


h=™N e CDC 
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ocedure, with two sets of angles involved. 

Mr. Sarpkaya has compared his curve for 8 = 90° in his Fig. 6 with the 
snissiat and Dow data and states: “The agreement between the theoretical 
d experimental results is excellent”, Inasmuch as the Génissiat data are 
r a converging valve body, far from satisfying the parallel walls of Mr. 
Tpkaya’s analyses, and the Dow data are not necessarily for B = 90°, the 
reement must be considered as accidental. Otherwise the value of the 
jalyses would be diminished by this rather loose conjecture. 

The comments by Mr. Sarpkaya on the distribution of stagnation points 
1eoretical) are very interesting and deserve consideration by future re- 
archers, The authors had begun to study the correlation of stagnation point 
th pressure distribution and hydrodynamic torque, but interest subsided as 
e torque problem diminished with the introduction of more modern actuators. 

With reference to free discharge versus submerged characteristics, the 
thors wish to reiterate: “----, it would be logical to assume that the flow 
aracteristics might not be identical even though the installations are quite 
milar.” The comment “advocated as a panacea” by Mr. Sarpkaya is hardly 
msistent with the above originally suggested hypothesis. In fact, the con- 
usions reached in his reference 11 were based upon an acknowledged “ex- 
oratory study”. Further, he implies a larger coefficient for the submerged 
se, when in Table II it is larger only for 10° to 30°. There is no argument 
latever on the need for “further study and a better explanation” of submerged 
lid flow characteristics. 

The discussion by Mr. Ronald E. Nece was mostly confined to a consider- 
ion of continuous pipe installations. The introduction of another factor, Ky, 
es not add particularly to the clarity of the data presented. 

The CDC, Rockwell, Netsch and Schulz, and Dickey and Coplen valves had 
n00th walls without “stops, recesses, or projections”. Wall conditions for 
e others were not given. Mr. Nece quite properly points out the possibility 
divergence in results, due to probable dissimilar wall conditions in proto- 
oe valves (particularly in the larger sizes) as opposed to the smaller valves 
rmally tested. 

In the paper it was stated that “a deduction from the overall loss had been 
ade for equivalent pipe friction to and from the valve body for purposes of 
rrelation”. For the Netsch and Schulz data presented in Table IV, one-tenth 
a velocity head was deducted for the flared conduit entrance, and a friction 
ctor of 0.014 for the constant Reynolds number of their tests was used to 
praise the conduit friction to and from the valve body. A sizable disparity 
uld therefore be found in the Netsch and Schulz values of Table IV convert- 

to Mr. Nece’s table. Mr. Nece’s Eq. (4) applies only to a circular valve; 
obtain Ky, for the square damper of Table IV the right side of this equation 
ist be multiplied by G 2, Should the reader check against Mrs. Cohn’s 


aphs, it is to be not i‘ at her plot for Disk II is incorrect inasmuch as she 

i not use the “open area” in interpreting Fig. 14 of Dickey and Coplen, as 

pulated. The head loss for the Dickey and Coplen data was stated to have 

en for the “damper only” and hence a deduction for conduit friction was not 

ide for their data in Table IV. 

In Table XIII are compared the theoretical “average C, from von Mises” 

- free discharge by Mr. Nece and similar values from Mr. Sarpkaya’s Fig. 
Mr. Nece assumed, as an expedient approximation, that the flow was divid- 
equally to each side of the blade. The exact theoretical values by Mr. 
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0° 
10° 
20° 
30° 
1,0° 
45° 
50° 
60° 
70° 
80° 


90° 


(* Value of 0.660 from p.38, "Engineering Hydraulics, by Rouse (22) 
of, =15°, fp =90°, where unbalanced distribution of flow on 
side of 2-dim. case was taken into consideration; used values 


from von Mises analysis). 
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TABLE XIIT 
COMPARISON=COEFFICIENT OF CONTRACTION 


"AVG. Co (Cey + Ceo)/2 

(VON MISES 90° 
2=DIM.)" 
(NECE) ~ _(SARPKAYA) 
1.000 1.000 
0.862 O.7k1 
0.773 0.711 
0.716 0.689 
0.675 0.670 
(0.660%) 0.662 
0.657 0.655 
0.632 0.60 
0.618 02629 
0.612 0.619 
0.611 0.611 
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irpkaya show this approximation to be too gross, particularly for the near- 
yen values. A value of 0.660 at 45° using the von Mises figures, for equal jet 
slocities and hence unequal distribution of flow, is compared with the 0.662 
om the Sarpkaya analysis. The Sarpkaya analysis is now superior to an 
laptation for butterfly valves of the von Mises analysis because the latter is 
sstricted to a 8 of 90°, while the former is solvable for any closure angle. 
oth have a failing as regards accounting for blade shape in determining theo- 
stical contraction coefficients. 

Mr. Nece’s reference D-3 deserves a careful review by anyone interested 
.a comparison of flow characteristics between various types of valves. 

Mr. Henry Johnson deplores the use of the two-dimensional tests because 
€ quantitative results as compared with the prototypes were not perfect, and 
sclares that “prediction of cavitation inception ------ on the basis of such 
sts would be meaningless”. The implications of Table IX have escaped Mr. 
yhnson. The assumptions of Table IX include usage of Eq. (7-b) with proto- 
pe flow data and usage of an average jet velocity rather than a local veloci- 
; the two-dimensional tests were not involved. The two-dimensional tests 
ere presented mostly to enhance confidence in the equations for Co, par- 
cularly Eq. (7-a) and (7-b), and to suggest the possible value of comparative 
yo-dimensional studies in the development or appraisal of a unique geome- 
y. In recently reported gate valve cavitation tests, 14) for positions up to 
1% of full open, audible cavitation was obtained, as would be expected, at a 
wer total head at Station 1 than anticipated “based on Borda loss”. Signifi- 
intly, the cavitation index for audible cavitation was roughly equivalent to 
e anticipated indices plus a constant. Similar descriptions of cavitation 
sts for butterfly valves are noticeably absent from the literature; recourse 
approximation methods is dictated by default. 

Unfortunately, the authors’ Eq. (1) for Cg is useful only for a circular 
ilve. An equitable comparison of the Dickey-Coplen circular Disk I and its 
yunterpart, the square damper, can be made only by multiplying the Ca 
ilues of the square damper by 7 as shown in Table XIV. This comparison 

4 
ids further evidence to the correlation of two-dimensional relationships with 
circular valve, despite the opinion of Mr. Johnson to the contrary. 

A description of the mechanical and operating features of butterfly valves, 
summary of the reports of 1,100 users, has recently been presented. (15 

The intense interest in the subject displayed by all of the discussers and 
e quality of the material added by the discussions has been very gratifying. 


REFERENCES 


.. Ball, J. W., “Cavitation Characteristics of Gate Valves and Globe Valves 
Used as Flow Regulators Under Heads up to About 125 Feet”, Trans. 
ASME, Vol. 79, No. 6, August 1957. 


. Mahon, Ross L., “Hydraulic Butterfly Valves”, Civil Engineering, pp. 52- 
57, October 1957. 


1558-16 HY 1 February, 1 


TABLE XIV 


FLOW COEFFICIENT, Co=CONTINUOUS PIPE 
COMPARISON OF VALVE SHAPE IN TERMS OF AREA RATHER THAN D2 


(DATA FROM TABLE IV) 


DICKEY=COPLEN  SQUARE-DAMPER 


DISK I (THIN) | DICKEY-COPLEN 
Xe (tres —) 
100 1.87 1.85 
20° 1.00 1.01 
30° 0.63 0.60 
0° 0.39 Ohh 
50° 0.22 0.21 
60° 0.120 0.116 
70° 0.062 0.05) 


80° 0.028 0.022 
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MECHANICS OF SEDIMENT-RIPPLE FORMATION2 


Discussion by Vito A. Vanoni and Norman H. Brooks 


VITO A. VANONI,! M. ASCE and NORMAN H. BROOKS,2 J.M. ASCE.—The 
thor has presented an interesting paper on an important subj ect, reviewing 
> various explanations for the occurrence of dunes and ripples on the beds 
streams, and developing an analogy between sand-ripple formation and in- 
bility of an interface between two liquids of differing density. His analogy 
pends on the basic assumptions that (a) a sand bed behaves like a dense 

id, and (b) sand waves and interfacial waves are basically caused by the 

me fluid forces. 

The writers find this analogy difficult to accept. At low velocities when 
ples are just forming, the sand bed does not move or deform like a fluid. 
gples are built up through the rolling and sliding of sediment particles over 
> surface, with the particles being entrained at some places and stopped at 
ier places. During this process only the grains at the very surface are 

e to move; hence, near the critical velocities the growth of ripples may be 
remely slow. 

By contrast, waves at an interface between two fluids may be quickly gener- 
1d because the denser fluid may deform as a continuous mass. A lowering 
the interface at one point results immediately in a corresponding rise else- 
ere because the whole dense fluid layer may flow in response to differential 
rmal stresses. On the other hand a sand bed, being firm, does not respond 
normal stresses at all, and the particles are constrained to tangential motion 
the surface only. (These conditions prevail, of course, only when the sedi- 
nt is just starting to move, and not at high velocities when the sediment 
nsport rate is large and the dunes are obliterated again.) Consequently, 
writers believe a theory for ripple formation must take into account the 
inular properties of the bed, which cannot be adequately represented as a 
ise liquid. 

To derive Eq. (10), it is actually not necessary to use stability theory for 
interface, as dimensional analysis will suffice. Whether or not ripples 

1 form in a flat sand bed depends on a critical relation between the follow- 
variables: T, the bed shear stress; p, the mass density of the fluid; d, 
meter of particles; 1, the fluid viscosity; and w, the settling velocity of the 
ticles. For simplicity, shape and size-distribution factors are not con- 
ered. A straightforward analysis immediately yields the dimensionless 
nbers U,/w and U a/v and Eq. (10) for the criterion of ripple formation. 
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The quantity U,d/v is simply a type of Reynolds number (sometimes called 
the bed Reynolds number) and need not be considered as an “instability inde, 
In fact, it is hardly appropriate to call this an instability index when the cri: 
terion for ripple formation is practically independent of it when U,d/v is 
greater than 100 in Fig. 9 or greater than 10 in Fig. 11. 

The five variables listed above embody all the important components of f 
problem, i.e., applied force (T), fluid properties (p, ), and sediment prope 
ties (w, d). In regard to the latter, it should be noted that w, d, p, and yw 
uniquely determine the submerged unit weight of the particles because it is” 
the only unknown in the drag formula 


| 
Sig 77 Sep ance me 


wherein Cp, the drag coefficient for spheres, is determined from the value ¢ 
wdp/wu (by an experimental curve found in any fluid mechanics textbook). In 
fact, if it is preferred to use (Vg - Y) as the fifth variable instead of w, the 
resultdn place of Eq. (10) is 


baler 
gerd =H? 


The author presents the experimental data in Fig. 9 in the form of Eq. (10) 
whereas Fig. 11 is in the form of Eq. (12). Although there is more scatter a 
points in Fig. 12, it does not necessarily mean that Eq. (12) is not as good as 
Eq. (10), as implied by the author in the last paragraph before the summa 
inasmuch as 


eT a ee 
(y-y)d 3/4 Cy pw" 3 Cy w* 


by Eq. (11) and the definition of the shear velocity U,. Thus the ordinate in 
Fig. 11 is essentially the square of that in Fig. 9, and the scatter of points 
should, accordingly, be about twice as great on a logarithmic graph, as is th 
case. 
The present writers are particularly interested in Figs. 10 and 11 which 
show the Shields’ curve for beginning of motion of sediment and the curve 
the beginning of ripples. These two curves coincide for bed Reynolds numbe 
in excess of 100, but for smaller Reynolds numbers the curves indicate that 
ripples start at higher values of U, than are required to start motion of the 
sediment. The writers have wondered if the two curves should, perhaps, no! 
coincide over the entire range, thus indicating that ripples always start as 
soon as motion of the sediment starts. ; 
In suggesting this, it is realized that there is the danger of being accuseé 
of flying into the face of fact in the form of experimental data. However, iti 
also realized that there is appreciable uncertainty in the results of sediment 
transport experiments, especially when the transport rate is very low, as it 
is in the experiments reported by the author. Under these circumstances it 
takes considerable time to move enough sediment to determine if a ripple 
pattern is developing, and the natural tendency is to err in the direction of | 
taining values of U, for beginning of ripples that are too large. An example 
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ch a difficulty is presented by the author in describing the experiments in 
lich the sand waves of Fig. 6 were obtained. As an example of the slowness 
th which bed ripples develop, the writers cite one of their experiments, (1) 
ade at low velocity in a flume 33.5 inches wide with a fine sand bed which 

is made level at the beginning of the run. When the flow was started, ripples 
arted forming slowly in several spots and finally after about 30 hours of 
ntinuous running the bed reached its steady state configuration of dunes. 
iring this development stage the slope, bed shear stress, and friction factor 
l gradually increased. Pertinent data for this run are given in two columns 

Table 2, one for 5.1 hours after the start of the run (when the first good 
ope determination was made), and the other at the end of the run, 55.3 hours 
ter starting. Although no concentration measurements were made during 
€ period of dune growth, visual observations indicated that the concentration 
creased as the dunes developed; after reaching equilibrium the total trans- 
rt concentration was only 3.3 ppm. 

In experiments of this kind, the shear velocity U, increases considerably 
om its initial value for the same depth and velocity because of the increased 
ughness. The run cited above (Tabie 2) is a good illustration, showing a 
ange in U, from 0.037 fps near the beginning of the run to 0.078 fps at the 
dofthe run. The corresponding changes in the dimensionless parameters 
/w, U,d/v , and T/(7s - ¥)d are also indicated in Table 2; if these data are 
atted on Figs. 9 and 11, it is apparent that the development of dunes in the 
urse of a run moves the data points considerably up and to the right. It 
ems likely that this effect may have contributed to a fictitious divergence 
tween the curves for initiation of sediment movement and ripple formation 
Figs. 10 and 11, especially in the case of data from the other sources. 
Most of the data plotted in Fig. 11 are not available to the writers, so it is 
ficult to judge their precision. It would be of considerable assistance to 
aders if the author would give the specific references from which these 
ta were taken. Without more knowledge of the data used, the writers would 
ll be inclined to move the curve for beginning of ripples, in Fig. 11, down- 
rd because of the natural tendency to obtain too large values of U, in the 
periments for the reasons discussed above. It might not be unreasonable 
move this curve down so it fitted the lowest points, in which case the curve 
¢ beginning of ripples and the Shields’ curve would be quite close. It is ob- 
us that this question cannot be settled with the information at hand, but that 
will be necessary to make more of the tedious and time-consuming experi- 
snts of the kind performed by the author before the matter can be settled. 
The formation of dunes from a flat bed is indeed a puzzling phenomenon. 

e writers do not concur with the author in the opinion that turbulence and 
all irregularities of the bed are not primary causes of ripple formation. 
though the pattern which develops may be quite regular, the natural turbu- 
ice in the stream may still be the instigator. When U,d/v becomes small, 
is believed that the frequency with which eddies penetrate the laminar sub- 
rer decreases but that the sublayer never becomes a perfect “insulator” for 
» bed even for U,d/v < 4 as suggested by the author. When an eddy sweeps 
wn to the bed, particles on the bed are disturbed. Where the particles are 
ed a few grains higher than the initial surface, a small wake develops which 
kes the particles immediately downstream more likely to be entrained. 

e effect of these leeside eddies carries downstream only several times the 
ight of the obstruction, so sand may accumulate to form another tiny ridge. 
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Another wake is formed and so on, with a regular ripple pattern developing 
from a random turbulent motion, as described by the author’s quotation from 
lis. 

ne is easily observed in a flume that the initial ripples have the shortest 
wave length, and that as time passes the wave length increases as the height 
increases, as would be expected according to the above explanation. After a 
sufficiently long time some equilibrium height and wave length are reached, 
the values of which probably depend principally on the grain size, depth, 
settling velocity and flow velocity. When the dunes are fully developed the ~ 
turbulence characteristics are certainly altered, but this fact does not a p 
exclude turbulence as a major cause, On the other hand, the writers concur 
with the author that the occurrence of ripples or dunes does not depend on ~ 
there being a free surface. This fact is conclusively demonstrated by obser 
vations of dunes in pipes by Ismail 2) and Craven. (3) 

In regard to the author’s experiments, several points are of interest. ; 
of all it should be noted that, with the exception of Run 1-4, the depths are ox 
ceeding small, ranging from 0.010 ft (1/8 in.) for Run 1-3 to 0.040 ft (1/2 ir 
for Run 1-6. Even apart from the problem of accurately maintaining and 
measuring such small depths, the validity of the results is open to question 
because of the very low Reynolds numbers. With the aid ofa Stanton-Moody 
pipe-friction diagram, it may be estimated that for Run 1-3 the friction facte 
f = 0,04, V = mean velocity = 0.75 fps and Reynolds number 4Vd/v = 3000; 
Run 1-6 the corresponding estimates are f = 0.05, V = 1.05 fps and 4Vd/v = 
15,000. If turbulence is an important factor in ripple formation, then the lad 
of ‘fully developed turbulence, as indicated by the low Reynolds numbers for 
all but Run 1-4, must exert an important but undetermined effect on the re-~ 
sults. 

The ripple photographs in Figs. 6 and 8 lead the writers to believe that thi 
ripples observed by the author may not have been typical. The shape of the ” 
dunes in Fig. 8 appear more like antidunes(4) than true ripples; the crests 
appear to be quite rounded and the profile rather symmetrical. On the con-~ 
trary, the writers have observed that small ripples are characterized by 
gentle upstream slopes, sharp crests, and steep downstream faces. From 
Fig. 8 the wave length appears to be about 2-1/4 inches or about 18 times th 
mean depth, a ratio which is more plausible for antidunes than ripples. If th 
Froude number is near or above the value one, then antidunes should in fact 
be expected. It is estimated that the Froude number for this illustration (Ru 
1-3) is V/g/gd = 1.3. It would not be surprising if the height of the antidunes 
was of the same order as the mean depth. 

Likewise, in Fig. 6 the oblique crests do not appear typical of ripples oc- 
curring at low velocities. They are perhaps similar to the sand fronts de- — 
scribed by Brooks(5) which occur at an advanced stage of transport when th 
dunes are about to disappear and the bed to become flat. 

The condition of beginning of ripples is an interesting one, but in fact it 
occurs seldom, if ever, under natural conditions. An alluvial river produces 
ripples or dunes on its bed at low stages and may completely obliterate or 
appreciably change them at higher stages. However, when the stage drops 
after a flood, dunes with re-form. If the stage drops to the point where sedi 
ment motion is about to cease, (i.e., to the point where dunes would be expec 
ed to start forming if the bed were initially flat), dunes already exist and 
represent the stable condition for the existing flow. If the flow or stage wer! 


CE DISCUSSION 1558-21 


TABLE 2 


Summary of Data for Run 2-12 (A) 


Hours since beginning of run 


Median sand size, 0.137 


Mean depth, £t 0.541 
Mean velocity, fps 0. 80 
Slope, 0.00039 


0.076 
0.078 
24.6 
0.012 
0.054 


Bed friction factor (adjusted for wall friction) 
J,, Bed shear velocity (adjusted for wall friction) fps 
Temperature, 5s © 


[, Bed shear stress, psf 


Settling velocity, fps 
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reduced still further, sediment motion would cease and the dunes would re- 
main unchanged. This idea is borne out by the presence of dunes in epheme 
streams during their dry period. Hence the writers do not concur with the 
author in the first paragraph on pp. 1197-19, where he indicates that dunes 
will disappear when the flow conditions fall below the curve in Fig. 9. 

Other dune problems, even more interesting than that of initiation, are (a 
how they are modified as the velocity increases, and (b) at what velocity a 
they obliterated. Such changes in the bed configuration of a stream appear f 
be responsible in large part for loops in stage-discharge curves for rivers 
(for example Ref. 6) and for large fluctuations in roughness coefficients, 
which are accompanied by variations in sediment transport rate. Hence the 
dune problems are intimately related to the whole problem of sediment 
movement. Their importance has been recognized widely as evidenced by tl 
attention devoted to it by workers in this field. (7 
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E, ROY TINNEY,! J.M. ASCE.—The physical explanation provided by Mr 
Liu for the mechanics of ripple formation is both interesting and helpful. 
Such an attempt to attack the complex problems of sediment transportation k 
the latest concepts of fluid mechanics is highly commendable. The writer 
feels, however, that the data and analysis presented by the author do not co 
vincingly verify his plausible verbal description of the mechanics of ripple 
formation. : 

In particular, the data presented in Fig. 9 shows that for all values of 
U,d/v greater than about 120 


1. Head, The R. L. Albrook Hydr. Lab., Div. of Industrial Research, State 
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U, = 0.13 (1) 
Ww 


Te = 0.0169 w2 (2) 


is equation is remarkable similar to formulae for critical tractive forces 
initial movement. Such a similarity between the formulae makes it ex- 
mely difficult to present a clear distinction between the two phenomena. 
fore a convincing argument for the mechanics of ripple formation can be 
de, it will be necessary to present data for the beginnings of ripple form- 
on at moderate values of the bed load rather than at the beginning of bed 
d movement. 
Furthermore, the value of U,d/v = 120 agrees extremely well with the 
ue d/é' = 10 (ainge U,d/v = 11.6 d/6') for the lower limit of “rough” 
indaries. The only conclusion that can be drawn from the data on the right 
f of Fig. 9 is that for “rough” boundaries ripples form as soon as motion 
the sediment begins. Just why ripples form on a “rough” bed remains as 
unverified hypothesis. 
One interesting point with regard to Fig. 9 is that the lowest value of U a/v 
- which data are given is 3.0. This corresponds to d/6' = 0.258 which is 
nost identical with the upper limit of “smooth” boundaries. If the author 
ild find no ripples for values of U a/v less than 3.0 it could be concluded 
t ripples do not form on “smooth” boundaries, a point that would strongly 
port the author’s laminar sub-layer instability explanation for incipient 
ple formation. 
With regard to critical tractive force formulae, the writer presents the 
lowing graphs of several critical tractive force functions for small sedi- 
mt sizes. Note that the critical tractive force concept is supported by 
isteins’ theory for bed load particularly for small grain sizes if the critical 
ctive force is defined as the tractive force required for negligibly small 
vement, say G = 1 to 5 lbs. per ft/hour. For sediment sizes over 0.03 ft. 
reral experimenters, namely Shields, Straub, Lane and Carlson, and others, 
i that the equation 


Te = 0.06 (7s - ¥)4 (3) 


; the available data reasonably well. 


M. L. ALBERTSON, ! M. ASCE, D. B. SIMONS, 2 A.M. ASCE, and E. V. 
‘HARDSON, ? J.M. ASCE.—Mr. Liu is to be commended for his excellent 
nmary of the salient literature on the subject of stability of an interface. 
rthermore, he has related very logically this literature and applied it to 
generation of ripples in alluvial streams. He has clearly demonstrated 
t there is a definite relationship between the formation of waves at the 
srface of two liquids, at the interface of sand and air, and at the interface 


e: Published with the approval of the Director of the U. S. Geological 
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sand and water. By this means considerable information already available 

S been brought to bear on the subject of ripple formation. In the opinion of 

> writers he has made a major contribution to this subject. . 

Mr. Liu has climaxed his research with the formation of a single plot which 

Ows Clearly the conditions under which movement of sediment is initiated 

d ripples are formed. His criterion for sediment ripple formation is a real 

Jestone in the struggle to understand the mechanics of sediment transport. 

r. Liu quite logically uses a ratio of the shear velocity to the fall velocity of 

> sediment as his dependent variable and the shear velocity Reynolds 

mber as his independent variable. The scatter of the data in this plot is 

ite small compared with the usual scatter of data involving sediment nee 

rt. Furthermore, data which have been taken more recently by Plate(1 

d a theoretical development which he has made under Mr. Liu’s direction 

ove even more conclusively the general validity of Fig. 9. 

Despite the excellence of Fig. 9 and the discussion preceeding its de- 

lopment, the writers believe further interpretation of this figure would be 

Ipful. The shear velocity Reynolds number for some purposes (and the 

iters believe this is one of them) can be represented more significantly as 

ratio of the sediment diameter d to the thickness of the laminar sublayer 

by the equation for the thickness of the laminar sublayer 6' = 11.6 v_ which 
U 


* 
n be inverted and both sides multiplied by d to give 11.6d _U,d which 


1 

sans that the magnitude of U,d/v is nearly 12 times as peest as the ratio 
6'. If the assumption is made that the sediment size d is approximately 
ual to the size of the sand grain roughness k used by Nikuradse in pipes, the 
tiation of ripple formation might logically be related in some way to the 
ansition from a smooth boundary to a rough boundary in rigid-boundary fluid 
>chanics. 

Fig. 9 is particularly significant in that it shows a transition from one re- 
ne of flow to a second regime of flow. The significance can be understood 
re completely from the following observations: 


1. Ripple formation is definitely dependent upon the ration U,/w. As the 
particle size increases, the shear velocity must increase in direct pro- 
portion if ripples are to be formed. Since the shear velocity is pro- 
portional to the square root of the shear, the shear must be increased 
4 times if the fall velocity of the bed material is doubled. 

2. For small U,d/v -values a much greater magnitude of the ratio U,/w 
is required for ripples to form than for the larger U,d/w-values. The 
entire range of U,/w-values is particularly significant when associated 
with the ratio d/é'. In other words, the formation of ripples is very 
closely related to the size of the sediment d relative to the thickness of 
the laminar sublayer. 

This concept is logical and quite in keeping with Dr. Liu’s hypothesis 
that the laminar sublayer is intimately involved in the phenomenon. 

3. At the minimum U,d/V -value for which there are date, U,d/v = 3 and the 
shear velocity U ,/ WwW required to form ripples is approximately 5 times 
that required at U,d/v = 50. The reason for this is understood more 
clearly when the ratio d/§' is used instead of U,d/v. At U,d/v = 3, 
d/d' is approximately 0.26, which is significantly approaching the point 
in rigid-boundary fluid mechanics where the transition begins from a 
smooth boundary to a rough boundary for Nikuradse’s data using sand 
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grains of uniform size. The data presented by Mr. Liu, however, we 
taken from experiments involving both sediment having a relatively 
narrow size-range and sediment having a wider size-range. q 
4. At U,d/v = 80, where U,/w becomes independent of U,d/v, d/6' =7 
which is approximately the value of k/§' where a rigid boundary be- 
comes completely rough based upon the data of Nikuradse and the cur 
of Colebrook-White. Apparently, the protective effect of the laminar 
sublayer is lost and the magnitude of U,/w necessary to form ripples: 
decreases from 0.6 to 0.12 as the sublayer is destroyed in the transit 
from smooth to rough boundaries. 
5. Once the laminar sublayer is destroyed (U,a/v > 80), the formation 6 
ripples no longer takes place. Evidently, then, Dr. Liu has demonst 
ed that ripples cannot form without a laminar sublayer existing and ~ 
Tison\6) has shown that ripples do not form if the flow is laminar. 
very significant conclusion can now be drawn that ripples form only ij 
the flow is turbulent and the boundary is smooth or in transition from 
smooth to rough—that is, only when a laminar sublayer is in existence 
This supports Dr. Liu’s contention that ripple formation is somehow | 
lated to instability at the laminar sublayer. 7 
6. In view of the foregoing analyses, the writers are inclined to believe 
that the formation of ripples is somehow associated not only with in- 
stability but also with the turbulence in the stream above the laminar 
sublayer and sand bed. Such turbulence in turn, is directly related tc 
the instability of the laminar sublayer. In other words, turbulence ani 
instability of the laminar sublayer are intimately related and cannot 0 
separated. In rigid-boundary fluid mechanics the boundary transition 
(that is, the transition from a smooth boundary to a rough boundary) i 
initiated when periodic fluctuations in velocity caused by overhead 
turbulence (or eddies caused by instabilities of the interface) penetra 
the laminar sublayer to the point where they encounter roughness pro 
jections. Under these conditions a part of the resistance to flow is be 
ing contributed by the form drag caused by the roughness projections. 
As the above turbulence becomes effective, the laminar sublayer is 
made thinner, and the boundary becomes rougher. Hence the local ve 
locity caused by the turbulence encounters the roughness projections 
greater percentage of the time until finally, when the boundary is com 
pletely rough, there is no laminar sublayer acting as a protective 
blanket and (in an alluvial bed) ripples can no longer be formed. 
7. When the boundary is completely rough, the drag is due entirely to 01 
drag around the roughness projections. The drag coefficient then be- 
comes independent of the viscous effects. With a smooth boundary (in 
contrast to a rough boundary) the laminar sublayer acts as a protecti 
blanket which completely covers the roughness projections and the ré 
sistance to flow is a function of the relative viscous effects as reflect 
by U,d/v or the ratio k/6' or d/é'. 


The fact that ripples form only in the smooth to rough boundary transitio 
region leads one to wonder what the difference might be in the nature of the 
ripples which form with a completely smooth boundary (if at all), with a 
boundary in transition, and with a boundary that is nearly rough. 

The fact that Figs. 9 and 10 show clearly the conditions under which sedi 
ment movement begins and ripples are initiated led the writers to speculate 
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the possibility of similar relationships for other types of sand waves such 
dunes, bars, antidunes, and the plane bed. Therefore, Fig. A was prepared 
ng data from several sources where adequate information was available. 
this figure it is clearly evident that, in general, dividing lines can be drawn 
separate the regimes of ripple formation, dune formation, transition, and 
idune formation. Furthermore, Fig. 9 is extended more than one cycle to 
left by means of the data of Kalinske and Hsia(2) and it is extended more 
n one cycle to the right by the data of Lane and Carlson.(3) For purposes 
comparison, the transition function from a smooth to a rough boundary is 

0 included in Fig. A. It should be noted that although the abscissa is the 
ne for both plots, the ordinate scales are entirely different and unrelated. 
netheless, plotting the boundary transition function permits direct compari- 
1 of the data and curves representing alluvial channels with the beginning 
the Colebrook- White transition, the beginning of the Nikuradse transition, 

1 the ends of each transition. In this figure a temperature of 20° C was used 
determine the fall velocity and kinematic viscosity where temperature was 
_available. The scatter among the points for Laursen, 4) Brooks, (5) 
rton,(6) and the USGS-CSU(7) data is the result of temperature variation. 
1es of constant temperature can be drawn for these data at an angle of 45 
srees. It is also noteworthy that for almost any transition line, a change in 
nperature for a given size of bed material could result in a change of re- 
ne. 

The following are observations regarding Fig. A which may have signifi- 
Ice. 


1. Kalinske and Hsia(2) reported ripple formation although they did not 
record the conditions under which movement began and the conditions 
under which ripples were initiated. The minimum value of d/6' for 
their data is approximately at the same point that the Colebrook- White 
transition begins—that is, where the Colebrook-White curve connects 
with the smooth boundary curve. To the left of this point the boundary 
is completely smooth and no influence of the overhead turbulence is felt 
at the boundary. In other words, the laminar sublayer is a completely 
protective blanket. This fact leads one to wonder if it is possible for 
ripples to form under these conditions. In any event, the writers have 
been unable to locate any data for smaller values of d/é'. It may be 
that the rather explosive nature of the fine materials makes it difficult 
to determine definitely when movement begins and ripples begin to de- 
velop. 

2. The curves describing the regimes of bed roughness are rather clearly 
defined and are nearly parallel for d/é' < 10. In Fig. A it is apparent 
that Lui’s curve is not parallel to the foregoing family of curves in this 
range, whereas the Shield’s curve turns up at the left end more nearly 
in accordance with them. Based on the rather definite and systematic 
orientation of the regime curves and Shield’s curve, it is possible that 
Liu’s curve should be lowered in the vicinity of d/6' =1. This results 
in a curve, see Fig. A, which more nearly conforms with the curves de- 
fining the other regimes of flow. Such a minor adjustment might also 
be justified by a careful study of the effect of size distribution of the 
sediment. . ; 

3. For increasing values of d/6' (the laminar sublayer becoming relatively 
thinner and thinner) the influence of viscosity becomes less and less for 
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each type of bed formation. Ripples do not form beyond d/6' ~102 
dunes cease to exist beyond d/6'~ 40. Apparently, the data are suf- 
ficiently complete to permit these approximate conclusions regarding 
both ripples and dunes. For the transition regime of bed formation, — 
however, the data are not sufficiently complete for any conclusion tol DE 
drawn about whether the transition ceases to exist beyond some large © 
value of d/5'. Gilbert(8) made observations which indicate that anti-— 
dunes can form in large size bed material, although as far as the w t- 
ers know only qualitative observations are available to support this 

conclusion. 
. Because of the extreme surface waves and irregularities for both the © 
antidune regime and certain types of transition flow, the Froude numk 
no doubt plays a very important part. No attempt has been made as yé 
to relate the Froude number beyond that done in the table in Fig. A. 
. Although the two variables U,/w and U,d are shown by Fig. A to be of 


primary importance, there is some question about the significant 
variable which will correlate the data in a quantitative manner as now 
done in a qualitative way by the terms: ripples, dunes, transition, and 
antidunes. The Froude number has some significance, as shown in the 
table, in defining the regimes of ripples and dunes, but there is con- 
siderable question whether it will prove to be as significant in these Fr 
gimes as in the antidune regime. The drag coefficient Cp is essentia 
a constant for each size of bed material and hence is not a logical thir 
variable except perhaps in establishing the points where the regimes, 
such as ripples and dunes, cease to exist. The remaining variables — 
which might be significant are the slope S, the size of bed material d 
relative to the depth of flow D, the concentration of total load Ct, the © 
Chezy discharge coefficient C/Ve , and the parameter used by Shields 
T./d Ads- Which of these would be most significant depends, at least’ 
some extent, upon the use of the plot. q 
. Fig. A shows that as d/6' increases the existence of two of the regime 

of flow are eliminated. When d/6' > 7 and sediment is moving, the © 

possibility of the first regime (movement without bed waves) is elimi- 
nated. When d/é' >10, ripples will no longer form and when d/6' > 
40, dunes cannot exist. For the combination of water and sediment ha 
ing specific gravities of 1.0 and 2.65 respectively, Fig. A shows that 
these changes in regime are related approximately to the size of 
sediment—the limiting size for the existence of sediment movement — 
without bed waves, ripples, and dunes being approximately 2 mm, 2.57 
mm, and 8 mm respectively. In the transition regime beyond d/6' =‘ 
several bed configurations are possible—such as a plane bed and a be 
with symmetrical undulations of various sizes. The limiting conditiol 
for the transition regime and antidunes are not established. In fact, 
they may not be in the realm of practical possibility. 

Although the causes of these limits are not certain, they are no dou 
related at least in part to the relative influence of viscosity as ex- 
pressed by d/6'. However, the limits may also be related to the grav 
tational forces relative to the forces creating the sand waves, the dra 
coefficient, the size of the sand waves relative to the size of bed ma- 
terial, and the Froude number. Perhaps both ripples and dunes dis- 
appear when the size of the grain roughness approaches the same ord 
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of magnitude as the size of the dunes and ripples. 

7. The form of the bed in the transition regime is very unstable and sub- 
ject to frequent changes with respect to both time and distance. At 
times two or more types of flow may exist simultaneously side by side 
within a given reach. The following types of bed roughness and water 
surface conditions have been observed by experimenters in the tran- 
sition regime: 


a. A plane bed with no visable sand waves of any kind and a very smooth 
water surface. Based upon flume studies, this condition is unstable 
and is difficult to reproduce, particularly as the length of flume is 
increased beyond 75 ft. 

b. An extremely rough pattern of sand waves of the form typical of the 
dune regime. The water surface above this type of bed is very 
turbulent and strong boils are visible at the surface. 

c. Symmetrical standing waves of low amplitude which form and gradu- 
ally disappear. Unlike antidunes these surface and bed waves have 
no tendency to migrate upstream or break and shear off. The bed 
may have symmetrical undulations; a diagonal dune pattern cross- 
laced like shoe strongs; a washed-out dune pattern; sand bars or any 
combination of these. 


The data indicate that the transition regime occurs when 0.6 < Fy, < 1.2. 
This variation of Fy may account for the instability of this regime. 

That is, a small local increase or decrease in depth changes Fy and 
may change the local conditions from those favorable for one bed con- 
figuration to another. In rigid-boundary flow, a marked instability with 
respect to depth and wave formations also exists at F, > 1.0. The in- 
stability of this regime, especially with the occurrence of sand bars and 
two distinctly different types of flow side by side, may be aggravated 
also by width depth ratios which are too large with respect to slope— 
that is, three dimensional flow. 

8. As d/6' increases beyond the transition regime, the stationary bed 
undulations (which are rather similar to sine waves) increase somewhat 
in size and antidunes, which travel upstream, are formed. For increas- 
ing values of d/6', the formation of antidunes is related to the size of 
the bed material, and as both d/é' and U,/w are increased the formation 
of antidunes is a function of the Froude number. As pointed out by 
Gilbert, the antidunes occur in coarse material of a size even greater 
than that for which data have already been taken. 


Fig. A has considerable usefulness for design. Typical examples of design 
es are tabulated as follows: 


1. As demonstrated by the data of Lane and Carlson, the horizontal line at 
U,/w = 0.12 is the limiting condition for which no movement of sediment 
will occur in a channel. In other words, this line can be used to com- 
pute the maximum flow which can be introduced in a given channel be- 
fore movement of bed material begins. Such a fact can be very useful 
in the design of stable channels and in the design of channel modifi - 
cations, such as cutoffs, for natural channels. If desired, a design 
could be made for somewhat larger values of U,/w without moving ap- 
preciable quantities of sediment. 
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2. 


Although Fig. A offers many possibilities for shedding light on the funda- 
mental mechanics of flow in alluvial channels and on the design of alluvial — 
channels, there are a number of aspects of these problems for which con- 
clusions cannot be drawn without further information. Therefore, further ré 
search and study apparently are needed as follows: 


ie 


2. 


1. Plate, E. J., “Studies of Sediment Ripple Formation,” Thesis (1957) Colo: 
rado State University, Fort Collins, Colorado. 


2. Kalinske, A. A. and Hsia, C. H., “Study of Transport of Fine Sediments by 
Flowing Water,” Bulletin 29, 1945, University of Iowa, Iowa City. 


. Studies should be made for various coarse materials to establish the 


By using the concentration of total load of sediment as the third vari 
able, it is possible (within reasonable limits of accuracy) to design a 
channel in a given size of bed material to carry a desired load of sedi 
ment. 
A channel can be designed to have the desired regime of sand waves, 
Where channel stabilization is the problem, and coarse non- cohesive 
materials are available, the size of the coarser material needed to col 
struct a stable armorcoat for the bed and/or banks can be estimated. 
The size of material required to stabilize the banks is a special 
problem, however, because of the influence of gravity and angle of re- 
pose on the stability of material placed on an inclined plane. 


Additional field data should be plotted on Fig. A to prove whether the 
conclusions reached using the flume data apply to field conditions. 
The influence of size distribution of bed material has not yet been in- 
cluded in Fig. A. It is possible this will have considerable influence i 
some types of flow—especially for small d/é'-values and for initiatio 
of movement. 

In an open channel the shear is not distributed uniformly across the 
channel. Hence the shear velocity U, as used in Fig. A is simply an 
average. For certain flow conditions this non-uniform distribution ma 
have appreciable influence upon the particular flow regime which will 
take place. Furthermore, certain shifting and sorting of the bed and 
bank material may occur. These factors need special study to obtain 
additional refinement. 

There is no doubt some relation between sediment transport by water 
and sediment transport by wind. Data from wind studies also should b 
applied to Fig. A. j 
There is considerable evidence that the height of a sand wave is limite 
not only by the depth of flow but also by other factors such as Cp. 
These should be investigated. 

Data should be taken for very fine material to extend those data of 
Kalinske and Hsia to the extreme limits—the initiation of movement an 
the formation of antidunes. 


transition and antidune regimes where d/6' > 40. 


The magnitude of d/é' at which ripples will no longer form and dunes 
will no longer form should be firmly established. 
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HE EFFICACY OF FLOOR SILLS UNDER DROWNED HYDRAULIC JUMP2 


Discussion by Mushtaq Ahmad 


MUSHTAQ AHMAD.!—The author has presented results of experiments on 
w head barrages working under drowned conditions of flow. He has tried to 
termine the efficiency of the sills against bed erosion by using the vertical 
locity distribution as a criterion. The writer had the opportunity of experi- 
enting for the last 15 years on a number of models of low head weirs, bar- 
ges, and canal falls. In the alluvial plains of Indus, in West Pakistan, hy- 
aulic structures are built on fine erodible sand and scour depths varying 
3m almost zero to as much as 60 feet have been experienced at the works. 
le writer has always felt and advocated(1) that apart from the length charac- 
ristics of hydraulic jump and visual observations of flow conditions, the ve- 
city distribution and comparative studies of scour should be taken as criteria 
r the proper design of stilling basins and scour prevention devices, especi- 
ly for those built on highly erodible material. 

If the vertical velocity distribution in and below the normal and drowned 
mp without and with scour control devices are examined (Fig. 1, a, b, and c), 
will be noted that: 


1. In anormal jump, the shooting flow continues as a live stream under 
the roller of the jump and then expands in the vertical direction. With 
velocity distribution, such as occurs in a natural stream, the velocities 
are maximum at or near the surface decreasing towards the bottom, but 
immediately below the normal jump, the distribution is just the reverse; 
that is, the maximum velocity filaments are near the bed and the veloci- 
ty decreases as we move upward. 

2. For maximum scour control, the designer must attempt to attain the 
normal velocity distribution within the confines of the stilling basin by 
moving the line of maximum velocity upward by means of sills or 
staggered baffle piers. 

3. For a drowned weir or barrage of the type tested by the author, the real 
hydraulic jump does not form at all; the fast moving jet continues to 
expand as it goes into the stilling basin, thereby determining the veloci- 
ty distribution pattern; The efficacy of sills or staggered blocks in the 
stilling basin in correcting velocity distribution or controlling scour 
has been found to be very low when the high velocity jet persists on the 
surface (see Fig. 1c) under conditions of a high drowning ratio. 2) No 
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(b). WITH STAGGERED BLOCKS - CLEAR JUMP 


Line of max. velocity --~. 
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hard and fast rules can be laid down about the best position of the sill 
for such conditions of flow, as variation in discharge, drop, downstream 
depth, etc., can change the length of the surface and bed rollers and 
thereby change the velocity distribution pattern. Such highly drowned 
structures shall always have to be tested on a model. 


The author has obtained a relation for the bed shear but has not actually 

ed it anywhere in the analysis of his data. Vmax, Vb, and L¢/Y, have been 
Vo Vo 

osen by the author as dimensionless variables for the study of scour control 

vices below low head hydraulic structures, but there are many others that 

ter the problem. A functional relationship of the type given below, with per- 

ps even more variables added, needs a thorough study. 


Past -f ( Vi, D(actual), Lg, dmax, Vs, V2D2 , Vmax, Vb, “1 ) 
92/3 gD; Do Dg Dg V2 Y V2 Vo H 


ere V, is the fall velocity of the sand particles and other symbols are 

own in Fig. 2. 

Unfortunately, the author has been working in a region of flow where there 
no hydraulic jump formation in the real sense and where there are many 
ficulties in generalizing results. The writer working on low head weirs and 
rrages with a real hydraulic jump has made comparative studies of the 

our depth in terms of Dg/92/3 (where Dg is scour depth below stilling basin 
or) against L¢/D9 (the optimum length of stilling basin), D(actual)/Dg (the 
timum depth of the stilling basin, h /H (the drowning ratio) and Vp/V2, with 
d without sills and baffle piers or staggered blocks. The conclusions are 
ren below: 


1. When a hydraulic jump forms, two rows of staggered rectangular or 
cubical blocks of height between 1/10 and 1/7 Dg placed at the point 
where the live stream under the roller starts expanding or at about 
2 . 5D» below the upper end of the jump gives best results to move the 
line of maximum velocity upwards without causing hurdling or extra 
disturbances.(3) Another row of staggered blocks of height between 1/5 
to 1/7 Dg near the downstream end of the stilling basin helps to create 
a bed roller and pile up bed material against the toe wall. 

2. A length of stilling basin greater than 6 Dg does not help in proportional 
reduction in scour. 

3. Scour control devices help to move the line of maximum velocity up- 
ward to obtain normal velocity distribution within the confines of the 
stilling basin of length 4 to 5 Dg. , 

4. Two rows of staggered blocks as recommended in (1) above are better 
scour control devices than continuous sills or streamlined blocks. 

5. As the weir becomes more and more drowned, the efficiency of 
staggered blocks as a device for changing velocity distribution and con- 
trolling scour decreases. The staggered blocks have to be made higher 
and have to be shifted downstream from their normal position in order 
to make them effective. 


chanics of Energy Dissipation in the Stilling Basin and Function of Scour 
ntrol Devices 


The author has discussed the changes in the size of the surface and bed 


1558-36 HY 1 February, 195 


vortices with the variation in the position of the sill and also the efficiency ¢ 
different types of sills as corrective measures for velocity distribution and 
scour control. It is admitted by the writer that the hydraulic jump is the bes 
energy dissipation device. The visualization of the mechanism or physical 
process of energy dissipation in a jump can explain some of the above and ~ 
more results obtained by the crates (hk It is well known that the energy per 


pound of flowing water at any point is measured by U_ where U is the mean: 

2g oe 
velocity. Apart from this, there is additional energy due to the fluctuating 
velocity vector and the total kinetic energy is equal to their sum. In a hy- 
draulic jump, at the interface of the high velocity jet and the overlying mas 
of rolling water, lies a region of high shear stress where the kinetic energy 
is transformed into heat, partly caused by the shear and partly by the creatit 
of vortices, eddies, and rollers. Each of the eddies formed and shed off fro 
the region of high stress carries a certain amount of rotational energy in ad 
dition to its energy of translation. These eddies diffuse into the body of the 
liquid in the direction of flow. These small rotating masses of fluid in cont 
with each other cause instantaneous changes in velocity and create surfaces 
of high local shear; thus, the energy of turbulence is slowly dissipated thro 2 
viscous action into heat. From the mathematical treatment of energy dissi- 
pation for isotropic turbulence, G. I. Taylor(5) showed that the rate of energ 
dissipation per unit volume is: 


__=15n U2 
Ww K2 


where U2 is the mean square root of fluctuating velocity component and A ig 
the scale of turbulence or the average size of the eddies. The equation thou 
not quantitatively applicable to nonisotropic turbulence, gives an important 
clue on the factors affecting rate of dissipation of kinetic energy. It shows 
that the rate of turbulent energy dissipation will increase as, ), the size of 
the eddies, decreases for a given intensity of turbulence. This means that 
get maximum energy dissipation, the high velocity jet should be broken inte 
as small vortices as possible and as early as possible in the stilling basin. 

This is what the hydraulic jump does and staggered blocks or sills, etc. 
are supposed to do. Anyone who has seen the process of scour formation in 
fine material below a model in glass flume must have noticed how sediment 
picked at irregular intervals and with varying force as vortices of different 
strength pass over the material. To control scour, therefore, the decay in 
vorticity should take place within the confines of the stilling basin to such 4 
extent that vortices moving along the bed are not strong enough to lift the Sé 
particle in excess of those tending to settle down. To achieve this, it is nee 
sary to push the line of maximum velocity near the surface and to decrease 
bed velocity. 

To summarize, the energy dissipation and scour control measures in a 
stilling basin would comprise: 


a. Form hydraulic jump 
b. Break the high velocity jet into a large number of as small vorticé 
as possible and as early as possible in the stilling basin to ensur E 
maximum energy dissipation 
c. Use scour control devices in the stilling basin to help obtain the 
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above objective by moving the line of maximum velocity towards the 
top and obtaining normal velocity distribution within the stilling 
basin. The bed velocity and strength of vortices moving near the 
erodible bed then decreases 


> above physical concepts of the phenomenon also incidently explain that: 


1. Staggered blocks are better scour control devices than a continuous sill 
which gives a continuous roller. 

2. Any scour control device in a supercritical jet which causes hurdling 
but no breaking of the jet into a large number of small individual eddies, 
is not effective as a scour control device. 

3. Streamlined baffle piers or those with the upstream face sloping or 
stepped are less effective than ones with vertical upstream faces. 

4. If the height of the blocks or sills is increased beyond a particular limit, 
in relation to the depth of flow in the stilling basin, the scour depth in- 
creases again because of additional disturbances and vortices. 
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IS THE WRITING OF FLOOD INSURANCE FEASIBLE?2 


Closure by John F. Neville 


JOHN F. NEVILLE.!—Messrs. Laursen and Toch question the writer’s 
iclusion that “insurance against the peril of flood applicable to fixed proper- 
cannot successfully be written” by suggesting means by which the amounts 
9e paid under a so-called insurance program could be reduced or limited. 
ne of these suggestions alters or contradicts the writer’s basic conclusion 
jarding the insurability of the peril of flood. In considering the insurability 
the peril of flood, certain fundamentals of insurance should not be over- 
ked, particularly the requirements that any insurance plan must be self- 
porting, that the protection offered must be sufficiently attractive to be 
eable, and that enough property owners would in fact buy the type of pro- 
tion offered to give the insurance company the spread of risk needed to 
‘mit the orderly working of the law of averages, and that the plan would not 
ourage selection against the insurance carrier. Another practical aspect 
nsurance is that over a long period premiums must necessarily be great- 
than losses due to the cost of acquiring and servicing business and of main- 
ling the insurance mechanism, notwithstanding the comment of Messrs. 
irsen and Toch that “there is no advantage in having insurance against 

yds of frequent occurrence since the premiums must necessarily be great- 
than the ‘losses’ due to those floods.” 

The discussion of damage by frequent floods as a factor in annual costs of 
intenance of the property presented by Messrs. Laursen and Toch has con- 
erable merit. Where the structure involved is located on the flood plain 
ject to frequent inundation, it may be feasible to fill in the site or place 
uable equipment at a higher elevation so as to eliminate damages from 

se frequent floods. This has actually been done in certain cases. The 

ual cost of such changes in location or design of structure may prove to 
less than the theoretical average annual flood loss. It would obviously be 
he advantage of the owner to incur such additional construction expense 
her than to buy flood insurance with its necessary expense loading, even if 
h insurance were commercially available. Prevention of flood damage ap- 
rs to be a sounder long-range approach to this problem than to seek to 
vide indemnity for flood losses by means financially unsound even though 
oretically quite attractive. 
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ATTENUATION OF SOLITARY WAVES ON A SMOOTH BED® 


Discussion by Yoshiaki Iwasa 


YOSHIAKI IWASA.—As a result of errors in the author’s calculation for 
rabolic flow, several equations in this paper should be corrected. 

As a result of the miscalculation in the shear distribution acting along the 
undary for parabolic flow in Part II, the following correction may be ex- 
essed. 


Location Error Correction 
ght hand of Eq. (18b) = asa 2 foe 
ght hand of Eq. (19b) Set + 60 
ght hand of Eq. (20) + fed ae + iia 
ght hand of Eq. (21) GE rene 
pression in p. 10 On = ea C= 2 
pression in Be rae 


Therefore, the rate of energy dissipation dE/dt for parabolic flow, calculat- 
by the laminar boundary layer thickness, should be divided by Vv 2. 
Finally, the relation (26) of attenuation becomes 


- ~z 
Ao afi 1- 0,8062A, -1,6348A2 -5.0970A~~-| AG kt 


0.8062Ai - 1.6348. A, - 50970A; Poe] = 0.10602R, Cinetee ig, } 


us, the ordinates of the theoretical curve in Fig. 2 should be multiplied by 


ind Fig. 3 may be corrected as follows. 

Evidently the corrected curves in Fig. 3 indicate the behaviours of attenu- 
on for low solitary waves with more accuracy. However, for high waves, 

: present first approximation for the influence of side walls will become in- 


mplete. 


Proc. Paper 1262, June, 1957, by Yoshiaki Iwasa. 
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STUDY OF BUCKET-TYPE ENERGY DISSIPATOR CHARACTERISTICS@ 


Discussion by Gale B. Dougherty 


GALE B. DOUGHERTY,! A.M. ASCE.—The authors and the sponsoring 
m are due the gratitude of the profession for making these data available 
* design use. Some recent studies in connection with the design of Weiss 
m for the Alabama Power Company provide data which confirm some of the 
sults obtained by the authors. Several different buckets were tested in con- 
ection with a 1:48 scale standard crest spillway model. These tests were 
ide in the Fluid Mechanics Laboratory of the University of Alabama and 
onsored by the Alabama Power Company. The model consisted of a center 
y flanked by two half bays. The bay width was 40 feet and the piers were 8 
st thick with ogive noses. The crest was 45 feet above the bucket invert. 
l tests were made with an upstream energy grade line 77 feet above the 
cket invert. Discharges were controlled by radial gates except for the larg- 
t which was free discharge. The erosion bed consisted of 1/4 to 3/8 inch 
a gravel and for most tests was placed 5 feet above the bucket invert. The 
” exit angle buckets tested had radii of 15, 25 and 40 feet. One 40-foot radi- 
bucket had a tangent lip added to place the end sill 10 feet above the bucket 
fert. 

Fig. 1 shows a comparison of the data obtained for the Weiss Dam to the 
thors’ curves (Fig. 3 of the original paper) for height of bucket roller. The 
fa were taken from plotted water surface profiles after the erosion had sta- 
ized. The correlation, while not conclusive, is indicative. No check was 
‘ained for the surge height relationship. As the discharge increases, the 
tted points fall below and to the right of the curves. This is attributed to 
: effect of the flow over the crest. It was noted that the surface of the 
ket roller was not horizontal, but was depressed on the upstream end. This 
sulted in decreasing the value hp. 

The authors indicate that a 45° exit angle was used because others had 
ind it to be the most effective under a variety of service conditions. At the 
iversity of Alabama, a 25-foot radius bucket with a 37-1 /2° angle was test- 
with lip elevations, 5, 7.3 and 10 feet above the bucket invert. The results 
these tests ‘indicated that in general the erosion results were comparable 
those for 45° exit angles. The addition of a tangent section to increase the 
vation of the bucket lip had a beneficial effect for the 37- -1/2° exit angle 
ket. For intermediate discharges, the 37-1/2° exit angle bucket had better 
ysion characteristics than the 45° bucket for the same lip elevation. 

The authors’ conclusion with regards to placement of the invert with re- 
ion to the channel is perhaps a good starting point, but should be used with 
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ution. From Fig. 2 (authors’) it appears that the conclusion was based on 
ving a turbulence chamber downstream of the bucket with a stream bed at a 
gnificant height above the bucket invert. A proper evaluation of the erosion 
itterns by the reader would require some information as to the size of ma- 
rial used in the study and the model scale. With a horizontal bed, in the 
niversity of Alabama tests, erosion to a depth lower than the bucket invert 
vation occurred in the relatively large material. A smaller size material 
.the erosion bed would have shown greater depths of erosion. Model studies 
in only give qualitative indications of erosion and erosive tendencies. The 
able roller formed under the jet may not-cause much erosion in the model, 
it it provides ample opportunity for ball mill action. This could produce a 
sep scour hole in any prototype stream bed material if given sufficient time. 

With regards to the surge height, the authors suggest that this be used as a 
lide in determining the top elevation of training or spray walls. This implies 
iat these walls will be above the surge height. The question is, must they be 
gh? A designer who has no recourse to model studies must arbitrarily, in 
-der to be conservative, answer, yes. The writer does not remember seeing 
1y experimental data that proved that the training walls had to be above the 
irge height. Model tests should be used to determine the required training 
all height with a definite possibility that the walls could be lowered with a 
aving in construction cost. The writer knows of model tests in a government 
iboratory which showed that the training walls for the project could be re- 
iced to one-third of the maximum tailwater depth without adverse erosive 
stion and the project was so constructed. While this was a baffle apron and 
ot a bucket, it is the writer’s conviction that similar results could be justified 
ith the bucket dissipator. 
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SYSTEMATIC CHANGES IN THE BEDS OF ALLUVIAL RIVERS2 


Discussion by T. Blench 


14 BLENCH,! M. ASCE.—This set of observations of dunes in a sand-bed 
iver is a valuable new addition to the collection of quantitative field studies 
f rivers (e.g. Refs. 1-9) that has been growing in the last few years. With 
1€ quantitative facts of canal and river regime now established from field 
bservations, in dynamically satisfying form,(10) the next step in field work 
as to be the collection and analysis of quantitative data on fluctuations about 
egime and on the details of transport of which bed-wave movement is a major 
emonstration; the laboratory, though of specialised value, can hardly provide 
1ore than the bottom point on an adequate correlation of natural phenomena, 
hile theoretical speculations based on rigid boundary occurrences seem un- 
kely to be of more than academic interest till observed facts are established. 

The writer initiated a set of similar observations in 1950 but was not ina 
osition to have them published; a library copy exists.(11) The main object 
f the observations was to provide a set of occurrences to be copied by a 
10del if the model was to be useful for predicting rates of scour. Daily sonic 
oundings were made on a couple of miles of navigation path for a dozen days 
efore and after a large flood of about 400,000 cusecs peak with a tidal addition 
f the order of + 100,000 cusecs. The river bed was sand, and the average 
ed-load charge was of the order of a couple of hundred thousandths by weight. 
t the low discharges of the observations bed waves were small and travelled 
1 groups at some 20 feet a day; at peak flood they had grown to 15 feet high, 
nd 500 feet long, travelling at 250 feet per day; they decayed very much ac- 
ording to the way they grew. Plots of wave number per thousand feet against 
ischarge at three locations in the run gave good and similar correlations. 
he results would hardly be typical of a nontidal river, since the water level 
as controlled closely by proximity to the ocean; with the nontidal river the 
ater surface would vary with discharge and the phenomena would probably be 
SS accentuated. 

In company with the authors the writer hopes their example will interest 
rganizations with the facilities in conducting further observations. He sug- 
2sts that quantitative science would benefit if observations aimed at correlat- 
i bed-wave dimensions and speed with fluid discharge, bed-load constitution 
nedian and dispersion about the median), and bed-load charge. A proper set 
observations would require a full description and record of relevant river 
roperties, viz. sectional form, slope, meander pattern in the zone of slope 
easurement, water levels relative to the sections, mechanical analyses of 
sd material, suspended load observations with load analyses, temperature of 
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water, and duration curve of discharge (i.e. frequency with which various dis: 
charges exceeded on long term); calculation from these data would give an 
approximate idea of bed-load charge, which is not measureable but is highly 
relevant to bed-wave form. 

Generally the writer finds himself in agreement with the opinions in the — 
paper, and able to endorse its facts from his own experiences. 
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SYNTHETIC STORM PATTERN FOR DRAINAGE DESIGN2 


Discussion by M. B. McPherson 


M. B. McPHERSON,! A.M. ASCE.—The Water Department of the City of 
hiladelphia is currently investigating the suitability of the hydrograph method 
r the design of some extensive storm sewer improvements. This investi- 
ition is relatively modest compared with the efforts expended and the results 
tained by Chicago. The authors are to be commended for their clear expo- 
tion of a new approach to a very controversial area of hydrology. The writer 
ids several major reservations insofar as the efficacy of the described syn- 
etic storm is concerned. However, most of this discussion will be restrict- 
i to a modification in the use of the data set forth and a comparison with 
milar modified data for Philadelphia and Baltimore. A comparison with 
altimore data is desirable in view of the extensive overall hydrologic study 
nich is in progress there. ‘1 

With regard to intensity—duration—frequency curves, the U. S. Weather 
ireau data(2) provides a common basis of comparison between the three 
ties. Two—and five-year frequency (or return period) values of intensity 
r sample durations are set forth in Table D-I. The values listed were taken 
om 3x enlargements of the U.S.W.B. graphs, by scaling, and are accurate to 
thin about + 0.05 inches per hour. The data from which the curves were 
tained (“maximum value for each year”) covered the following periods: 


Baltimore - 1903 to 1951, incl. 
Chicago —1905 - 1912 and 1926-1951, incl. 
Philadelphia - 1903 to 1951, incl. 


1e Eltinge and Towne equation was adopted “about 1944” by Chicago accord- 
z to reference 4 of the authors. Independently determined curves were ob- 
ined by Chow(3) with Chicago Weather Bureau data for 1913 to 1947 in- 
usive; the values he obtained on the basis of annual maximums have also 

en listed in Table D-I. For a city without a first-order station, the compa- 
ble approximate values taken from Yarnell’s (4) graphs may be of interest 
i available records to and including 1933, for most of the U. S.). Appraisal 
means of the curves given by Hathaway, (5) based upon the one-hour graphs 
Yarnell, would yield similar results. The Philadelphia and Chicago values 
e quite similar. The Baltimore intensities are somewhat higher. There is 
significant difference between the Weather Bureau, Eltinge and Towne, and 
ow values for Chicago. The equations given by the authors could be used 

r the Philadelphia U.S.W.B. 5-year curve. 

The authors have used 81 station-rainfalls in their Table I, for 27 storms. 
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The intensities for the various durations represent a wide range of frequen- 
cies. In an effort to segregate the data, six storms for Station 14 were select 
ed which have 15-minute intensities closest to the 5-year values: June 20-2 
1937; June 12, 1944; August 9, 1946; July 3, 1954; August 18, 1954; October { 
1954, In Table D-II are given the computed values of “r” for the four du- 
rations of Table I, and for only the 15-, 30-, and 60-minute durations, on the 
basis of both antecedent time and antecedent mass. (The reason for excludir ; 
the 120-minute will be given later.) The selected six-storm group values of 
“r”, for all four durations, are not appreciably different from those computer 
by the authors for all 81 station-rains. The comparative values, with and 
without the 120-minute duration are also reasonably consistent, considering 
the type and extent of the samples. 

As stated previously, the six storms for Station 14 were selected as repr 
senting most closely the 15-minute, 5-year frequency. In Table D-III are 
given the approximate frequencies of all four durations (based on the U.S.W.J 
curves for Chicago). As would be expected, no single Storm has a consister 
frequency for more than two durations. 

Intensity —duration—frequency curves are derived for each duration by tz 
ing maximum excessive rainfalls, regardless of where the critical duration 
occurs in each record, and arranging intensities by rank for statistical analy 
sis. In this way the various durations of a frequency curve represent a seri€é 
of unrelated values from a variety of storms; a frequency curve thus gives 
intensities on the average, which are for a rarer occurrence than indicated. 
The use of frequency-curve data for drainage design, by any method, thus be 
comes highly empirical. The method of transposition of actual storms, usual 
ly employed for larger watersheds, would be preferable. However, two re- 
strictions, or reservations, confound usage of that method: (1) the difficulty 
in arriving at criteria for selecting a “typical” or design storm based on 
either infrequency, pattern or resulting discharge, and (2) even if an actual 
storm is used, point rainfall remains the source of data and there would ben 
particular correlation with the design area in the transposition. If the autho 
cannot adapt the transposition technique to Chicago, with the unique availz 
ity of seventeen private automatic recording rainfall gages, the prospect 2 
pears rather poor for other cities. In any event, the attempt by the authors: 
to divorce themselves as much as possible from the illogical facets of the s¢ 
called “rational” method is indeed noteworthy. 

In Table D-II the antecedence values of “r” for Chicago are similar with 
regard to both time and mass of rainfall. The authors stated that the Chicag 
results “might not be suitable in other localities.” The authors have used 
data which is not publicly available. Also, they did not define the basis for 
selecting the antecedent mass. It is assumed that “excessive rainfall” as 
used follows the U.S.W.B. definition, of a depth in inches equal to or greater 
than the sum of one one-hundredth of the duration in minutes plus twenty- 

‘hundredths of an inch. All U.S.W.B. data published since 1936 have included 
the maximums for given durations but not the actual record of excessive rat 
with respect to time, as was the practice prior to 1936. 

In order that data for the three cities could be equitably compared, the 
records for 1896 to 1935 were searched for each of the three cities, in the 
corresponding issues of the “Annual Report of the Chief of the W.B.” Nine 
storms for each city were selected on the following arbitrary bases: (1) 1 
minute duration precipitations nearest the 5-year U.S.W.B. value, and (2) a 
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eport tabulation including at least 30 minutes of rainfall. Only 15, 30 and 
0-minute durations were considered; there are very few storm records which 
xtend as long as 120-minutes. The writer could see no practical value in 
xtending the record by the “extended duration” procedure. The data for these 
torms are given in Tables D-IV-B, D-IV-C and D-IV-P. The calculated 
alues of the advance proportion “r” for these storms are as follows: 


City By Peak 5-Min. By Antec. Mass 
Baltimore 0.399 Approx. 0.071 
Chicago 0.294 0.546 
Philadelphia 0.494 0.667 


he 0.294 for Chicago based on the peak 5-minutes compares with the corre- 
ponding three duration value of 0.299 in Table D-II, but the 0.546 based on 
1e antecedent mass far exceeds all values given in Table D-II. The figures 
yr Baltimore show the greatest disparity. It was noticed in searching the 
altimore data that large antecedent masses occurred with storms of higher 
‘equency (and lesser intensities). The Baltimore value of 0.071 is approxi- 
vate since a modified Eq. (13) was used. Philadelphia would appear to have 
n “r” of about one-half; the high value of 0.667 was influenced substantially 
y storm P-7. The approximate frequency in years for each duration of each 
torm using U.S.W.B. curves is given in Table D-V. The usual inconsistency 
1 frequency values for different durations, as noted previously in Table D-II, 
3 apparent. 

From the above data there does not appear to be any outright correlation 
etween an “r” calculated on the basis of an antecedent time and on the basis 
f an antecedent mass. In the opinion of the writer, the values 0.386 and 0.376 
f the paper are characteristics of the specific sample used and their near- 
ess in magnitude should be considered fortuitous. In the case of Baltimore, 
lere was no conceivable combination of storms that would have raised the 
antecedent mass value of “r” much above one-tenth. Therefore, the premise 
f the authors is correct: the Chicago values are definitely “not suitable for 
ther localities.” Although the writer has used the method of the authors in 
stermining the value of “r” from mass antecedent rainfall, he can find no 
isis or reason given by the authors to justify the assumption that the ante- 
edent portion of their synthetic hyetograph should have the same trend or 
lape as the part after the peak. 

The authors have advised the writer, by correspondence, that they are de- 
sloping a program for storm sewer runoff analysis by means of a digital 
ymputer. Although evaluation and allocation of the factors involved in the 
yerland flow method of surface routing would become much more complicat- 
1, why could not actual storms be used? Programming would become more 
volved, but the basic program would remain unchanged. Using the same 
orm pattern for both branch sewers and main sewers seems to be somewhat 
consistent. Perhaps design could be based upon the routing of several 
oderately infrequent storm patterns, in which different storms might well 
2 critical for different sewers in the system. It is recognized that a basis 
r selecting the proper storms presents a formidable stumbling block. How- 
rer, the method of storm evaluation proposed by the authors retains too 
any of the fallacies and empiricisms inherent in the “rational” method to 
scommend its principle for general use by others. The pioneering effort of 
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the authors is not to be underrated; perhaps no better solution will be achiey 
for some time to come. 
The value 0.410 on p. 1332-12 should read 0.416. 
The writer was capably assisted in this study by Robert E. Wenzinger, ~ 
Drexel Institute of Technology, Student Chapter ASCE, and recent cooperatiy 
program employee of the Philadelphia Water Department. 


1, 
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TABLE D-I: Comparison of Frequency Data. 
(Tabulated Figures are Intensity in Inches Per Hour). 


DURATION - MINUTES 
City and Data Source 15 30 60 120 
For Five-Year Return Period: 


Baltimore: Weather Bureau 4.70 93.20 1.95 1.25 
Yarnell ho 3.0 1.9 ae 
Chicago: Eltinge & Towne B00 2578) 43477 1.05 
Ven Te Chow 4.0h 2.82 1.81 0.99 
Weather Bureau 4.05 . 2.80 1.80 1.10 
Yarnell 3-9 2.6 ew 1.0 
Philadelphia:Weather Bureau 3-95 2.70 1.70 1.00 
Yarnell 4.2 2.8 1.9 aliege 
For Two-Year Return Period: 
Baltimore: Weather Bureau 3-45 2.30 1.42 0.83 
Yarnell 325 2.4 25 0.8 
Chicago: Ven Te Chow Sele. 2022) (1529 0.71 
Weather Bureau lO T2575 = a .35 0.80 
Yarnell Se Pa eee amily Oc7 
Philadelphia:W eather Bureau 3-05 2.05 1.30 0.75 
Yarnell 3-5 2.3 1.4 0.8 
TABLE D-II: Chicago Data - Paper 1332, Table I, 
Comparison of Values of "r". 
By Peak By Antec. 
5-Min. Mass 
81 Sta.-rains; 15, 30, 60, 120-min. 0.376 0.386 
81 Sta.-rains; 15, 30, 60-min. 0.387 0.350 
Sta. 14, Six storms; 15, 30, 60, 120-min. 0.356 Rosas 
Sta. 14, Six storms; 15, 30, 60-min. 0.299 0.25 


TABLE D-III‘ Chicago Data - Paper 1332,Table I; Six Storms, Station 14 


Equivalent Return Period, or Frequency, Years 


Duration zB ae 

-Min. 30-Min. -Min. -Min. 
Max. ee. Max. Mass Freq. Max. Mass Freq. Max.Mass = 

. 1.56 8.3 1.83 5.0 1.92 i 

ae ce rh 1.7 1.04 1.2 1.43 1.4 

0.92 3.4 1.22 ae 1633 a8 Ey a 

se oO. 1. ° e 

ao An a 9.8 2.49 35. eee as 
0.90 Sot naa ly ¢ 2.6 1.25 LT 1.5 Off 
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TABLE D-IV-B: Nine Approximately 5-year Return Period Storms - Baltimore - 
U.S. WeB. Records 1896 to 1935 Inclusive, 


15-MIN. DURATION 30-MIN. DURATION 60-MIN. DURATION 


Ante- Peak Ante- Peak Ante- Peak 
Max.cedent Five Max.cedent Five | Max. cedent Five 
Mass Mass Min. |Mass Mass Min.|Mass Mass Min. 


Bel : 
7/18/07 1-43 0.41 2 2.20 0.01 ho ne - = 45 mi 
Be2 : 

7/25/01 1.18 Oo 2 1.85 0 2 -- -—- = 35 
B- 1.04 0.04 i. 0.04 -- -- «<< y 
8/e/i9 3 30 3 2T 
B-4 

7/7/05 1.03 0-21 2 1.28 0.14 ieee eh ole va 
Boom ip 

8/27/02 1.02 0.10 2 1.24 0 3 - st ae 25 
B-6 

8/12/34 1.01 0.09 2 1.05 0.09 2 -- -—- = 12 
B-7 

9/1/10 0.98:0.01 2 1.39 0.01 2 -- -—- - 36 
B-8 4 
6/13/15 0.94 0.01 3 1.25 0.01 3 - —— == 38 


B-9 
8/25/05 0.91 0.01 2 1.20 0.01 2 -- -- Ho 


& 
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‘BLE D-IV-C: Nine Approximately 5-Yr. Return Period Storms - Chicago - 
U.S. W.B. Records 1896 to 1935 Inclusive 


15-MIN. DURATION 30-MIN. DURATION 60-MIN. DURATION DURAT. 


Ante- Peak Ante- Peak 
Max. cedent Five | Max. cedent Five Max. cedent Five 
Mass Mass Min. Mass Mass Min. Mass Mass Min. 


C-1 
5/96 1.17 0.18 2 1.30 0.10) 3 -- == oe 17 min. 
c-2 

120/28 1.16 0.23 1 1370". 0,23 2 2.30 0.01 2 60 

C-3 

7/21 1.10 0 1 2.03 0 a -- eee 34 

c-4 

Mm5/0O6 0.97 0.23 2 7526 0-018 i -- -- -- 30 

C-5 

14/09. «0.90 2.00 2 1-20: 40" 1.53 1.79 4 62 

c-6 

3/26 0.87 0.82 3 1255 (0640) 15 -- -- -- ho 

C-7 

12/35 0.87 0.03 1 1.05 0.03 1 1507 £0.03,85-1 2h 

c-38 

18/35 0.86 0.11 2 0.96 0.03 2 -- = = 19 

C-9 


122/24 0.81.0.25 2 1.02 0.25 2 Det O. Ol G 60 
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TABLE D-IV-P: Nine Approximately 5-Yr- Return Period Storms - Philadelphia - 
U.S. W.B. Records 1896 to 1935 Inclusive. 


D -MIN- DURATION 30-MIN. DURATION 60-MIN. DURATION DURAT. 
A Ante- Peak Ante=- Peak n Peak 
©. Max. cedent Five Max. cedent Five Max.cedent Five 
Mass Mass Min. Mass Mass Min. Mass Mass Min. 
P-1 
Tili/1Le weet 0.20" | 2 1.65 0.01 3 “os -- 32 
p-2 
8/16/17 1.09 0.80 3 2.07 0.09 5 we) eee Sew 


P=3 
9/14/0h 1.04 1.29 3 1.63 0.70 6 2.260.30 8 


p-4 

6/26/30 1.02 0.04 1 1-41 0.04 1 2.01 0.04 1 75 
1S) 

8/10/99 0.99 0.40 2 1.48 0.00 4 1.97 0.00 4& ' 55 
P-6 : 

10/4~5/06 0.89 1.42 2 1.44 0.94 mp i — hes a 37 


P=7 
9/15-16/000.87 2.69 3 1.51 2.29 5 2.23 1.57 11 85 
P-8 
9/6/26 0.84 0.01 2 1.16 0.01 2 -- = -- 34 


P=9 
6/19/28 0.83 0.27 3 1.20 0.02 5 - = -- 31 
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TABLE D-V: Equivalent Return Period, or Frequency, Storms Given in 
Tables D-IV-B,-C, & -P. 


DURATION 
15 - minutes 30 = minutes 60 - minutes 
Max.Mass Frequency Max.Mass Frequency Max.Mass Frequency 
Baltimore - 
1.43 nies 2.20 20. -- -- 
1.18 5.2 1.85 9. == == 
1.04 3.4 1.30 2.7 -- -- 
1.03 3.3 1.28 2.6 1.41 2.0 
1.02 3.2 1.24 2.4 -- -- 
Por Saul 1.05 el -- -- 
0.98 2.8 1.39 Sick -- -- 
0.94 2.5 1.25 2.4 -- -- 
0.91 2.3 1.20 2.3 eo Sa ae SS 
Chicago - 
1.17 il. 1.30 3-7 a ae = 
1.16 aes 1.70 13% 2.30 18. 
1.10 8. 2.03 45. -- -- 
0.97 nea 1.26 323 -- -- 
0.90 3.1: 1.20 2.8 1.53 2.8 
0.87 2.7 1.55 8. -- -- 
0.87 2.7 1.05 1.9 1.07 1.3 
0.86 2.6 0.96 1.6 2S =e 
0.81 2.2 1.02 Tee 1.44 2.3 
Philadelphia - 
1.24 18. 1.65 Avs -- -- 
1.09 9. 2.07 60. -- -« 
1.04 417 1.63 ile 2.26 20. 
02 6 aeeet 6. 2.01 10. 
0.99 5.0 1.48 8. 1.97 10. 
0.89 3.2 1.44 Whe -- -- 
0.87 2.9 1.51 9. 2.23 20. 
0.84 2.6 1.16 2.9 -- -- 
0.83 2.5 1.20 3.2 ee oe 
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CHARACTERISTICS OF FLOW OVER TERMINAL WEIRS AND SILLS2 


Discussion by F. V. A. Engel 


F. V. A. ENGEL.!—Both authors have advanced a very interesting thesis, 
e. that under certain circumstances, a sharp-crested weir will not operate 
} such, but behaves like a free overfall. For very large area ratios, exceed- 
g say 0.95, the critical depth which is normally a special feature of Venturi 
* broad-crested weirs, will occur upstream of the sharp-crested weir itself. 
he observed phenomenon deviates considerably from the usual conception of 
scharge characteristics over sharp-crested weirs. The authors do not ap- 
sar to have explained, under what conditions and for what particular reasons 
occurs. There may be several reasons. In the following an attempt will be 
ade to put forward an explanation, and the authors are requested to comment 
i these remarks. This is especially necessary, since the diagrams present- 
| do not permit the arguments advanced to be checked by the writer. It is 
sumed that the authors have complete sets of tables covering their experi- 
ental investigations which would be available to amplify their interesting 
per. 

Fig. 3 shows that a decrease in the discharge coefficient may occur for 
‘ea ratios larger than 0.8, which corresponds to a h/w value of 4. Both the 
shbock and von Mises relationship indicate a discharge coefficient charac- 
ristic asymptotically approaching infinity. As the rate of flow is finite, the 
tual discharge coefficient (excluding the velocity of approach factor) should 
erefore approach zero. Both authorities, however, assume a basic and 
mstant discharge coefficient of 0.605 and 0.611 respectively. The second 
rm in Eq. (8) represents in general the velocity of approach. As previously 
entioned, there are reasons why the discharge coefficient of 0.605 should 
tually decrease rapidly when approaching an area ratio of unity. This may 
plain the deviation in the discharge coefficient characteristic, as shown in 
g. 3, up to a value of 10 for h/w. 

Besides the marked decrease in the discharge coefficient other features 
ay account for the change from the ordinary weir type flow to the “free over- 
1” characteristic. From the description of the experimental procedure one 
ay conclude that, for the very low weir crests, a considerable boundary or 
splacement layer may occur. This may affect the geometry of the sharp- 
ested weir, which projects only slightly over the channel invert so creating 
tificial nozzle entrance conditions, or again a sharp-crested weir may de- 
nerate into a broad-crested one. However, this feature could only be es- 
lished by a very careful investigation of the surface of the water level ap- 
oaching the weir over a considerable length. The profiles given by the 


Proc. Paper 1345, August, 1957, by P. K. Kandaswamy and Hunter Rouse. 
Consultant, 8 Seaton Road, Workington, Cumberland, England. 
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authors in their paper are unfortunately not sufficiently detailed in this re- 
spect. Furthermore, boundary layer effects may explain the considerable dif 
ferences in test results in the two independent investigations by each of the ~ 
authors, as illustrated in Fig. 3. 

The writer2 has recently established another criterion of weir flow assum 
ing a constant basic discharge coefficient, as most authorities do, in separat 
ing the velocity of approach characteristic as a specific factor. This criteri 
was developed with analogous arguments comparing the broad-crested or 
Venturi weir with the sharp-crested weir. In both the former mentioned wei 
types, critical depth ratios related to the area ratio can be established mathe 
matically. For sharp-crested weirs, the ratio of the depth over the weir 
crest to the upstream depth over the weir crest, should lie approximately be 
tween 0.87 and unity for the complete range of area ratios from zero to unity 
This means that the depth over the weir crest is greater than the critical 
depth. ’ 
From the authors’ investigations for very large area ratios, and particul: 
ly when approaching unity, the depth at the crest will be considerably less — 
than the critical depth. This would appear to contradict the criterion ofa ~ 
sharp-crested weir mentioned above. The question arises, at what area rati 
does this change from one condition to the other occur? What is the area 
ratio for which the sharp-crested weir ceases to act as such and conforms t 
the performance of a free overfall (or broad-crested weir) ? 

A diagram presenting the depth ratio as a function of the area ratio 
h/ (h + w) may reveal characteristic features of weir flow dynamics. For © 
instance, some discontinuity may exist. The type of flow may change abrupt 
from ordinary weir flow where the critical section is situated downstream 
the weir crest, to the type of flow with distinct free overfall characteristics. 
This results in the critical depth being reached by the water surface well up 
stream of the weir crest. If a boundary layer exists there may be considera 
ble scale effects and the change from one type of flow to the other may not b 
related to a well defined area ratio. The authors could enhance the value of 
their important results if they would prepare from the test results of their 1 
investigation a diagram based on the outlines mentioned above. ) 


2. Engel, F. V. A. and Stainsby, W., “Velocity-of-approach factors in unifie 
weir equations.” Paper No. 6247 to be published in the February 1958 — 
issue “Journ. Inst. Civil Engineers” (London), 
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100 FREQUENCY CURVES OF NORTH AMERICAN RIVERS2 


Discussion by H. C. Riggs 


H. C. RIGGs,! A.M. ASCE.—Mr: Kuiper’s study of sampling from a log- 
rmal population of 1000 items emphasizes the possible variation of a sample 
om the population in slope and position as well as in shape. Frequency 
irves of annual floods on different streams (after adjustment for differences 
basin characteristics) would not be expected to vary among themselves by 
lance alone as greatly as random samples because of a certain amount of 
spendence associated with their being based on data collected within the last 
) years or so. If chance variation among a group of frequency curves could 
> removed, the adjusted curves would be more representative of the period 
time spanned by the combined records. They would not necessarily be 
oser to their “true” positions. On the other hand they might represent con- 
tions in the next 60 years as well as the “true” curves. 

Many attempts to define a frequency distribution which is generally appli- 
ible to annual flood data have been made. The author has demonstrated that 
e log-normal is an approximation to flood records on 100 American rivers. 
ost of these rivers are large. Drainage basins of large streams usually 
ymprise sub-basins which differ greatly in flood-producing characteristics. 
tegration of these effects and attenuation by channel storage tend to result 
a certain uniformity between flood characteristics of large streams. Hence 
e shapes of frequency curves of large streams might be expected to vary 
ss than the shapes of a group of frequency curves of small streams. Many 
90d-frequency curves of small rivers are curved sharply upward (on paper 
‘signed for the theory of extreme values). It has not been determined 
nether this is a characteristic of small rivers or a function of the relatively 
ort period of record. Few records of 40 years or more are available on 
nall streams. 

Certain streams produce floods in response to more than one cause. Some 
estern streams have snow-melt floods in the spring and rain floods at other 
mes of the year. Similarly, floods near the East Coast may be caused by 
ontal storms or by hurricanes. If the floods resulting from these unlike 
uses actually belong to different distributions then the distribution of annual 
ods is unlike either one. Therefore, it would seem desirable to examine 
e physical and meteorological characteristics of drainage basins and be as- 
red of their reasonable similarity before forcing their frequency curves in- 
statistical agreement. 

The writer agrees with the author’s statement “--- that frequency curves 
maximum annual flood peaks for one particular river station can be drawn 


Proc. Paper 1395, October, 1957, by E. Kuiper. 
-Hydr. Engr., U. S. Geological Survey, Washington 25, D. C. 
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with greater accuracy when due attention is paid to other frequency curves 0) 
the same or comparable streams.” If the individual frequency curves are — 
drawn as straight lines then differences in the slopes of those lines are found 
These differences result from variation in basin characteristics and from 
chance variation in weather. If the variation due to basin characteristics wel 
removed, the adjusted slopes could be averaged to remove the chance wea 
differences. For this approach, selection and measurement of significant 
basin parameters would be a first requirement. 

The author selected 5 parameters and, by a priori reasoning, determined 
the general effect of each. It has been pointed out in this discussion tha 2 
uniformity in drainage basins is the rule rather than the exception. Such non 
uniformity may tend to make the effects of certain parameters compensating 
and to make others difficult to measure. A parameter describing the integra 
ed effect of soil types in a drainage basin of several thousand square miles 
would be so difficult to determine that it would be of little practical use. 

Nevertheless one must begin with a certain concept, test it with actual dat 
and accept, reject, or modify it as required. Differences not readily de- 
scribed quantitatively may be allowed for by grouping streams according to — 
certain conditions. For instance, with the exception of regulated ones, the ~ 
rivers studied by the author might be divided into two groups, (1) those derit 
ing their principal floods from rainfall, and (2) those receiving significant — 
flood volumes from melting snow. 

For use in an example of regional analysis, the writer selected 17 of the 
river-stations according to grouping (1) above. Slopes of the frequency curvy 
for these stations are measured by the author’s R factor. However, relating 
slope of the frequency curve to basin characteristics does not allow slope to 
vary with discharge. This difficulty is avoided by relating the one per cent 
flood to the mean flood and to basin characteristics. Preliminary graphical 
multiple correlation showed a better relation with average discharge of the 
stream than with drainage area so only the former was used. 

The computed regression equation based on data for the 17 stations is ] 


1% flood = 2.57 (mean flood)!-31 (ay, disch.)7-40 


The standard error is -10% and +12%, about half as large as would be obtaine 
by using a mean value of slope for all stations. Data used in the analysis, — 
and computed values of the “one-per-cent flood” are given in Table 1. The 
relation is presented only to demonstrate a method of reducing chance vari- 
ation in slopes of flood-frequency curves. It may not apply to small streams 
Use of additional or more appropriate parameters should give better results 
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PROCEEDINGS PAPERS 


The technical papers published in the past year are identified by number below. Technical- 
division sponsorship is indicated by an abbreviation at the end of each Paper Number, the 
Symbols referring to: Air Transport (AT), City Planning (CP), Construction (CO), Engineering 
Mechanics (EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), Pipeline (PL), 
Power (PO), Sanitary Engineering (SA), Soil Mechanics and Foundations (SM), Structural (ST), 
Surveying and Mapping (SU), and Waterways and Harbors (WW), divisions. Papers sponsored 
by the Board of Direction are identified by the symbols (BD). For titles and order coupons, 
refer to the appropriate issue of “Civil Engineering.” Beginning with Volume 82 (January 
1956) papers were published in Journals of the various Technical Divisions. To locate papers 
in the Journals, the symbols after the paper numbers are followed by a numeral designating 
the issue of a particular Journal in which the paper appeared. For example, Paper 1449 is 
identified as 1449 (HY 6) which indicates that the paper is contained in the sixth issue of the 
Journal of the Hydraulics Division during 1957. 


VOLUME 83 (1957) 


FEBRUARY: 1162(HY1), 1163(HY¥1), 1164(HY1), 1165(HY1), 1166(HY1), 1167(HY1), 1168(SA1), 
11G9(SA1), 1170(SA1), 1171(SA1), 1172(SA1), 1173(SA1), 1174(SA1), 1175(SA1), 1176(SA1), 
See 1178(SA1), 1179(SA1), 1180(SA1), 1181(SA1), 1182(PO1), 1183(PO1), 1184(PO1), 
1185(P01)°. 


MARCH: 1186(ST2), 1187(ST2), 1188(ST2), 1189(ST2), 1190(ST2),1191(ST2), 1192(ST2)c, 1193 
(PL1), 1194(PL1), 1195(PL1). 


APRIL: 1196(EM2), 1197(HY2), 1198(HY2), 1199(HY2), 1200(HY2), 1201(HY2), 1202(HY2), 1203 

(SA2), 1204(SM2), 1205(SM2), 1206(SM2), 1207(SM2), 1208(WW1), 1209(WW1), 1210(WW1), 

x 1211(WW1), 1212(EM2), 1213(EM2), 1214(BM2), 1215(PO2), 1216(PO2), 1217(PO2), 1218 

(SA2), 1219(SA2), 1220(SA2), 1221(SA2), 1222(SA2), 1223(SA2), 1224(SA2), 1225(PO)°, 1226 
’ (WW1)e, 1227(SA2)¢, 1228(SM2)*, 1229(BM2)°, 1230(HY2)°. 


_ MAY: 1231(ST3), 1232(ST3), 1233(ST3), 1234(STS), 1235(IR1), 1236(IR1), 1237(Ww2), 1238 
(WW2), 1239(WW2), 1240(Ww2), 1241(WW2), 1242(WW2), 1243(WW2), 1244(HW2), 1245(HW2), 

_-—s«:1246(8w2), 1247(HW2), 1248(WWw2), 1249(HW2), 1250(HW2), 1251(WW2), 1252(ww2), 1253 
(IR1), 1254(ST3), 1255(ST3), 1256(HW2), 1257(mR1)°, 1258(HWw2)¢, 1259(ST3)¢. 


JUNE: 1260(HY3), 1261(HY3), 1262(HY3), 1263(HY3), 1264(HY3), 1265(HY3), 1266(HY3), 1267 
(POS), 1268(PO3), 1269(SA3), 1270(SA3), 1271(SA3), 1272(SA3), 1273(SA3), 1274(SA3), 1275 
(SA3), 1276(SA3), 1277(HY3), 1278(HY3), 1279(PL2), 1280(PL2), 1281(PL2), 1282(SA3), 1283 
(HY3)¢, 1284(PO3), 1285(PO3), 1286(PO3), 1287(PO3)¢, 1288(SA3)°. 


_ JULY: 1289(SM3), 1290(EM3), 1291(EM3), 1292(EM3), 1293(EM3), 1294(HW3), 1295(HW3), 1296 
(HW3), 1297(HW3), 1288(HW3), 1299(SM3), 1300(SM3), 1301(SM3), 1302(ST4), 1303(ST4), 1304 

_  (ST4), 1305(SU1), 1306(SU1), 1307(SU1), 1308(ST4), 1309(SM3), 1310(SU1)¢, 1311(EM3)¢, 1312 
2 (ST4), 1313(ST4), 1314(ST4), 1315(ST4), 1316(ST4), 1317(ST4), 1318(ST4), 1319(SM3)¢, 1320 
(ST4), 1321(ST4), 1322(EM3), 1323(AT1), 1324(AT1), 1325(AT1), 1326(AT1), 1327(AT1), 1328 
. (AT1)¢, 1329(ST4)¢. 


_ AUGUST: 1330(HY4), 1331(HY4), 1332(HY4), 1333(SA4), 1334(SA4), 1335(SA4), 1336(SA4), 1337 
7 (SA4), 1338(SA4), 1339(CO1), 1340(€01), 1341(CO1), 1342(CO1), 1343(CO1), 1844(PO4), 1345 
(HY4), 1346(PO4)¢, 1347(BD1), 1348(HY4)°, 1349(SA4)°, 1850(PO4), 1351(P04). 


SEPTEMBER: 1352(IR2), 1353(ST5), 1354(ST5), 1355(ST5), 1356(ST5), 1357(ST5), 1358(ST5), 
1359(IR2), 1360(IR2), 1361(ST5), 1362(IR2), 1363(IR2), 1364(IR2), 1365(WW3), 1366(WW3), 
i367(WWS), 1368(WW3), 1369(WW3), 1370(WW3), 1971(HW4), 1372(HW4), 1373(HW4), 1374 
(HW4), 1375(PL3), 1376(PL3), 1377(IR2)¢, 1878(HW4)°, 1379(IR2), 1380(HW4), 1381(wws)c, 
1382(ST5)©, 1383(PL3)°, 1384(IR2), 1385(HW4), 1386(HW4). 


OCTOBER: 1387(CP2), 1388(CP2), 1389(BM4), 1390(EM4), 1391(HY5), 1892(HY5), 1393(HY5), 
1394(HY5), 1395(HY5), 1396(PO5), 1397(PO5), 1398(PO5), 1399(EM4), 1400(SA5), 1401(HY5), 
1402(HY5), 1403(HY5), 1404(HY¥5), 1405(HY5), 1406(HY5), 1407(SA5), 1408(SA5), 1409(SA5), 
. 1410(SA5), 1411(SA5), 1412(EM4), 1413(EM4), 1414(PO5), 1415(EM4)°, 1416(PO5)¢, 1417 
s (HY¥5)°, 1418(EM4), 1419(PO5), 1420(P05), 1421(PO5), 1422(SA5)°, 1423(SA5), 1424(EM4), 
ei 1425(CP2). 


& c 

_ NOVEMBER: 1426(SM4), 1427(SM4), 1428(SM4), 1429(SM4), 1430(SM4)°, 1431(ST6), 1432(ST6), 
ie: 1433(ST6), ner 1435(ST6), 1436(ST6), 1437(STS6), 1438(SM4), 1439(SM4), sae) 
A 1441(ST6), 1442(ST6)°, 1443(SU2), 1444(SU2), 1445(SU2), 1446(SU2), 1447(SU2), 1448(8U2)°, 


MBER: 1449(HY6), 1450(HY6), 1451(HY6), 1452(HY6), 1453(HY6), 1454(HY6), 1455(HY6), 
4 yassuayer®, 1457(PO6), 1458(PO6), 1459(PO6), 1460(PO6)©, 1461(SAG), 1462(SAG), 1463(SA8), 
"_-1484(SA6), 1465(SAG), 1466(SAG)¢, 1467(AT2), 1468(AT2), 1469(AT2), 1470(AT2), 1471(AT2), 
. 1472(AT2), 1473(AT2), 1474(AT2), 1475(AT2), 1476(AT2), 1477(AT2), 1478(AT2), 1479(AT2), 
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1480(AT2), 1481(AT2), 1482(AT2), 1483(AT2), 1484(AT2), 1485(AT2)¢, 1486(BD2),1487(BD2), 
1488(PO6), 1489(PO6), 1490(BD2), 1491(BD2), 1492(H¥6), 1493(BD2). 


VOLUME 84 (1958) 

* Y: 1494(EM1), 1495(EM1), 1496(EM1), 1497(IR1), 149B8(IR1), 1499(IR1), 1500(IR1), 

ater psa eeeanay 1504(IR1), 1505((R1), 1506(IR1), 1507(IR1), 1508(ST1), 
1509(ST1), 1510(ST1), 1511(ST1), 1512(ST1), 1513(WW1), 1514(WW1), LESSEN WH), zeletwls: 
1517(WW1), 1518(WWi), 1519(ST1), A ler 1521(IR1)°, 1522(ST1)°, 1523(WW1)°, 
1524(HW1), 1525(HW1), 1526(HW1)°, 1527 5 

FEBRUARY: 1528(HY1), 1529(PO1), 1530(HY1), 1531(HY1), 1532(HY1), 1633(8A1), 1534(SA1), 

1535(SM1), 1536(SM1), 1537(SM1), 1538(PO1)°, 1539(SA1), 1540(SA1), 1541(SA1), 1542(SA1), 

1543(SA1), 1544(SM1), 1545(SM1), 1546(SM1), 1547(SM1), 1548(SM1), rata iiects AESOANSL 

1551(SM1), 1552(SM1), 1553(PO1), 1554(PO1), 1555(PO1), 1556(PO1), 1557(SA1)¢, 1558(HY1)°, 

1559(SM1)°. 

Ssion of several papers, grouped by divisions. 
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